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PART  NUMBERS:  12585710-275,  Compound  Actuator  Assembly 
DESCRIPTION:  RECOIL  SYSTEM 

The  recoil  system  consists  of  two  recoil  cylinders,  a  counter¬ 
recoil  cylinder,  an  energy  recovery  cylinder,  and  a  load  position 
actuator-  There  are  also  two  counter-recoil  accumulators  and  a 
reservoir  accumulator. 

The  two  recoil  cylinders  and  the  energy  storage  cylinder  have  a  6 
inch  free  recoi 1  length  to  keep  retardation  forces  to  a  minimum 
until  the  project  exits  the  barrel . 

The  recoi 1  system  is  designed  to  give  a  constant  78,000  lbs. 
retaring  force  profile  when  fired  from  load  position  with  a 
nitrogen  volume  of  4800  cu.  in.,  pressure  of  2166  psi ,  and 
temperature  of  —25  F. 

In  normal  operation  (firing  from  battery  position)  the  recoil 
force  profile  consists  of  6  inches  of  free  recoil,  then  hard 
recoil  of  75,000  lbs  until  recoi 1  travel  is  3.25  feet.  Finally, 
the  recoil  force  is  gradually  reduced  as  recoil  travel  approaches 
full  recoil  stroke.  This  trailing  off  of  retardation  force 
during  the  last  5.5  feet  of  recoil  travel  is  to  improve  firing 
stability. 

STATUS: 

All  dimensions  and  mounting  requirements  have  been  determined. 
The  final  orifice  profile  was  pending  completion  of  weight,  CG 
and  inertia  analysis  to  provide  optimum  firing  stability. 
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VC-1147. 

UM-16.675 

DB-6.15 

WP-96. 

WC-28. 

WA-4 . 

PRINT  5 

5  F0RMAT(//,1X, 'MUZZLE  VELOCITY( FT/SEC)-?'  ) 

READ (  6  ,  *  )  VM 
VM-2710 . 

RHOC-.06 
RTB-11563400 . 

GAMMA-1.238 

CK-1./7. 

RL-8.75 

DRP-3. 

DRR-2 . 

DCRP-3 . 

DCRR-2.25 
P0-3000 . 

V0-7000 . 

CK2-1.3 
WR-3800 . 

BETA- . 7 
BULK-150000. 

RHOOIL— .0308 
G-32.174 
PI-3.141592 
AB-PI *DB*DB/4 . 

D— VC+AB*UM*12 . 

AR-PI * ( DRP*DRP-DRR*DRR ) /4 . 

ACR-PI*(DCRP*DCRP-DCRR*DCRR)/4. 

GAMAM- GAMMA- 1 . 

GAMAP— GAMMA+1 . 

GAMAl  -GAMAP/GAMAM 
GAMA2-2 . *  GAMMA/ ( 1 . -GAMMA) 

RAT I 0-1 . +WC/( 6 . *WP ) 

RTO-RTB-GAMAM* ( 1 . /6 . + ( 1 . +CK ) *WP/ ( 2 . *WC ) ) * VM* VM 

SIG-D*SQRT( ( ( GAMAP/2 . ) *  * GAMAl )/( GAMMA*RTO* RATIO ) )/( AB*GAMAM*6 . ) 
PM-WC*RTO*RATIO*12 ./( G*D ) 

C  WRITE  ( *, ' ( 5F13.4 ) ' )  PM, RATIO, RTO, SIG,D 

WEFF-WP+ ( WC+WA )/2 . 

A— VM/ ( 1 . - ( PM*AB*G*UM )/( VM*VM*WEFF ) ) 

B— ( ( A/VM ) — 1 . ) *UM 
TM-2 . *UM/VM 

UO-EXP ( LOG ( UM ) -2 . -UM/B ) 

XR-0 . 

XRD-0. 

T-0. 

U-UO 
DT- .001 
XRK-0 . 

XRDL-0 . 

50  Cl  — 1. 

C2-1 . 

C3--1 . 

DELXR-0 . 

DELXRD-0 . 

DO  300  1-1,4 
Cl-Cl+2 . *C2 
C2-C2-.5 
Dl-1 . 5-Cl* . 5+C3 


•  L 


C3-0. 

XREFF-XR+Dl *XRK 
XRDErr-XRD+Dl *XRDL 
XRK-DT*XRDEFF 
IF(T.GT.TM)  GO  TO  100 
DELT-Dl*DT/2. 

VP-A*U/{B+U) 

BB-a+U- ( VP+A ) *DELT 

CC— (  ( VP+A)  *U+VP*B  )  *DELT 

DELU— ( -BB+SQRT ( BB*BB-4 . *CC } )/2. 

DELT2- ( B*  LOG ( 1 . +DELU/U ) +DELU ) /A 

UEFF-U+DELU+ ( D1 *DT-DELT2 ) *A* ( U+DELU )/( B+U+DELU ) 

FGAS— (WEFF/G) *A*A*B*UEFF/( B+UEFF )  **3 . 

GO  TO  200 

100  TEFF-T-TM+D1*DT 

FGAS-PM*AB* ( 1 . -BETA) *  ( 1 . +TEFF/SIG ) *  *GAMA2 
200  I F ( XREFF . LT . 1 . )  GO  TO  250 
IF(XREFF.LT. 1.5)  GO  TO  260 
I F ( XREFF . LT .3.25)  GO  TO  270 
I F ( XREFF . LT .4.5)  GO  TO  280 
I F ( XREFF . LT .8.5)  GO  TO  290 
GO  TO  500 

250  FR-40000 . *XREFF*XREFF+5000 . 

GO  TO  295 

260  FR-(-160.*(XREFF**3. ) +520 . *XREFF*XREFF-480 . *XREFF+165 . ) *1000 . 

GO  TO  295 
270  FR-75000. 

GO  TO  295 

280  FR* ( 3 . 84* ( XREFF* *3 . )-48 . 64*XREFF**2 .+194 . 48*XREFF-175 . 12 ) *1000 . 
GO  TO  295 

290  FR— 10000. *XREFF+110000. 

295  XRDL-DT*G* ( FGAS-FR ) /WR 
DELXR-DELXR+C1 *XRK 
DELXRD-DELXRD+C1 *XRDL 
300  CONTINUE 

XR-XR+DELXR/6 . 

XRD-XRD+DELXRD/6 . 

U-UEFF 

t-t+dt 

I F ( XR . LT . 1 . )  GO  TO  350 
I F ( XR . LT . 1 . 5 )  GO  TO  360 
I F ( XR . LT .3.25)  GO  TO  370 
IF ( XR . LT .4.5)  GO  TO  380 
I F ( XR . LT. 8.5)  GO  TO  390 
GO  TO  500 

350  FR-40000.*XR*XR+5000. 

DFDX-80000 . *XR 
GO  TO  400 

360  FR-(-160.*(XR**3. )+520 . *XR*XR-480 . *XR+165 . )*1000. 

DFDX« ( -480 . *XR*XR+1040 . *XR-480 . )*1000. 

GO  TO  400 
370  FR-75000. 

DFDX-0. 

GO  TO  400 

380  FR- ( 3 . 84* ( XR* *3 . ) -48 . 64*XR*XR+194 . 48*XR-175 . 12 ) *1000 . 

DFDX- (11. 52*XR*XR-97 . 28*XR+194 . 48 ) *1000. 

GO  TO  400 

390  FR-( -10000 . *XR+110000  .  ) 

DFDX— 10000. 

400  ER-1. 

IF ( XR . LT . . 5 )  ER-0 . 

DV-ACR* ( XR+ER* ( XR- . 5 ) )*12. 


o  n 


FCR-PO*ACR*( 1 .+ER)/( 1 .-DV/V0)**CK2 
DFCRDX-FCR*CK2*12 . *ACR* { 1 . +  ER )/( VO-DV ) 

RECOIL- ( FR- FCR ) /2 . 

I F ( RECOIL . LE . 0 . )  RECOIL-1. 

P-RECOIL/AR 
DFRDX-DFDX-DFCRDX 
DPDX-DFRDX/ ( 2 . *AR ) 

DPDT-DPDX*XRD 
V— AR* ( RL-XR ) *12 . 

AO-(-DPDT*V/BULK+AR*XRD*12. )/( 1 1 5 . 5*SQRT( P ) ) 

I F ( AO . GT .1.75)  AO-1.75 
IF( AO. LT. 0 . )  AO-0. 

PRINT  450 , T, XR , AO , XRD , P , FR 

450  FORMAT ( 3X , F4 . 3 , 4X, F5 . 3 , 4X , F9 . 7 , 4X , F5 . 2 , 4X , F6 . 0 , 4X, F7 . 0 ) 
WRITE  (*, » (1X.F7.0) ' )  FR 
WRITE  (*, ' ( IX , F5 . 3 ) ' )  XR 


I F ( XRD . LT . 0 . )  GO  TO  500 
I F ( XR . GT .8.75)  GO  TO  500 
GO  TO  50 
500  STOP 
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VC-1147. 

UN-16 .675 

OB-6.15 

WP-96. 

WC-28. 

WA-4. 

Vh-2710 . 

RHOC-.06 
RTB-1 156 3400 . 

GANNA- 1.2 38 
CK-1./7. 

RL-8.75 
DRP-3 . 

DRR-2 . 

DCRP-3 . 

DCRR-2.25 

P0-3000. 

V0-7000 . 

CK2-1 . 3 
WR-3800  . 

BETA- . 7 
BULK-150000. 

RHOOIL- .0308 
G-32.174 
PI-3.141592 
AB— PI*DB*DB/4 . 

D— VC+AB  *  UN  *12. 

AR-PI * ( DRP*DRP-DRR*DRR )/4 . 

ACR-PI * ( DCRP*DCRP-DCRR*DCRR )/4 . 

GAMAM-GAMMA-1 . 

GANAP-GANNA+1. 

GAHA1-GAHAP/GAHAK 
GANA2-2 . * GANNA/ ( 1 .-GANNA) 

RATIO-1 . +WC/( 6 . *WP ) 

RTO-RTQ-GANAN* ( 1 . /6 . + ( 1 . +CK ) *WP/( 2 . *WC ) ) *VH*VN 

SIG-D*SQRT( ( ( GANAP/2 . ) **GANAl )/(GANNA*RTO*RATIO) )/( AB*GANAN*6 . ) 
PM-WC*RTO*RATIO*12 . /( G*D ) 

C  WRITE  ( *  ,  M  5F1 3 . 4 ) ' )  PN , RATIO, RTO, SIG ,  D 

weff«wp+ ( wc+wa ) /2 . 

A— VN/( 1 .-( PN*AB*G*UH)/( VN*VN*WEFF) ) 

B- ( ( A/VN ) - 1 . ) *UN 
TH-2 . *UH/VN 

UO-EXP ( LOG ( UN ) -2 . -UN/B ) 

XR-0 . 

XRD-0. 

T-0  . 

U-UO 

DT-.001 

XRK-0. 

XRDL-0. 

50  Cl  — 1. 

C2-1 . 

C3--1 . 

DELXR-0. 

DELXRD-0 . 

DO  300  1-1,4 
Cl-Cl+2 . *C2 
C2-C2-.5 
Dl-1 . 5-C1* . 5+C3 
C3-0 . 

XREFF-XR+Dl *XRK 
XRDEFF-XRD+Dl *XRDL 


XRK-DT*XRDEFF 
IF(T.GT.TM)  GO  TO  100 
DELT-Dl*DT/2. 

VP-A*U/(B+U) 

BB-B+U- ( VP+A ) *DELT 

CC— (  ( VP+A)  *U+VP*B )  *DELT 

DELU- ( -BB+SQRT( BB*BB-4 .  *CC ) )/2. 

DELT2- ( B*LOG ( 1 . +DELU/U ) +DELU ) /A 

UEFF-U+DELU+ ( Dl*DT-DELT2 ) *A* ( U+DELU ) / ( B+U+DELU ) 

FGAS- ( WEFF/G ) *A*A*B*UEFF/{ B+UEFF) **3 . 

GO  TO  200 

100  TEFF-T-TM+Dl*DT 

FGAS— PM* AB* ( 1 . —BETA ) * ( 1 . +TEFF/SIG ) **GAMA2 
200  I F ( XREFF . LT . 1 . )  GO  TO  250 
I F ( XREFF . LT . 1.5)  GO  TO  260 
IF ( XREFF . LT .3.25}  GO  TO  270 
I F( XREFF.LT. 4.5)  GO  TO  280 
IF( XREFF.LT. 8 .5)  GO  TO  290 
GO  TO  500 

250  FR-40000.*XREFF*XREFF+5000. 

GO  TO  295 

260  FR- ( — 160 . * ( XREFF** 3 . ) +520 . *XREFF*XREFF-480 . *XREFF+165 . )*1000. 

GO  TO  295 
270  FR-75000. 

GO  TO  295 

280  FR-( 3. 84* (XREFF** 3. ) -48 . 64*XREFF**2 . +194 . 48*XREFF-175 . 12 ) *1000 . 
GO  TO  295 

290  TR— 10000. *XREFF+110000. 

295  XRDL“DT*G* ( FGAS-FR )/WR 
DELXR-DELXR+C1 *XRK 
DELXRD-DELXRD+C1 *XRDL 
300  CONTINUE 

XR-XR+DELXR/6 . 

XRD-XRD+DELXRD/6 . 

U-UEFF 

T-T+DT 

I F ( XR . LT . 1 . )  GO  TO  350 
I F ( XR. LT .1.5)  GO  TO  360 
I F ( XR . LT. 3.25)  GO  TO  370 
I F ( XR . LT .4.5)  GO  TO  380 
IF ( XR . LT .8.5)  GO  TO  390 
GO  TO  500 

350  FR-40000 . *XR*XR+5000 . 

DFDX-80000. *XR 


GO  TO  400 

FR-(-160.*(XR**3. )+520 . *XR*XR-480 . *XR+165 . )*1000. 
DFDX- ( -480 . *XR*XR+1040 . *XR-480 . )*1000. 

GO  TO  400 
FR-75000. 

DFDX-0 . 

GO  TO  400 

FR— ( 3 . 84 * ( XR* * 3  . )-48 .64*XR*XR+194 . 48*XR-175 . 12 ) *1000  , 
DFDX-( 11 . 52*XR*XR-97 . 28*XR+194 . 48  ) *1000  . 

GO  TO  400 

FR- (-1 0000. *XR+1 10000. ) 

DFDX— 10000  . 

ER-1 . 

I F ( XR . LT . . 5 )  ER-0 . 

DV-ACR*(XR+ER*(XR-.5) ) *12  . 

FCR-P0*ACR* { 1 . +ER )/( 1 . -DV/V0 ) **CK2 
DFCRDX-FCR*CK2*1 2 . *ACR* ( 1 . +ER )/( VO-DV ) 

RECOIL— ( FR-FCR)/2 . 


o  n 


IF { RECOIL. LE. 0. )  RECOIL-1. 

P-RECOIL/AR 
DFRDX-DFDX-DFCRDX 
DPDX— DFRDX/( 2 . *AR) 

DPDT-DPDX*XRD 
V-AR* ( RL-XR) *12 . 

AO— ( — DPDT*V/BULK+AR*XRD*12 . )/( 115 . 5*SQRT( P ) ) 
I F ( AO . GT .1.75)  AO-1.75 
I F ( AO . LT . 0 . )  AO-O. 

IF ( XR . LT . 2 . 9 )  GO  TO  475 
PRINT  450, XR 
450  FORMAT ( IX , F5 . 3 ) 

PRINT  460, AO 
460  FORMAT ( IX, F5. 4) 

WRITE  (*, ' (1X,F7.0) ' )  FR 
WRITE  (*, ' (1X,F5.3) ' )  XR 
475  IF( XRD. LT. 0 . )  GO  TO  500 
IF(XR.GT. 8.75)  GO  TO  500 
GO  TO  50 
500  STOP 
END 
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DESCRIPTION: 


RECOIL  SYSTEM  (WITH 
ABANDONED  DECEMBER, 


FIXED  ORIFICE, 
1986) 


Design  analysis  of  this  system  over  Christmas  break  in  1986 
revealed  extensive  complications  due  to  the  combination  of 
temperature  f 1 uctuat i ons ,  energy  recovery  levels,  and  possible 
cookoff  at  latch  (load)  position  CB/700  pg  143.  Review  of 
cylinders  and  rods  designed  and  detailed  indicated  changing  to 
conventional  orifice  rods  was  a  superior  path. 


STATUS: 

Abandoned  in  December  1986. 


AUTHOR:  Jeff  Ireland 
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3395. 

5570. 

6559. 

3403. 

b  5  3  5 . 

6535  . 

3411. 

b  5  00  . 

6499  . 

3419  . 

b4b  3 . 

6462. 

3427. 

b  4  2  5  . 

6425. 

3435  . 

b387. 

6387. 

3443. 

6348  . 

6346. 

345  1. 

6309. 

6309  . 

3459  . 

6271  . 

6270  . 

3467. 

5232. 

6232  . 

3474  . 

6193. 

6193. 

3482  . 

6154. 

6154. 

3490. 

b  1 1 6  . 

6116. 

349e  . 

fa  0  7  8  . 

6078  . 

3505. 

5041  . 

6041  . 

3513. 

6003. 

6003. 

3520  . 

5957. 

5957. 

3526  . 

5931. 

5931. 

3535  . 

5895. 

5895. 

3543  . 

5850  . 

5050  . 

3550. 

5825  . 

5625  . 

3558  . 

5791. 

5791  . 

3555. 
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.0590 

16.71 

3  .  A  7  A 

.  06  0  0 

16.71 

3.527 

.0610 

16.71 

3.530 

.0620 

16.71 

3.533 

.0636 

16.71 

3.685 

.  06A  0 

It.  71 

3.737 

.  0656 

16.71 

3.7a9 

.0660 

16.71 

3. 3  A0 

.  0670 

It.  71 

3.8  90 

.  068  0 

16.71 

3. 9 A0 

.  0690 

16.71 

3.969 

.  C700 

16.71 

8.039 

.0710 

16.71 

A  .027 

.  0720 

16.71 

8.135 

.073° 

16.71 

8.133 

.  C7ao 

16.71 

8.230 

.0750 

lb.  71 

A  .  277 

.0760 

16. 71 

9.323 

.  0770 

lb.  71 

8.359 

.  078  6 

16.71 

A  .  A  1  5 

.0790 

15.71 

A  .  A  5  0 

.  0800 

16.71 

8. 50  A 

.0816 

lo.  71 

8 , 5  A  3 

.  0920 

lb. 71 

8 . 5  ->2 

.083  0 

16.71 

8.535 

.  0  8  A  o 

lb. 71 

8.678 

.0*50 

16.71 

8.721 

.  0866 

16.71 

8.753 

.  0970 

16.71 

8.30A 

.0886 

16.71 

8.8A5 

.  0890 

16.71 

8.336 

.090  6 

16.71 

8.927 

.09in 

16.71 

8.957 

.0920 

16.71 

5.005 

.  0 930 

16.71 

5. 0A5 

•  09  A  0 

16.71 

5 . 0  8  A 

.095  0 

16.71 

5.122 

•  09  b  0 

16.71 

5.150 

.  0970 

16.71 

5.193 

.0960 

16.7  1 

5.235 

.0«96 

16.71 

5.271 

.  1006 

16.71 

5  .303 

.1010 

16.71 

5 . 3  A  3 

.1020 

16.71 

5.3  79 

.1036 

16.71 

5.  A1  A 

.  10A0 

16.71 

5.  A  A9 

.  1050 

16.71 

5  .  A  3  3 

.1060 

16.7  1 

5.517 

.1070 

16.71 

5.551 

.1060 

16.71 

5 . 5  5  A 

.  1090 

16.71 

5.516 

.  11C0 

16.71 

5. 5  A9 

.1110 

16.71 

5.531 

.1120 

16.71 

5.712 

.1130 

16.71 

5 . 7  A  A 

.  1 1A0 

16.71 

5.77a 

.1150 

16.71 

5.305 

.1166 

16.71 

5.835 

.1170 

16.71 

5.855 

.  neo 

16.71 

5 . 8  9  A 

53.9 

5757. 

5757 

53.5 

5728  . 

5728 

5  3.0 

5592  . 

5592 

52.5 

5650. 

5650 

52.0 

5529. 

5629 

5  1.5 

5596. 

5598 

51.1 

55b7 . 

5567 

50.6 

5538  . 

5538 

50.2 

5508  . 

5503 

A  9.  7 

5880  . 

588  0  * 

89.3 

5851. 

585  1, 

A  8. 8 

5828. 

5a2  8  , 

A  8  .  A 

5397. 

5397, 

A  7 . 9 

5370. 

53  70, 

A7 .5 

5383. 

5383, 

A  7 . 0 

5318. 

5313, 

A  6 . 5 

5292  . 

5292  , 

A  5  .  1 

5267. 

525  7. 

A  5  .  7 

5283. 

5283. 

85.2 

5219. 

5219. 

88.3 

5  195. 

5195. 

88.  A 

5172. 

5172  . 

A3. 9 

5189. 

518  9, 

A3. 5 

5127. 

5127. 

83.1 

6105. 

5105. 

82.7 

5088. 

50  B  A  . 

82.2 

5062. 

5052  . 

81.8 

5082  . 

50A2. 

81.  A 

5021  . 

5021  . 

81.0 

5001  . 

5001  . 

80.6 

8931  . 

8  9  81. 

80.1 

8962  . 

8952. 

39.7 

8  9  A  3  . 

898j. 

39.3 

8  92  A  . 

8928  . 

36.9 

A  9  06  . 

8905. 

36.5 

A  8  6  9  . 

8886. 

38.1 

A 8  70  . 

8870. 

37.7 

8853. 

8853. 

37.3 

86  35. 

8835  . 

35.9 

8319. 

8619  . 

35.5 

8802  . 

8802  . 

35.  1 

8766. 

8765. 

35.7 

A  7  70 . 

8770. 

35.3 

A  7  5  A  . 

8  758. 

38.9 

8739  . 

8739  . 

38.5 

8728  . 

8728  . 

38.  1 

8709. 

8709  . 

33.7 

8698  . 

869A. 

33.3 

8660. 

8oB0. 

32.9 

8  6  66. 

8666. 

32  .  b 

8  d  5  2  . 

8652. 

32.2 

8636  . 

8638  . 

3  1.3 

8625. 

8625. 

31.8 

8612. 

8612. 

31.0 

8599. 

8599. 

30.6 

8566. 

85  8  5  . 

30.3 

8578  . 

8578  . 

29.9 

8562. 

8562  . 

29.5 

8550. 

85  50. 

29.1 

853b. 

8536  . 

ft 


3572. 

3579. 

3567. 

359A  . 

3501  . 

3503  . 

3615. 

3622. 

3629  . 

3536. 

36  A  2  . 

36A  9  . 

3656 . 

366  2  . 

3669. 

3676. 

35S2  . 

3589. 

3695. 

3702  . 

3706. 

371  A  . 

3720. 

3727. 

3733  . 

3739. 

37A5. 

3751  . 

3757. 

3  7  b  3  . 

3769. 

3775. 

3760. 

3786. 

3792  . 

3798  . 

3803. 

3809. 

3  8 1  A  , 
3820. 
3825. 
3830  . 
3636. 
38A1. 

3  8  A  b  . 
3851  . 
3857  . 
3862  . 
3857. 
3e72. 
3877  . 
3881  . 
3886. 
3e*l  . 
3896. 
3900. 
3905. 
3910. 
391  A  . 
3919. 
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«  1  1  9 

It. 71 

5.923 

.  12C  0 

lb.  71 

5.951 

.121° 

lb.  71 

5.930 

.  1?2'> 

lb.  71 

b  .  0  0  7 

.  1230 

It,. 71 

b.035 

.124  0 

16.71 

b.  0b2 

.  1  25  0 

16.71 

b  .  069 

.1260 

16.71 

6.115 

.12  7  0 

16.71 

6.  141 

.  126^ 

lb.  71 

b.  167 

.1290 

16.71 

b.  142 

.1300 

lb.  71 

b.  217 

.1310 

16.71 

b.  241 

.  1  ?  2  r 

It. 71 

b.  2d5 

.13  3  0 

16.71 

b.  269 

.1340 

16.7  1 

b.  313 

.135  0 

16.71 

6.336 

.13  6  0 

Id.  71 

6.358 

.1370 

lo.  71 

b.  361 

.136  0 

lb. 71 

b.  403 

.1390 

lb.  71 

6.424 

.14  0  0 

lb. 71 

6.44b 

.  1  4  1  C 

16.71 

6.467 

.  1420 

16.71 

6.487 

.1430 

16.71 

6.507 

.1440 

16.71 

6.527 

.  145C 

16.71 

6.547 

.1460 

16.71 

b  .  5  5  6 

.1470 

lb.  71 

6 . 5  b  5 

.146  0 

16.71 

6.603 

.14  4'? 

Id. 71 

6.621 

.  15C 

16.71 

6.639 

.1510 

16.71 

6.657 

.1520 

16.71 

6.674 

.15  3  0 

16.71 

6. 6  90 

.15  4  0 

16.71 

6.707 

.1550 

Id.  71 

6.723 

.15  6  0 

16.71 

6.738 

.1570 

16.71 

6.754 

.156  0 

16.  71 

b  .  7  d  9 

.15  9  0 

16.71 

6.  76  3 

.  16C0 

16.71 

6. 798 

.1610 

16.71 

6.512 

.  162  o 

16.71 

6.925 

.  1630 

Id. 71 

6.538 

.  1640 

Id.  71 

6.951 

.  1650 

16.71 

b  .  86  4 

.166  0 

It. 71 

6.976 

.1670 

16.71 

6.338 

.  16en 

16.71 

6.  900 

.16  9  0 

16.71 

6.911 

.  17C0 

16.71 

6. 922 

.1710 

lb. 71 

6.932 

.1720 

lb.  71 

6.942 

.1730 

lb. 71 

b.  952 

.1740 

lb. 71 

6.962 

.  1750 

lb. 71 

6.971 

.1760 

lb.  71 

6.960 

.1770 

lb.  71 

6.936 

.1760 

lb.  71 

6.996 

28. B 

4527. 

4527 

28.4 

4515. 

4515 

28.0 

4504. 

4504 

27.6 

4493. 

4493 

27.3 

4482. 

4482 

26.9 

4472. 

4472. 

2b.  5 

44b2. 

4462. 

26.2 

4451. 

4451. 

25.8 

4442  . 

4442. 

25.4 

4432. 

4432. 

25.1 

4422  . 

4422. 

24.7 

4413. 

4413. 

24.3 

4  4  (J4  . 

44  04  . 

24.0 

4395. 

4395  . 

23.6 

4386. 

433b. 

23.3 

4377. 

4377. 

22.9 

4368. 

4366. 

22.5 

4360. 

4360. 

22.2 

4352  . 

4352. 

21.8 

4344  . 

4344  . 

21.5 

4336. 

4336. 

21.1 

4326. 

4326. 

20.9 

4321  . 

4321. 

20.4 

4313. 

4313. 

20.0 

4306. 

4306. 

19.7 

4299 . 

4299  . 

19.3 

4292  . 

4292  . 

19.  0 

4285  . 

4235  . 

18.6 

4279. 

4579. 

18.3 

4272  . 

42  72  . 

17.9 

4  2  6  b  . 

4266. 

17.6 

4260. 

4260. 

17.2 

4254. 

4254. 

16.9 

4  2  45  . 

4246  . 

16.5 

4242. 

4242  . 

16.2 

4236. 

4236. 

15.8 

4231. 

423  1. 

15.5 

4225. 

4225. 

15.2 

4220. 

4220. 

14.8 

4215. 

4215. 

14.5 

4210. 

4210. 

14.1 

4205. 

4205. 

13.8 

4201  . 

4201. 

13,4 

4196. 

419  b. 

13.1 

4191  . 

4191  . 

12.7 

4  1  o  7  . 

4187. 

12.4 

4183. 

4183. 

12.1 

4179. 

4179. 

11.7 

4175. 

4175. 

11.4 

4171  . 

4171. 

11.0 

4167. 

4157. 

10.7 

4  lb  4  . 

4164. 

10.4 

4160. 

4160. 

10.0 

4157. 

4157. 

9.7 

4154  . 

4154. 

9.3 

4150. 

4150. 

9.0 

4147. 

4147. 

8.7 

4145. 

4145. 

8.3 

4142. 

4142. 

8.0 

4139. 

4139. 

3923. 
3927. 
3932  . 
393b. 
394  0. 
3944 . 

3  949. 
3953. 
3957. 
39b  1 . 
39b5  , 
3968  . 
3972. 
397b. 
3960. 
3963  . 
3967. 
3991  . 
3994. 
3996. 

4  001. 
4004  . 
4008  . 
4011. 
4014. 
4017. 
4020. 
4023. 
402b  . 
4029. 
‘.032  . 

4035  . 

4036  . 
4041  . 
4043  . 
404b  . 
4046. 
4051. 
4053  . 
405b  . 
4056. 
4061  . 

4  Ob  3  . 
4065  . 
40b7  , 
4069  . 
4071  . 
4073. 
4075. 
4077. 
4079  . 
4061  . 
4062. 
4084  . 
406b  . 
4067. 
4069  . 
4090. 
4091  . 
4093. 


lb. 71 

7.03* 

7.7 

lb.  71 

7.312 

7.3 

16.71 

7.  319 

7.0 

16.71 

7.026 

6.6 

16.71 

7.032 

6.3 

16.71 

7.338 

6.0 

lb. 71 

7. 04  4 

5.6 

16.71 

7.363 

5.3 

It.  71 

7.066 

5.0 

16.71 

7.369 

9.6 

lb.  71 

7.069 

9.3 

16.71 

7.063 

3.9 

16.71 

7.  372 

3.6 

16.71 

7.376 

3.3 

It.  71 

7.073 

2.9 

16.71 

7.  Oo  1 

2.6 

16.71 

7.069 

2.3 

16.71 

7 . 3  a  6 

1.9 

16.71 

7.337 

1.6 

16.71 

7.  339 

1.3 

16.71 

7.  390 

0.9 

lb.  71 

7.391 

0.6 

16.71 

7.091 

0.3 

lb. 71 
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10039. 

7607 

1. O4  5  1  75  89  215  . OQOu 

0.0905 

7091 . 11 3b 

.0010 

0. 1  ? 

0.  001 

1.2 

2500  . 

10*. 

.  0020 

0. 29 

0.003 

3.7 

2500. 

1 4  9  . 

.0030 

0. 49 

0.  009 

8.2 

2502  . 

2  4  7 . 

.  0090 

0.96 

0.020 

15.5 

2507. 

4  4  7. 

.005" 

1 .79 

0.090 

25.3 

2519. 

795. 

.00b" 

2.39 

0.071 

35.9 

2  5  4 1 . 

1323. 

.  00  70 

4 . 4  3 

0.112 

4  5 .  8 

25  71  . 

2032. 

.  0  "  5  " 

6.29 

0.  162 

53.9 

2690. 

2617. 

.00*0 

8 . 4  2 

0.219 

60.2 

3161. 

3252. 

.0100 

10.76 

0. 2o2 

65.0 

3  5  0  2  . 

3652. 

.  1  1  " 

13.2  6 

0.399 

53.6 

3  9  08 . 

402  8 . 

.  0120 

15.69 

0. 419 

71.3 

4322  . 

43  06 . 

.012? 

16.71 

G  .  4  4  1 

72.0 

4  4  32  . 

44  7  3  . 

.0130 

16.71 

C.  4  9  1 

72.3 

4&  52 . 

45  3  3  • 

.019" 

16.71 

0 . 5  s  3 

72.5 

4  9 1 7 . 

4  *  4  5  . 

.015" 

lb.  71 

C.636 

72.7 

5  1 4  1 . 

5155  . 

.  Olt  " 

lb. 71 

0.709 

72.3 

5331  . 

5350. 

.0170 

16.71 

0.731 

72.8 

59*1  . 

5507  . 

.0160 

16.71 

0.  3  5  4 

72.8 

5o26 . 

5b  4  0  • 

.01*" 

16.71 

0.  927 

72.6 

5  7  39  , 

5750. 

.020" 

16.71 

0.99* 

72.5 

5833. 

58  4  3 . 

.0210 

16.71 

1.072 

72.3 

5910. 

591  o  . 

.0 72  0 

16.71 

1  .  I44 

72.0 

5973  . 

5*7*. 

.  023  0 

16.71 

1.216 

71.6 

b  0  2  2  . 

602b  . 

.  O?4" 

16.71 

1.237 

71.5 

5060  . 

6055. 

.  025" 

16.71 

1.359 

71.1 

6066  . 

6092  . 

.0260 

16.71 

1 . 4  30 

70.  8 

5  10  6. 

6110  . 

.  02  70 

16.71 

1.500 

70.9 

5117. 

612  0. 

.  0  2  c  C 

lb. 71 

1.571 

70.0 

6121  . 

612  j. 

.  029  " 

lb.  71 

1.690 

6  9,6 

511b  . 

b  1 2  0  . 

.031" 
.  0320 
.033  0 
.034  " 
.  0350 
.03  oO 

•  037" 
.036" 

•  039" 
.  0  4  C  o 
.  DM  0 

•  0  4  2  ~ 
.  0  4  3  " 
.099  " 
.  O4  5  " 

.04b" 
.  0A70 
.  o^e  o 

.  0  4  9  O 
.  050" 
.0510 
.  052  0 
.  053  " 
.059" 
.  0  5  5  C 
.056" 
.0570 
.0580 


Is.  71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.7  1 

16.71 

16.71 

16.71 

16.  71 

16.71 

16.71 

16.71 


1.710 
1.773 
1.8-i7 
1.916 
1.963 
2.051 
2.117 
2.  16  4 
2.250 
2.315 
2.360 
2  .  4  4  5 
2.509 
2.5  72 
2.635 
2.698 
2.760 
2.522 
2.563 
2 . 9  4  4 
3 . 0  0  4 
3 . 0  6  4 
3.12  3 
3.152 
3 . 2  4  1 
3.299 
3.356 
3 . 4  1  3 
3.970 


6  9.2 

65.3 
66.  4 
67.9 

67.5 
67.0 

66.6 
66.1 
65.  7 

65.2 
b  4 . 7 

64 . 3 

63.5 

63.3 
6  2.9 

62. 4 
b  1 . 9 

61.5 
61.0 

60.5 
60.  1 
5  9.6 

59.2 
55.7 

55.3 

57.6 
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|  S//OT  €>S> 

hutzuc  enArffL  &£T4  =  -~1 

■ 

~f  S°  ?<-f  I/A  -norJ  At^O-U^ 

3  £££ 17 ?1-  t\OCS  f^occ £  DU—  OF  &A7r&ry 

O' <-  *<-OuJS  /fiy-ro  p*s-iAC>e_  C  \Cf/y=s 
Keco'CS 

Ao  -  l-ZSo  itS1-  £/tuc><-c_  or?i/c/CF 

V<2  =■  C?OOQ  /(\j~* 

Po  -  1500  /^/ 


g- 32. 1 /a 
P  I  -3 • 1 A1 592 
VC-11A7. 

UM-16.675 

OB-b.15 

AB -P I *DB*D  B/ A . 

WS-2A.*PI/(20.*DB) 

D-VOAB*UM*12. 

RL-8.75 
PIP-. 1081 
WP-96. 

VM-2710. 

PMAX-A750Q. 

RHOC-.Ob 

RTB-11563AOO. 

GAMMA-1. 233 
WC -2  6 . 

WA  -A. 

CK-1./7. 

GAMAM-GAMMA-1. 

GAMAP-GAMMA+1. 

GAMA 1-GAMAP/GAMAM 
GAMA?«2.*GAMMA/(1. -GAMMA) 

RATI0»l.*WC/(5.+WP ) 

RTORTB-GAMAM*  (  l./6„*(l.+CK)*wP/(2.**C))*VM*VM 

SI  G-D*SQRT  (  {  (  GAMAP/2  •  )  **GAMA1 )  /  (  GAMMA*RT0*RAT  10  )  )  /  (  AB*GAMAM*f>.  ) 
PM=.^C  +  RT0*RATI0M2./(3*D) 

WRITE  (♦»  '  (  5F13  .  A  )  • )  P M , R A T 1 0 , RTO , S I G , D 
WEFF*WP*(kC*WA)/2» 

AmVM/(l,-(PM*AB*G*UM)/(VM+VM*tiEFF)) 

B* ( ( A/VM)-1. ) ♦  UM 
TM=2 •  +UM / V  * 

U0«tXP(L0G(JM)-2.-JM/6) 

DR-2  . 

DP-2.75 
OS-1.75 
DP  ASS-1. 5 
WR-A000. 

BULK-150003 . 

BETA-. 7 
RH03IL-.0303 
THETA-75. ♦PI/180. 

UB-. 25 
DT». 001 

APR-PIMDP*DP-DK*3R)/A. 

ASP»PI*(DS*OS-DPASS*DPASS)/A. 

PE  0-2500 . 

PE-PEO 

PR-100. 

WO  I l«R  HDOI L*APR*Rl*12 .♦ ( A. ♦( PE*3.^PK l/BJLK ) 

US-.05 

W1-.  5 

W2-.25 

S-AOO. 

FGF«WR*(SlN(TMETA)-JB*COS(THETA)) 

PRLF-3000. 

E2-0. 

A02-0. 

A01-1.25 

VQ-6000. 


LK  2 • 1 • 3 
Pl-PE 
P2*PR 
XR*0. 

XRD*0. 

U«UD 
T*0 . 

P  1  K  *  0  . 

P2K-0. 

XRK«0. 

XR0L*0. 

TM1-TM+DT 
TM2*TM-0T 
50  Cl— 1, 

C2  *  1 . 

C  3  *  - 1  • 

DE  LP  1*0 . 

DE  LP  2*  0 . 

DELXR*0. 

DELXRD*0. 

IF(T.GT.TM)  G 3  TO  60 
IEIT.LE.TM2)  GO  TO  70 
IF(T.LT.TH)  DT»Tn-T 
IF  (  T.LT.  T*l  )  GO  TO  70 
DT*Tni-Tn~DT 
GO  TO  70 
60  DT»Tni-TK, 

70  DO  300  1*1, A 
C1=C1+2.*C2 
C2  =C  2-  .  5 
D1»1.5-C1*.5*C3 
C  3  *0 . 

P 1 E F  F  *P 1 4  Q l *9 ik 

P2EFF=P24C1*P2K 
XREFF=Xk+Dl*XRK 
XRDEFF*XRD*Dl*XkDL 
IF (P2EFF .GE .PE )  E2*.75 
El-1 .-E2 
OP 1*P1EFF-P= 

IFOPl.5E.PRLF)  A32  *0 . 

DP2*P2EFF-Pt 
IF (DP2.LE. G. )  DP2*0. 

AO  *A0 1  + A02 
B1*115.5*A3 

D31- <61*E1/  ( 1 2 . * A P R  )  )  **2  . 

A3* ( (RL-XREFF )/( DT*3ULK)  )**2  . 

B31»-(2.*XRDEFF*SCRT(A3)^D31/2.) 

C31*XRDEFF**2.-D31*DP1 

P1K*(— B31-S0RT ( B31+B31-4.*A3*C31 ) )/( 2  «  *  A  3 ) 

D32*(Bl*E2/(3b.*APR) )**2. 

B32»-(2.*XR3EcF*S0RT(A3)*D32/2.) 

C32»XRDEFF**2.-D32*DP2 

FLO»IND«B32*B32-<».*A3»C3  2 

I F ( F_  0WN3 . LE . 0. )  FL3WN0*0. 

P2K*(-B32-$3RT(FL0irfx3  )  ) / ( 2 • *  A3  ) 

XRK«DT*XRDEFF 

IF(T.GE.TM)  GO  TO  100 

0ELT-Dl*DT/2. 

VP -A*U/ ( B ♦ J ) 

BB-B*U-< VP»A)*DELT 


CL»-( (VK^A)*U>VK*d)*D;LT 
DELU-(-BB*S0RT<BB*B8-<i.*CC)>/2. 

0ELT2«<B*l3G<1.*DELU/U)*DELJ>/A 

UEFF*U*DELU»{Dl*DT-3EtT2)*A*<U*DELU)/(B«-U»DELU> 

FGAS*(WEFF/G)*A*A*B*UEFF/<B«'UEFF)**i. 

GO  TO  200 

100  TEFF»T-TMOl*0T 

FGAS-P«*AB*  l  l.-BETA>*  l  l.*TE  =  F/Sl  G  )  **GAi1A2 
200  FSEALl*PI*US*3R*(Wl*S*<i.*W2*{PlEFF*3»*P2EFF)  ) 
FSEAL2«PI*US*3S*{Wl*S*2.*w2*<PE*PR>) 

FSEAL*FSEA_1+FSEAL2 

XR3L*3T*(5/rtR)*(FGAS-APR*(PlEFF*3.*P2EFF)*ASP«,(PE*PR)-FSEAL'»FGF) 

0ELP1*DELP1*C1*P1K 

DELP2«DELP2*C1*P2K 

DELXR«DELXR*C1*XR< 

delxrd»delxrd*ci*xr3l 

3 C 0  CONTINUE 

P 1 *P 1  +  DE  LP  1/6. 

P2»P2*DEi.P2/6. 

XR*XR*DELXR/6. 

XRD«XRD-*DELXR3/b. 

U«UEFF 
T*  T ♦ 0 T 

DkOIL«RHOOlL*APR* ( Rl-XR>*12.*( A. ♦ ( PI +  3 . *P 2 ) / B JL K  ) -W3 1 L 
RH0E«RH001L*»l.+PE/bULK  ) 

DV«DW3Il/RH3E*p 1 +DS  *DS*XR*3. 
pE*PE0/( 1. ♦OV/VQ  )**CK2 
PRINT  AOO, T, U, XR» XR3,  PI ,P2 *  PE 

AOO  F0RMAT(3X,F5.4,5X,F5.2»5X»F5.3,5X»F5.1»5X*Fb.0»5X,Fb.0*5X,F6.0) 

IF ( XRD.LE . 0. )  GO  TO  500 
IF IXR.GE.8.75)  GO  TO  500 
GO  TO  50 
500  STOP 
END 
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82  9 86  LIGHTWEIGHT  TOWED  HOWITZER  DEMO 

PARTIAL  PHASE  2  VOLUH  (U>  FNC 
NORTHERN  ORDNANCE  DIV  R  RATHE 
SSIFIED  FHC-E-2041-VOL-C-PT-2  DAAAZi-88- 

4STRATOR  PHASE  i 
:ORP  MINNEAPOLIS 
:T  AL  APR  87 
-C-0047  F/G 

wjd 

1INN 

9/6 

2/1 

NL 

1 

- 

• 

G«32 .17a 
PI *3 .1^15^2 
V011A7. 

1^*16.675 

08*6.15 

AB*P I*Db*D3/A.  ' 

WS*2A.*PI / ( 20. *DB) 

D*VC*AB*U*1*12  • 

RL-S.75 
PIP*. 1081 
WP-96. 

Vr-2710. 

PKAX*A7500 . 

RHOC*. 05 
RTB*1 15b  3A  00 . 

GAHv,a  =  i.232 
WC-26. 

W  A  *  A  . 

CK*l./7, 

GA  WA**GA*:*.A-1. 

GAklAP*GAf'MA*l. 

GAHA 1*GAKAP/GAMAK 

GA  * A2  =  2  .*GAHNA/(  l.-GANNA  ) 

RATin*l.*wC/(5.*WP ) 

R  T  0  *  r  T  fc  — G  A  ‘'l  A  *1  +  (  1  •  /  b  .  ■♦■  ( 1  •  ♦  C  K  )  *  W  P  /  (  2  •  *  W  C  )  )  *  V  M  *  V  K 
SIG=n+SQPT(((GAMAP/2.H‘*GA'1Al)/(GAi*tA*RT0*RATI0))/(AB*GAMAf1*b 
P.m»WC*RTP*RaTI0*12./  (3*0) 

WPITF  (*,  •  (5F13.A)  * )  PM,RATI0,R7O»SI3»D 
a1  EFF=*P4-(*C+tfA)/2. 

A*  V H / (  l.-(Pvl*A9*G*J*1)/(VM*V'1*wEf-F  )  ) 

P=( ( A/V»i)-1. )+UM 
TM=2  ,*UK/ VN 

U0=S*P  (L0G(  J*)-2.-Ur./B) 

0R  =  2. 

DP  *  2 . 7  5 
DS  *  1 . 7  5 
DPASS*1.5 
WR  * A  000 . 

8ULK*1 50000 . 

BETA*. 7 
R  HOD  I L  * .  03  08. 

THETA=75.*PI /10O. 

UB*. 25 
DT  * . 00 1 

APR*P I *(DP*0P -DR *DR) /A. 

AS*P  I*DS*DS/A. 

AS*0. 

PE  0*  25  00 . 

PE*PE0 

PR*100. 

t40lL*RH00IL*A?R*RL*12.*(A.-f(Pt^3.*PR)/BJLK) 

US*. 05 
W1-.5 
W2-.25 
S«A00. 

FGF«WR*(S IN (THETA) -J3*CJS( THETA ) ) 

PRLF-3500. 

E2-0. 

A02-0. 

A01-.875 


V0-5000. 

CK2-1.3 

Pl-PE 

P2«PP 

XR«0. 

XRD*0. 

U«J3 

T«0. 

P  1 K  *  0 . 

P2<*0. 

XR<*0. 

XRDL*G . 

Tr 1*TM*DT 
TM2*TM-DT 
50  Cl  — 1. 

C2  *  1 . 

C3  =- 1 . 

delpi-o. 

DELP2*0. 

DELXR«0. 

DELXRD*0 • 

IF(T.GT.TM)  GO  TO  60 
I F  (  T  .LE  •  T*1 2  )  GO  TO  70 
IFtT.LT.TM)  dt-tm-t 
IF ( T  .LT.Tm  )  GO  TO  70 
OT  “TM-TM-DT 
GO  TO  70 
60  OT *TM1-TM 
70  00  300  1*1,4 
C1«C1*2.*C? 

C2*C2-  .5 

D1*1.5-C1*.5+C3 

C3*0. 

PlEFF»Pl-*-Dl*PlK 

P2EFF*P2+D1*P2K 

XREFF»Xk+01*XRK 

XRDEFF*XR0-Ol*XRDL 

IFtP2EFF.GE.PE)  E2-.75 

E1-1.-E2 

0P1-P1EFF-PE 

IF  (OPl.GE.PkLF)  A32*0 • 

DP2*P2EFF-PE 
IF (0P2.LE. 0. )  DP2«0. 

A0*A01 ♦ AC2 
B1«115.5*A0 

D31«(Bl*El/f 12.*APR) )**Z. 

A3«( (RL-XREFF )/ ( DT*6JLK ) )**Z. 

B31»-(2.*XR0E?F*S3RT(A3)+031/2.) 

C31«XRDEFF**2.-031*0P1 

P1K*(-B31-S3RT(B31*sJ31-A.*A3*C31))/(2.*A3) 
032»(Bl*E2/t36.*APR) )**2. 

B32— t  2.*XR3EFF*SQRTt  A3  l*3  32/2.  > 
C32*XR0EFF**2 .-332*0P2 
FL  OWN'D  «B  32*  B  3  2  -  A .  *A3*C32 
IFtFLOWNO.LE.O.)  FLOWND-0. 

P2K«(-B32-S3RT( FLOWND ) )/( 2.*A3 ) 
XRK«0T*XRDEFF 
IF  (T.GE.T1 )  GO  TO  100 
DELT«Dl*DT/2. 

VP«A*U/ <  B ♦ J ) 


9B*B-Mj-(VP*A)*D£tT 
CO-((VP*4)*U*VP*3»*DELT 
DELJ«(-b84$3RH  bB*BB-A.*CC> )/2. 

DELT2*(e*L0G< l.+DELU/ JJ+DELUi/A 

UEFF«U+DELU*(D1*DT-3ELT2)*A* ( U*DE LU ) / ( B DE LU ) 
FGAS*U'EFF/G)*A*A*B*UEFF/<  3*U£FF  >**3. 

GO  TO  200 

100  TEFF=T-TV*D1*0T 

FG  AS  xP*1*Afc*  (  l.-BETA)*(  1 . ♦ TE F F / S 1 G > *♦ GAM A2 
200  FSFAL1»PI*JS*0R*(W1*S*A.4«2*<P1EFF*3.*P2EFF>> 
FSEAL2*=PI*US*DS*(Wl*S*2.4i»2*(PE4PR)) 

FSEAL-FSEAL1 

XRDL=0T-MG/rtR)*(FGAS-APR*{PlEFF43.*P2EFF)4AS*(PE4PR)-FSEA^4FGF) 

r)ELPl*DELPl4Cl*PlK 

DELP2*DELP24C1*P2< 

DELXR=DELXR4C1*XR< 

DELXRD=DELXR04C1*XR3L 
300  CONTINUE 

Pl-PlOELH/6. 

P2*P2+DELP2/6. 

XR*XP+DELXR/b» 

XRDxXRD+DELXRO/c. 

U*UEFF 

T=T+DT 

Dft'3IL*RHDGlL*ADR*{  Rt-XP. )  ♦  12  .  *  (  A.  ♦  (  PI  +  3  .  *P  2  )  /BULK )  ~ri  D I  L 
RH0E*RHC0IL*C1.4pE/d'JLK) 

DV=DWIL/RH0E4AS  +  XR*12. 
pE*PEO/(  l.OV/V0)**CK2 
PRINT  AGO* T,U» *k * XR 3 ♦ P 1 ♦ P2 * P E 

A  00  F0RKAT(3X,  =  5.A*5X,F5.2oX,F5.3«5X,F5.1»5X,F&.0,5X,Fb.0,5X»Fb.0> 
IFIXRD.LE.O.)  GO  TO  500 
IF ( XR.GE.5. 75)  GO  TO  500 
GO  TO  50 
500  STOP 
END 
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16.71 

„  .1160 

16.71 

aw  .ii7o 

16.71 

«  .1180 

16.71 

3.240 

43.7 

6947, 

3.263 

43.1 

6827. 

3.326 

42.5 

6709. 

3.368 

42.0 

6594. 

3  ,  A 10 

41.4 

6482. 

3.451 

40.8 

6374. 

3.491 

40.3 

6269. 

3.531 

39.7 

6167. 

3.571 

39.2 

6068. 

3.610 

38.7 

5972. 

3.649 

38.2 

5880. 

3.  .,85 

37.6 

5791. 

3.723 

37.1 

5705. 

3.760 

36.7 

5622. 

3.7  97 

36.2 

5542. 

3. 933 

35.7 

54b5. 

3.869 

35.2 

5390. 

3.903 

34.7 

5319. 

3.938 

34.3 

5250. 

3.972 

33.8 

5163. 

4.005 

33.4 

5118. 

4.039 

32.9 

5056. 

4.071 

32.5 

499b. 

4.  104 

32.1 

4938. 

4.135 

31.7 

4882. 

4.167 

31.2 

4828. 

4.198 

30.8 

4776. 

4.228 

30.4 

4725. 

4.259 

30.0 

4676. 

4.289 

29.6 

4628. 

4.318 

29.2 

4582. 

4.347 

28.8 

4538. 

4.375 

28.4 

4495. 

4.404 

28.0 

4453. 

4.432 

27.7 

4412. 

4.459 

27.3 

4372. 

4.486 

26.9 

4334. 

4.513 

26.5 

4297. 

4.5  39 

26.2 

4260. 

4.565 

25.8 

4225. 

4.591 

25.4 

4191. 

4.616 

25.1 

4157. 

4.641 

24.7 

4125. 

4.665 

24.4 

4093. 

4.690 

24.0 

4062. 

4.713 

23.7 

4032. 

4.737 

23.3 

4003. 

4.780 

23.0 

3975. 

4.783 

22.6 

3947. 

4.805 

22.3 

3920. 

4.828 

22.0 

3894. 

4.849 

21.6 

3868. 

4.871 

21.3 

3843. 

4.892 

21.0 

3819. 

4.913 

20.7 

3795. 

4.933 

20.3 

3772. 

4.953 

20.0 

3749. 

4.973 

19.7 

3727. 

4.993 

19.4 

370b. 

5.012 

19.1 

3685. 

6946. 

2785. 

6827. 

2791. 

6709. 

2796. 

6594. 

2801. 

6482. 

2806. 

6374. 

2811. 

6269. 

2815. 

6167. 

2820. 

6068. 

2825. 

5972. 

2830. 

5880. 

2834. 

5791. 

2839. 

5705. 

2843. 

5622. 

2843. 

5542. 

2852. 

5465. 

2856. 

5390. 

2860. 

5319. 

2865  . 

5250. 

2859. 

5183. 

2873. 

5118. 

2877. 

505b. 

2881. 

499b. 

2885. 

4938. 

2889. 

4832. 

2892. 

4823. 

2896. 

477b. 

2900. 

4725. 

2904. 

4675. 

2907. 

4623. 

2911. 

4532. 

2914. 

4538. 

2918. 

4495. 

2921. 

4453. 

2925. 

4412. 

2928. 

4372. 

2931. 

4334. 

2935. 

4297. 

2938. 

4260. 

2941. 

4225. 

2944. 

4191 . 

2947. 

4157. 

2950. 

4125. 

2954. 

4093. 

2956. 

4052. 

2959. 

4032. 

2962. 

4003. 

2965. 

3975. 

2968. 

3947. 

2971. 

3920. 

2974. 

3694. 

2976. 

3866. 

2979. 

3843. 

2982. 

3619. 

2984. 

3795. 

2967. 

3772. 

2969. 

3749. 

2992. 

3727. 

2994. 

370b. 

2997. 

3685. 

2999. 
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XK 


mmv 


.1190 
.1200 
.1210 
.1220 
.  1230 
•  124  0 
.1250 
.1200 
.1270 
.1280 
.  129C 
.1300 
.1310 
.1320 
.1330 
.1340 

.1350 

.1360 

.1370 

.1380 

.1390 
.  1400 

.1410 
.1420 
.1430 
.1440 
.1450 
.1460 
.  1470 
.1480 
.1490 
.1500 
.1510 
.1520 
.1530 
.1540 
.1550 
.1560 
.1570 
.1580 
.1590 
.  1600 
.1610 
.1620 
.1630 
.  1640 
.1650 
.1660 
.1670 
.1680 
.1690 
.1700 
.1710 
.1720 
.1730 
.1740 
.1750 
.1760 
•  1770 
.1780 


16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 

16.71 


5.031 

5.050 

5.063 

5.086 

5.103 

5.  121 

5.138 

5.155 

5.171 

5.137 

5.203 

5.219 

5.234 

5.249 

5.254 

5.273 

5.292 

5.306 

5.320 

5.333 

5.346 

5.353 

5.371 

5.383 

5.395 

5.406 

5.413 

5.429 

5.439 

5.450 

5.460 

5.470 

5.480 

5.489 

5.493 

5.507 

5.515 

5.524 

5.532 

5.540 

5.547 

5.554 

5.561 

5.568 

5.574 

5.580 

5.585 

5.592 

5.597 

5.603 

5.607 

5.612 

5.616 

5.620 

5.624 

5.628 

5.631 

5.634 

5.637 

5.639 


18.8 

3665. 

18.4 

3645. 

18.1 

3625. 

17.8 

3606. 

17.5 

3586. 

17.2 

3570. 

16.9 

3552. 

16.6 

3535. 

16.3 

3519. 

16.0 

3502. 

15.7 

3487. 

15.4 

3471. 

15.1 

3456. 

14.8 

3441. 

14.5 

3427. 

14.3 

3413. 

14.0 

3400. 

13.7 

3366. 

13.4 

3374. 

13.1 

3361. 

12.8 

3349. 

12.5 

3337. 

12.3 

3325. 

12.0 

3314. 

11.7 

3303. 

11.4 

3293. 

11.1 

3282. 

10.9 

3272. 

10.6 

3263. 

10.3 

3253. 

10.0 

3244. 

9.8 

3235. 

9.5 

3227. 

9.2 

3219. 

9.0 

3211. 

8.7 

3203. 

8.4 

3195. 

8.2 

3188. 

7.9 

3181. 

7.6 

3174. 

7.4 

3168. 

7.  1 

3162. 

6.8 

3156. 

6.6 

3150. 

6.3 

3144. 

6.0 

3139. 

5.8 

3134. 

5.5 

3129. 

5.2 

3125. 

5.0 

3120. 

4.7 

3116. 

4.5 

3112. 

4.2 

3109. 

3.9 

3105. 

3.7 

3102. 

3.4 

3099  . 

3.2 

3096. 

2.9 

3093. 

2.6 

3091. 

2.4 

3089. 

3665. 

3001 

3645. 

3  00*. 

3625. 

3006 

3606. 

3008 

3538. 

3010, 

3570. 

3013, 

3552. 

3015. 

3535. 

3017. 

3519. 

3019. 

3502. 

3021. 

3437. 

3023. 

3471. 

3025. 

3456. 

3027. 

3441. 

3029. 

3427. 

3030. 

3413. 

3032. 

3400. 

303*. 

3386. 

3036. 

337*  . 

3037. 

3361. 

3039. 

3349. 

30*1. 

3337. 

30*2. 

3325. 

30**. 

331*. 

30*5. 

3303. 

30*7. 

3293. 

30*8. 

3282. 

3050. 

32.72. 

3051. 

3263. 

3052. 

3253. 

305*. 

324*. 

3055. 

3235. 

3056. 

3227. 

3058. 

3219. 

3059. 

3211. 

3060. 

3203. 

3061. 

3195. 

3062. 

3186. 

3063. 

3131. 

306*. 

317*. 

3065. 

3163. 

3066. 

3162. 

3067. 

3135. 

3068. 

3150. 

3069. 

31**. 

3070. 

3139. 

3070. 

313*. 

3071. 

3129. 

3072. 

3125. 

3073. 

3120. 

3073. 

3116. 

307*. 

3112. 

307*. 

3109. 

3075. 

3105. 

3075. 

3102. 

3076. 

3099. 

3076. 

3096. 

3077. 

3093. 

3077, 

3091. 

3078. 

3089. 

3073. 

•  179  0 
.1800 
.1810 
.1820 
.1830 

•  1  8  A  0 
.1850 
.1860 
.1870 
.1880 

FORTRAN 


lb. 71 
16.71 
16.71 
16.71 
16.71 
16.71 
16.71 
16.71 
16.71 
16.71 


STOP 


5.6*2 

5.6*4 

5.6*5 

5.6*7 

5.6*8 

5.6*9 

5.650 

5.650 

5.650 

5.650 


2.1 

1.9 

1.6 

1.3 

1.1 

0.8 

0.6 

0.3 

0.1 

>0.2 


3067. 

3085. 

306*. 

3083. 

3081. 

3081. 

3080. 

3060. 

3079. 

3079. 


3087. 

3085. 

308*. 

3083. 

3061. 

3061. 

3060. 

3080. 

3079. 

3079. 


3076. 

3078. 

3079. 
3079. 
3079. 
3079. 
3079. 
3079. 
3079. 
3079. 
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C/220 -f 

PART  NUMBER:  12585720-003,  Reservoir  Accumulator 

DESCRIPTION: 

The  reservoir  accumulator  size  and  pressure  requirements  were 
determined  -from  the  following  oil  displacement  analysis. 

Nitrogen  volume:  8000  cu.  in. 

Pre-Charge  Pressure:  120  psi  at  70  F 


STATUS: 

Mounting  and  dimensional  requirements  have  been  provided  to  York 
for  finalized  design. 


AUTHOR:  Jeff  Ireland 
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DESCRIPTION:  TOWING  STABILITY 

STATUS: 


net  change 


The  following  section  is  exerpted  from  the  Phase  I  Dynamic 
Analysis  Report.  All  towing  stability  analysis  was  performed  in 
Phase  I  since  it  was  felt  that  the  results  of  the  analysis  would 
si gni f i cant  1 y  with  the  configuration  changes. 


AUTHOR:  Scott  Dacko 
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Section  9 


Towing  Stability 

Extensive  work  has  been  performed  in  the  area  of  towing  stability.  The 
objective  is  to  maximize  LTHD  stability  and  ensure  that  its  perfomance 
in  towing  la  equal  to  or  greater  than  that  of  the  M198.  Both  static  and 
dynamic  responses  are  examined,  with  considerable  effort  placed  on 
dynamic  analysis.  Advanced  simulation  and  animation  software  was  used 
for  the  dynamic  analysis,  and  the  results  currently  exist  in  the  form  of 
graphs  in  Appendix  F  and  animations  on  videotape.  The  analysis  is 
presented  in  the  following  sections. 

9.1  Static 

Static  stability  for  both  the  LTHD  and  M198  were  analyzed.  The  center  of 
gravity  heights  for  both  vehicles  were  estimated  and  compared  and  the 
distances  between  wheels  were  also  compared.  The  results  are  shown  in 
Table  9-1.  The  data  used  to  determine  these  values  can  be  found  Appendix 
F.  The  important  result  for  comparison  purposes  is  the  tipping  angle. 
This  is  the  angle  at  which  each  vehicle  would  begin  to  tip  over  sideways, 
since  the  center  of  gravity  would  be  over  the  wheel  centers.  As  we  can 
see,  the  M198  tipping  angle  is  better  by  3.7  degrees.  This  is  attributed 
to  its  wider  wheel  base.  The  LTHD  purposely  had  a  narrower  wheel  base  to 
allow  it  to  track  in  the  tow  vehicle's  track  (M813).  The  LTHD  track  was 
increased  to  about  87  inches  (from  73  inches)  to  eliminate  this  problem 
on  June  3. 

9.2  Dynamic 

For  comparison  purposes,  more  Important  system  characteristics  exist  than 
those  described  in  the  static  analysis.  These  include  the  responses  of 
the  vehicles  at  various  speeds  to  various  road  conditions.  Before  these 
are  presented,  however,  lunette  loads  are  briefly  covered.  Then  the 
responses  to  bumps  on  flat  road  and  bumps  on  roads  while  rounding  a 
corner  are  described. 

9.2.1  Lunette  Loads 

The  maximum  loads  to  be  experienced  by  the  lunette  were  determined  for 
designing  a  lunette  with  the  best  size  and  shape  and  to  determine  the 
impact  on  barrel  deflection.  Load  conditions  were  estimated  by 
accounting  for  the  LTHD  moment  of  inertia  for  yaw  and  pitch,  maximum 
deceleration  by  the  M813,  jack-knifing  situations,  and  the  effect  of 
large  bumps  on  a  typical  truck  suspension.  It  was  determined  that  the 
barrel  can  easily  withstand  these  forces,  which  are  considered  worst-case 
conditions. 


Lpgfl  Condition 
Forward-back  load 
Vertical  load 
Lateral  load 
Vertical  fc  lateral 


M«  Lflfld 
9,000  lbs 
2,200  lbs 
6,750  lbs 
7,100  lbs 


FftCL.gr  Pf  Safety  Design  to; 

2  18,000  lbs 

4  8,800  lbs 

2  13,500  lbs 

2.27  16,115  lbs 
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TOWING  STABILITY  -  STATIC 

Mi  98 

TOW  LTHD 

C.G.  HEIGHT.  IN.  53.84  48.29 

DIST.  BE7W.  WHEELS 

CENTERS.  IN  93.0  73.0 

TIPPING  ANGLE.  DEG'S  40.8  37.1 


9.2.2  Dynamic  Analysis  Model 

This  section  documents  dynamic  towing  stability  analyses  performed 
comparing  the  FMC  LTHD  and  the  M198.  The  work  was  performed  by  Don 
Cronquist,  FMC  Central  Engineering  Laboratories.  Important  dynamic 
parameters  are  found  in  Table  9-2. 

The  purpose  of  the  dynamic  analysis  and  simulation  is  to  prowide 
preliminary  towing  stability  comparisons  of  the  LTHD  and  M198  using 
simplified  models.  The  results  also  demonstrate  PMC's  capabilities  for 
performing  such  analyses. 

9. 2. 2.1  Computer  Models 

Towing  stability  simulations  were  run  with  computer  models  of  the  LTHD 
and  M198.  The  simulations  were  run  using  version  8A  of  the  ADAMS 
computer  program,  a  program  for  simulation  of  three-dimensional  non¬ 
linear  mechanical  systems  of  rigid  bodies.  This  program  is  currently  in 
use  by  a  number  of  companies  in  the  automotive  industry  and  elsewhere. 

Figures  9-1  and  9-2  illustrate  the  computer  models  of  both  towed 
howitzers  obtained  using  the  graphics  provided  by  the  ADAMS  program. 

Each  towed  howitzer  model  consists  of  the  following: 

.  A  model  data  file  which  providea  the  mass  and  inertial  properties, 
geometric  data  including  force  location  and  direction,  part-to-part 
connection  data,  and  graphics.  Some  of  the  terrain/obstacle  data  is 
also  provided  in  the  model  file. 

.  A  set  of  user-written  subroutines,  written  in  FOETEAH,  which  Include 
the  tire/ terrain  model,  and  a  subroutine  used  to  guide  the  trailers  in 
the  cornering  simulations. 

Standard  ADAMS  elements  were  used  in  most  portions  of  the  vehicle  model. 
The  tire/terrain  model  is  not  a  standard  element,  however.  This  model  is 
a  point  contact  tire  capable  of  running  in  any  direction  on  a 
cross-country  terrain  surface.  A  tire  obstacle  of  rectangular  cross 
section  may  be  superimposed  on  the  of  the  terrain  surface.  The  tire 
model  computes  approximate  force  of  interaction  between  the  tire  and 
obstacle  in  both  the  horizontal  and  vertical  directions.  The  tire 
obstacle  is  a  new  feature  which  was  added  to  the  tire  model  for  this 
dynamic  analysis  project. 

The  effect  of  the  obstacle-enveloping  capability  of  the  tire  is  included 
in  the  tire/obstacle  interaction  model,  but  is  not  included  in  the 
tire/terrain  interaction  model.  In  this  way,  vehicle  maneuvering  may  be 
done  with  efficient  use  of  computer  resources,  but  obstacle-enveloping 
characteristics  may  be  included  when  the  tire  is  interacting  with  an 
obstacle  of  rectangular  cross-section. 

Previous  simulations  using  the  tire  model  used  a  wash-board  shaped 
terrain  in  which  the  terrain  height  varied  only  along  one  dimenaion  of 
the  model.  Future  simulations  may  be  performed  in  which  terrain  height 
is  a  function  of  both  horizontal  coordinates. 
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TOWING  STABILITY  -  DYNAMIC 


M198 

TOW 

LTHD 

TIRE  SPRING  CONSTANT. 
LB/IN 

*4,507 

1.312 

DAMPING  COEFFICIENT 
ROLLING*  NEGLIGIBLE 

.019 

.019 

TOW  WEIGHT.  LBS 

15.600 

8.982 

MOMENT  OF  INERTIA 

ABOUT  ROLL  AXIS. 
FT-LB-SEC*2 

2096.5 

562.7 

H 


TOW  SPEED*  25  MPH 


Uli 


In  the  present  computer  simulations,  the  base  terrain  was  flat,  but  an 
obstacle  of  rectangular  cross  section  vas  superimposed  on  the  terrain 
surface  to  test  the  vehicle  stability  limits.  7or  this  phase  of  the 
project,  approximate  tire/obstacle  force  aquations  vere  used.  However, 
tire  test  data  may  be  used  to  improve  the  accuracy  of  the  tire/obstacle 
forces  for  future  work. 

Representative  model  data  files  for  the  LTHD  and  M198  are  ahovn  in 
Appendix  L.  The  input  data  is  included  in  the  appendix  with  each  tire 
model. 

The  tire/terrain  model  is  not  shown  in  this  report  due  to  its  proprietary 
nature.  This  program  was  already  in  use  at  IWC  before  its  use  in  the 
LTHD  program.  As  part  of  the  LTHD  and  M198  study,  the  tire  model  vas 
enhanced  to  compute  horizontal  and  vertical  force  of  interaction  of  the 
tire  with  an  obstacle  of  rectangular  cross-section.  The  portion  of  the 
tire  model  is  shown  in  Appendix  L. 

It  should  be  noted  that  an  error  vas  detected  in  the  equations  for  the 
horizontal  force  between  the  tire  and  the  obstacle.  This  error  allowed  a 
rearward  force  between  the  tire  and  the  trailing  corner  of  the  bump.  The 
error  is  documented  in  Appendix  L. 

The  effects  of  the  error  in  horizontal  obstacle  force  vere  generally  very 
small  in  the  cases  where  the  simulations  vere  re-run.  The  largest  effect 
occurred  in  the  LTHD  straight  road  simulations  in  which  the  roll 

angle  increased  from  roughly  23.5  degrees  to  25.5  degrees  in  one  case  and 
In  another  case  the  maximum  roll  angle  decreased  from  roughly  13  degrees 
to  11.2  degrees.  These  comparisons  are  documented  in  Appendix  H. 

9. 2. 2. 2  Assumptions 

The  following  assumptions  vere  made  in  this  analysis: 

.  All  parts  are  rigid  bodies. 

.  All  joints  sre  ideal  kinematic  joints  without  looseness  or 
flexibility. 

.  The  tires  have  constant  spring  rates  up  to  the  rim- radius,  at  which 
the  spring  rate  increases  by  a  factor  of  5.  The  tire  forces  vere 
computed  based  on  the  position  of  the  center  of  the  tire  spindles 
relative  to  the  terrain  surface  and  relative  to  the  obstacle. 


.  These  simulations  do  not  Include  the  effects  of  rolling  resistance  or 
lateral  spring  rate. 

.  The  tov  hitch  of  the  towing  vehicle  had  the  following  Idealized 

notion:  constant  speed  in  either  a  straight  path  or  a  straight  path 
followed  by  a  constant  radius.  There  is  a  30  inch  transition  between 
the  straight  path  to  the  radius.  The  tow  hitch  was  represented  by  a 
spherical  joint  without  friction  or  looseness. 

.  Zero  tire  damping  was  assumed  in  the  straight  path  bump  simulations. 
Damping  in  the  tire/terrain  forces  was  1.9%  of  critical  in  the 
cornering  simulations  (based  on  the  appropriate  tire  mass  and  spring 
rate).  In  all  cases,  no  tire  damping  was  applied  during  the 
tire/obstacle  Interaction. 

.  The  terrain  surface  was  rigid  and  flat.  A  rigid  rectangular  bump  was 
superimposed  on  the  terrain  surfaces. 

9. 2. 2. 3  Results 

The  simulations  show  that  the  LTHD  tolerates  larger  terrain  obstacles 
that  the  M198  when  being  towed  along  a  straight  road.  However,  in 
cornering  simulations  the  M198  tolerated  larger  tire  obstacles  that  the 
LTHD  at  18  and  20  mph,  while  the  LTHD  tolerated  larger  obstacles  at  15 
mph. 

A  modified  version  of  the  LTHD  with  the  width  between  tire  centers 
increased  from  73  Inches  to  87.25  inches  was  able  to  tolerate  a  higher 
bump  than  the  M198  in  a  cornering  simulation  at  20  mph. 

Figures  9-3  and  9-4  summarize  the  simulation  results.  Figure  9-3  shows  a 
plot  of  maximum  roll  angle  versus  bump  height  for  the  straight  road 
towing  simulations.  This  plot  compares  the  performance  of  the  LTHD  and 
KL98  and  reveals  that  for  equally  sized  bumps  the  the  LTHD  does  not  roll 
as  far  as  the  M198.  The  simulations  also  show  that  the  M198  rolls  over 
for  a  12  inch  obstacle  while  the  LTHD  returned  to  an  upright  position. 

Figure  9-4  susnarizes  the  results  for  cornering  maneuvers.  This  is  a 
plot  of  towing  speed  versus  maximum  obstacle  height  without  rollover. 
Where  an  obstacle  height  is  Indicated  on  this  plot  it  means  that  a 
simulation  was  run  at  this  bump  height  and  at  the  next  higher  bump  height 
and  that  the  towed  howitzer  rolled  over  only  at  the  higher  bump  height. 
The  bump  height  Increment  for  these  simulations  was  2  Inches,  except  that 
simulations  were  also  run  with  a  bump  height  of  1  inch. 

Figure  9-4  shows  that  at  18  and  20  miles  per  hour  the  Ml 9 8  trailer 
tolerated  higher  bumps  than  the  LTHD.  However,  at  15  miles  per  hour  the 
LTHD  tolerated  a  higher  bump  than  the  M198.  These  simulations  were 
performed  with  a  500  inch  radius  at  the  tow  hitch.  Figures  L-l  and  L-2 
in  Appendix  L,  are  plots  of  roll  angle  vs.  time  for  the  LTHD  and  M198 
straight  road  towing  simulations.  Each  figure  shows  the  response  for 
three  different  obstacle  heights.  Figures  L-3  through  L-5  are  plots  of 
roll  angle  vs.  time  for  the  cornering  simulations.  Figure  L-3  shows  the 


BUMP  HEIGHT  (INCH) 


•table  response  of  tbe  LTHD  to  an  8  inch  obstacle  and  an  unstable 
response  to  a  10  inch  obstacle  at  15  Biles  per  hour.  Figure  L-4  shows 
the  stable  response  of  the  M198  to  a  6  inch  obstacle  and  an  unstable 
response  to  an  8  inch  obstacle  at  15  aph.  Figure  L-5  shows  the  stable 
response  of  the  M198  to  a  1  inch  obstacle  and  the  unstable  response  to  a 
2  inch  obstacle  at  20  aph. 

Two  aodifications  of  the  initial  LTHD  configuration  were  also  sianlated: 

.  Lowered  CG  6.0  inches  from  48.29  inches  to  42.29  inches.  This  is  a 
12. 4%  reduction  in  CG  height.  Note  that  these  CG  heights  are  measured 
relative  to  the  bottom  of  the  tires  at  their  unloaded  radii. 

.  Width  between  tire  centerlines  increased  from  73  inches  to  87.25 
Inches.  This  is  a  19. 5X  increase  in  width.  This  modification  was 
simulated  both  with  the  original  tire  spring  rate  of  1500  lb/in  per 
tire  and  with  a  higher  spring  rate  of  3000  lb/ in  per  tire. 

The  LTHD  with  the  first  modification  (lowered  CG)  was  found  to  roll  over 
with  1  and  2  inch  obstacle  heights  in  the  20  mph  cornering  maneuvers. 

This  is  the  same  as  with  the  initial  LTHD  configuration.  The  Increment 
between  theae  obstacle  heighta  (1  inch)  was  thus  too  large  to  detect  an 
improvement  due  to  this  modification. 

The  LTHD  with  the  second  aodification  (greater  width  between  tire 
centerlines)  was  found  to  be  stable  with  a  2  inch  obstacle  height  at  20 
mph.  Note  that  the  original  LTHD  rolled  over  with  a  1  inch  obstacle  and 
the  M198  rolled  over  with  a  2  inch  obstacle.  Thus  the  cornering 
stability  of  the  LTHD  with  Increased  width  between  tires  was  superior  to 
both  the  original  LTHD  and  the  M198. 

Appendix  L  documents  the  error  in  the  horizontal  component  of  force  while 
the  tire  is  interacting  with  the  tail  end  of  the  bump.  In  moat  of  the 
cases  checked,  this  error  did  not  significantly  affect  the  maximum  roll 
angles.  The  largest  effect  was  found  in  the  case  of  the  LTHD  straight 
road  towing  siaulations  in  which  one  aaxlmum  roll  angle  was  decreased 
from  10.7  degrees  to  9.2  degrees  and  in  another  case  the  maximum  roll 
angle  was  increased  from  23.5  degrees  to  25.5  degrees.  (See  Figures  L-8, 
L-9  and  L-10  in  Appendix  L.)  However,  in  the  M198  straight  road 
siaulations  there  were  much  smaller  changes  in  the  aaximum  roll  angle 
simulations  in  the  cases  where  both  tire  aodels  were  simulated  and  the 
results  compared.  (See  Figures  L-6,  L-7,  L-ll  and  L-12  in  Appendix  L. 

The  aaxiaua  roll  angle  versus  bump  height  plot  of  Figure  9-3  is  based  on 
the  corrected  bump  force  subroutine  except  for  the  H198  simulation  at  the 
6  inch  bump  height.  At  the  data  in  the  cornering  maneuver  plot  of  Figure 
9-4  was  based  on  the  original  subroutine.  However,  it  is  not  likely  that 
the  results  of  Figure  9-4  will  be  affected  by  the  error  because  of  the 
•mall  effect  this  error  has  on  aaximum  roll  angle  in  the  cornering 
simulations  checked  and  because  of  the  large  increment  in  obstacle 
heights  in  Figure  9-4. 

It  should  be  noted  that  both  the  LTHD  and  M198  became  highly  sensitive  to 
tire  obstacle  as  the  cornering  speed  approached  their  respective 
stability  limits. 
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9. 2. 2. 4  Considerations  for  Additional  Analysia 
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The  following  ideas  surfaced  upon  performing  the  towing  stability 
analysis.  The  ideas  are  not  necessarily  recommended;  rather,  they  could 
be  considered  as  options  for  further  study. 

.  Strong  consideration  could  be  given  to  reducing  the  ratio  of  the  CG 
height  to  width  between  tire  centerlines  for  the  LTHD,  such  as  in  the 
case  of  the  LTHD  with  87.25  inches  width  between  tire  centers.  Rote 
that  a  6  inch  (12.4%)  reduction  in  CG  height  did  not  produce  a 
significant  improvement  in  stability,  while  a  14.24  inch  (19. 5%) 
increase  in  width  between  tire  centerlines  did  produce  a  significant 
Improvement.  This  recommendation  la  intended  to  increase  the  speeds 
and  obstacle  heights  required  to  overturn  the  vehicle. 

.  Consideration  could  be  given  to  uaing  the  tow  hitch  to  restrain  the 
LTHD  from  rolling  to  large  angles.  The  writer  believes  that  a  small 
restraining  torque  at  the  tow  hitch  could  prevent  the  LTHD  from 
overturning  in  cases  vhen  small  obstacles  are  encountered  at  speeds 
near  the  stability  limits.  If  this  restraint  is  provided  by  means  of 
a  slip  clutch,  a  suitable  maximum  torque  could  be  selected  which  will 
prevent  damage  to  either  vehicle  and  which  would  be  unlikely  to 
overturn  the  tow  vehicle. 

The  second  option  is  intended  to  make  the  LTHD  less  prone  to  rollover 
without  warning  as  the  turning  stability  limit  is  approached.  Hot  only 
will  larger  obstacles  be  required  to  overturn  the  towed  howitzer,  but  the 
driver  will  be  somewhat  able  to  feel  the  torsional  response  of  the  towed 
howitzer  in  instances  in  which  it  almost  turned  over.  The  driver  will 
thereby  be  able  to  learn  about  the  howitzer's  towing  stability  limits 
without  a  howitzer  rollover  actually  occurring. 

In  this  section,  towing  stability  concerns  have  been  examined.  A 
detailed  analysis  has  been  performed  on  the  LTHD's  towing  stability  and 
by  estimating  the  H198's  CG  location  and  inertial  properties,  a  similar 
analysis  was  performed.  More  complex  models,  which  build  upon  the  models 
presented  here,  can  be  developed  for  further  work  in  the  area  of  towing 
stability. 
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Appendix  L 


Tovlng  Stability  Data,  Models  and  Results 
All  figures  are  Boll  Angle  vs.  Time: 


Figure  L-l  LTHD,  25  aph, 
Figure  L-2  M198,  25  aph, 

Figure  L-3  LTHD,  15  aph, 
Figure  L-4  M198,  15  aph, 

Figure  L-5  M198,  20  aph, 

Figure  L-6  M198,  25  aph, 

Figure  L-7  M198,  25  aph, 

Figure  L-8  LTHD,  25  aph. 
Figure  L-9  LTHD,  25  aph, 
Figure  L-10  LTHD,  25  aph, 
Figure  L-ll  M198,  15  aph, 
Figure  L-12  M198,  15  aph, 

LTHD  Model  . 

LTHD  Data  . 

LTHD  Prograa,  ADAMS  Model 

M198  Model  . 

Ml 9 8  Data  . . 

M198  Program,  ADAMS  Model 
FORTRAN  Subroutines  for  Til 


straight  direction  .... 
straight  direction  .... 

cornering  . 

cornering . L- 

comering  . L- 

straight,  corrected  .  L- 

straight,  corrected  .  L- 

straight,  corrected,  10  inch  .  L- 

straight,  corrected,  12  inch  .  L-l 

straight,  corrected,  14  inch  .  L- 
comering,  corrected  .  .  . 
straight,  corrected  .... 


Model  . 
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CG  AND  MMI 

OF  LTHD  TOWED  - 

PART  A  ONLY 

TOTAL  WBIGHT  .  8379 

C.G.  X  COORD  (IN)  - 

.859291E-01 

C.G.  T  COORD  (IN)  . 

3.11609 

C.G.  Z  COORD  (IN)  . 

91.7685 

JYZ  (FT-LB- 

■S*2)  . 

24701. 

,7 

JXY  (FT-LB- 

S*2)  . 

362.078 

JXZ  (FT- LB- 

■S'2)  . 

24729. 

.2 

DESCRIPTION  WEIGHT 

X 

CANNON 

2600 

0 

BREECH 

500 

.5 

BAND 

86 

0 

AUTOPR 

45 

0 

MBRKASSEM 

137 

0 

ASSTGUN 

59 

0 

GUNNER 

115 

0 

SLIDEUNIT 

396 

0 

SLUNITLP 

182 

0 

PLATFORM 

265 

0 

SPADES 

190 

0 

SPLKBLTS 

8 

0 

TSHFTUPR 

10 

0 

TSHFTLWR 

30 

0 

TRLLKPR 

16 

0 

WRSKPLTFM 

20 

0 

LFTGEYES 

14 

0 

GIMBALASSY 

256 

0 

PULLYASSY 

16 

0 

LTRYLASSY 

18 

0 

TRLASSYLH 

215 

-36 

TRLASSYRH 

215 

36 

MNFASSYG 

95 

0 

TRVCYL 

50 

31 

HDLNSL 

40 

-18 

HOSETSG 

20 

-18 

MNPOLD 

190 

0 

LFTEYSMF 

16 

0 

TRLLKPF 

14 

0 

ELVCYLLH 

25 

-14 

ELVCYLRH 

25 

14 

COMPCBL 

25 

0 

IBCUSH 

26 
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ACCBNKASSY 

300 

0 

HFLUID 

150 

0 

RCYLREC 

440 
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RCYLNREC 

160 

0 

CRCYLREC 

220 

0 

CRCYLNREC 

160 

0 

FVDYK 

100 

0 

REARYOK 

100 

0 

LTBUNGY 

10 

0 

DMF 

68 

0 

DCTVB 

38 

0 

DLVB 

48 

-27 

DRVB 

48 

27 

DTPIN 

8 

0 

DTHANDL 

3 

0 

DPVDLPIN 

4 

0 

DBRKCAL 

44 

0 

AOVEROILAC 

15 

0 

-26 

52 

DBKROTOR 

42 
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-26 

52 

DHUBS 

36 
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-26 

52 

DBRGSBLTS 

16 

0 

-26 

52 

BQUILCYLLH 

225 

-21 

12 

192 
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12 

192 
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LTHO  -  440  INCH/SEC  ON  10"  BUMP 
FILE  LTHD-OAT 

INITIAL  VELOCITY  =  440  INCHES/SEC  3° 

NOTE:  EULER  ANGLE  DISPLACEMENT  OUTPUT  IS  YAW#  PITCH#  ROLL 
PA/1#GR 

A/9#  EU=  90.#  90.#  0. 

ARKER/11#EU=90#90#0  Z-AXIS  PAR.W/X-AXIS 

MARKER/12#EU=0#-90#0  ..7-AXIS  PAR.V/Y-AXIS 

MARKER/ 13 #EU=0#0#0  Z-AXIS  PAR.W/Z-AXIS 

GRAPHICS  OF  TERRAIN: 

MA/5001#  QP=0. ,  36.5#  0. 

MA/5011#  QP  =  20  00. #  36.5#  0. 

GR/14#OUSS5001#  5011 
MA/5002#  flP=0. ,  -36.5#  0. 

MA/5004#  S)P=440.#  -36.5#  0. 

MA/5006#  QP=44 0. #  -36.5#  10. 

MA/5008#  «P=470.#  -36.5#  10. 

MA/5010#  fiP=470. #  -36.5#  0. 

MA/5012#  QP=20 00. #  -36.5#  0. 

GR/15#  OU=50  01#  5002#  5004#50  06#5008#5010# 5012# 5011 
HOWITZER  BODY 

PA/2#MA=8484.#CM=200#IP=  53.48E6#  53.31E6#  2.33E6 
#QG=0.#0. #0. #EU=90#90#0#VX=440 
MA/200#QP=0.#48.29#-216.5 
MA/223#  QP=0. #  48.29#  -216.5#  EU=0.#  -90.#  -90.  INITIALLY  PARALLEL  W/  GLOBAL 

ON  GROUND#  INITIALLY  AT  VEHICLE  C.G. 


MA/10#  PA=1 #  OP=  -216.5#  0.#  48.29 
MA/201#QP=12.5#76.5#-51.6 
MA/203#QP=12.5#76.5#-265.5 
MA/205#OP=12.5#99.#-265.5 
MA/207#flP=12 .5#99. #-326.35 
4/2 09# QP=1 2. 5 #33. #-326. 3 5 
A/211#QP=12.5#33.#-51.6 
A/202#«P*-12.5#76.5 
A/204#QP=-12. 5 #76. 5 
A/206#QP=-12. 5 #99.# 

A/208#  QP=-12. 5#99.# 

A/210#  QP=-12. 5#33»# 

A/212#QP=-12- 5#33«# 

MA/215#8P=12#33.#-51 
MA/217#QP=12#15-#-51 
MA/219#QP=12#15.#-26 
MA/221#QP=12#33.#-26 
MA/216#QP=  -12#  33.# 

MA/218#  QP=  -12#  15.# 

MA/220# QP=  -12#  15.#  -267.6 
MA/222#«P*  -12#  33.#  -267.6 
*"*  GR/1# OU=2 01# 203# 205 # 2 07# 209# 211# 2 01 
#202#204#206#208#210#212 
#202/203# 204/205# 206/20 7# 208/ 209 #21 0/211 
#212/215# 217#218#216#222#220# 219# 221 
#215/217# 218/219#220/221#22 2 
MA/213#«P=0.#35.75#-51.6 

GR/2#CY#CM=213#LENGTH=28.85#RADIUS=3.75#SIDES=8 


8 

¥ 


$ 


1 


& 


#76.5# 

-51. 

6 

#76.5# 

-265 

.5 

#99.#- 

265. 

5 

#99.#- 

326. 

35 

#33.#— 

326. 

35 

#33.#— 

51.6 

.#-51. 

6 

-#-51. 

6 

.#-267 

.6 

.#-267 

.6 

33.# 

-51. 

6 

15.# 

-51. 

6 

& 


LEFT  SIDE  WALKING  BEAM: 

»  PA/3#MA*244.#CM=300#IP»  2.78E5#2.78E5#4.64E4 
m  >G=  -243.6#  21.5#24.75#  EU=90#90#0#VX=440 
WALKING  BEAM  GRAPHICS: 

MA/300#fiP*ll. 34#  -5.38#  0. 

5S  MA/301#  «P*  9.#  -9.#  36. 
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£ 

& 

t 


w 

$ 


k> 

■*  4W 

H< 

3 

c 


i 


& 

5j 

v3 


s 


& 


HA/303# 

8P= 

9.  #  —  3. #  36 • 

HA/305# 

QP= 

9.#  7.#  0. 

HA/307# 

8P= 

9.#  -3.#  -36. 

HA/309# 

0P= 

9.#  -9.#  -36 

MA/311# 

8P= 

-9.#  -9.#  36. 

A/313# 

8P= 

-9.#  -3.#  36. 

A/315# 

8P= 

-9.#  7.#  0. 

HA/317# 

«P= 

-9.#  -3.#  -36. 

HA/319# 

QP= 

-9.#  -9.#  -36 

GR/3#OU 

=301# 

303#305#307#309#301# 

#303,313/305 

#3 15/307# 31 7/309# 319 

RIGHT 

SIDE  WALKING  BEAM: 

PA/4#HA 

=2  44. 

#CM=400#IP=2.78E5#2. 

#QG=  -243.6# 

-21.5,24.75 

#  EU=90# 

90#0# 

VX=44  0 

WALKING  BEAH  GRAPHICS: 

HA/400# 

8P= 

-11.34#  -5.38#  0. 

HA/401# 

8P= 

-9.#  -9.#  36. 

HA/403# 

8P= 

-9.#  -3.#  36. 

HA/405# 

OP= 

—9.#  7.#  0. 

HA/407# 

8P= 

-9.#  -3.#  -36. 

HA/409# 

QP= 

-9.#  -9.#  -36 

HA/411# 

8P= 

9.#  -9.#  36. 

HA/413# 

8P= 

9.#  -3.#  36. 

HA/415# 

8P= 

9.#  7.#  0. 

HA/417# 

8P= 

9.#  -3.#  -36. 

HA/419# 

8P= 

9.#  -9.#  -36 

GR/4#0U 

% 

H 

O 

«# 

tl 

403#405#407#409#401# 

#  403 #413/405 #415/40 7# 41 7/40  9#  419 
TOWING  VEHICLE: 

PA/5#HA=10.#CM=500#IP=1.E3#1.E3#1.E3#8G=0#0#0 

#EU=90#90#0#VX=440 

HA/500#  8P=0#0#0 

JO/l#RE#  1=1001#  J=1002 

HA/1001#  PA=3#  8P=0#0#0#  EU=9 0.#  90.#  0. 

HA/1002#  PA=2#  8P=21. 5#  24.75#  -243.6#  EU=  90.#  90.#  0. 
GR/9#  Cl#  CH=1002#  RADIUS=3. 

JO/2#RE#  1=2001#  J=2002 

HA/2001#  PA=4#  8P=0#0#0#  EU=90.#  90.#  0. 

HA/2002#  PA=2#  8P=-21.5#  24.75#  -243.6#  EU=  90.#  90.,  0. 
GR/10#  Cl#  CH=2002#  RADIUS=3. 

JO/3#  SP#  1=3001#  J=3002 
HA/3001#  PA=5#  8P=0#  32#  0 
HA/3002#  PA=2#  8P=0#32#0 
GR/11#  Cl#  CH=3002#  RADIUS=3. 

JO/4#  TR#  1=4002#  J=4  00 1 

HA/4001#  PA=1#  QP=10.#  0.#  32.#  EU=90.#  90.#  0. 

HA/4002#  PA=5#  QP=0.#  32.#  10. 

GR/12#  Cl#  CH=4001#  RADIUS=3. 

GR/13#  Cl#  CH=4002#  RADIUS=3. 

„£/l#  CO#  J0=4#  PAR=  440.#  0. 


LEFT  FRONT  TIRE 
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MA/1#  PA=3#  8P=15. #  -6.75#  36. 

HA/5#  PA*3#  QP=14. #  -6.75#  36 

FO/11#  AOF#  1=1#  J=ll#  PAR=  1500.#  10.#  1#  5#  0. 

FO/21#  AOF#  1=1#  J=12;  PAR=  1500.#  10.#  1#  5#  0, 

F 0/31#  AOF#  1=1#  J=13#  PAR=  1500.#  10.#  1#  5#  0. 

TIRE  GRAPHICS: 

A/320#  PA=3#  8P=  9.#  -6.75#  36.#  EU=90#90#0 
GR/5#  CY#  CH=320#  LENGTH=12#  RADIUS=18. #SIDES=8 


RIGHT  FRONT  TIRE: 

HA/2#  PA=4#  fiP=-15. #  -6.75#  36. 

MA/6#  PA=4#  <1P=-16. #  -6.75#  36. 

FO/12#  AOF#  1=2#  J=ll#  PAR=  1500.#  10.#  2#  6#  0. 

FO/22#  AOF#  1=2#  J=12#  PAR=  1500.#  10.#  2#  6#  0. 

FO/32#  AOF#  1=2#  J=13#  PAR=  1500.#  10.#  2#  6#  0. 

TIRE  GRAPHICS: 

HA/420#  PA=4#  8P=  -9.#  -6.75#  36.#  EU=90#-90#0 


GR/7#  CY#  CM=420#  LENGT H=12 #  RADIUS=18. #SIDES=8 


LEFT  REAR  TIRE: 

HA/3#  PA=3#  GP=15.#  -6.75#  -36. 

HA/7#  PA=3#  8P=14.#  -6.75#  -36. 

FO/13#  AOF#  1=3#  J=ll#  PAR=  1500.#  10.#  3#  7#  0. 

FO/23#  AOF#  1=3#  J=12#  PAR=  1500.#  10.#  3#  7#  0. 

F 0/33#  AOF#  1=3#  J=13#  PAR=  1500.#  10.#  3#  7#  0. 

TIRE  GRAPHICS: 

HA/321#  PA=3#  8P=  9.#  -6.75#  -36.#  EU=90#90#0 
GR/6#  CY#  CH=321#  LENGTH=12#  RADIUS=18. #SIDES=8 


RIGHT  REAR  TIRE: 

A/4#  PA=4#  QP=-15. #  -6.75#  -36. 

A/8#  PA=4#  8P=-16. #  -6.75#  -36. 

FO/14#  AOF#  1=4#  J=ll#  PAR=  1500.#  10.#  4#  8#  0. 

FO/24#  AOF#  1=4#  J=12#  PAR=  1500.#  10.#  4#  8#  0. 


;y  FO/34#  AOF#  1=4#  J=13#  PAR=  1500.#  10.#  4#  8#  0. 
‘V  TIRE  GRAPHICS: 


HA/421#  PA=4#  8P=  -9.#  -6.75#  -36.#  EU=90#-90#0 
GR/8#  CY#  CM=421#  LENGT H=12#  RADIUS=18. #SIDES=8 


TIRE  HOHENTS#  LEFT  FRONT: 
F0/41#A0T#I=1# J=ll 
FO/51#AOT#I=l#J=12 
FO/61#AOT  #1=1# J=13 


TIRE  HOHENTS#  RIGHT  FRONT: 
FO/42# AOT  #1=2# J=ll 
FO/52#AOT  #1=2# J=12 
F 0/6 2# AOT  #1=2# J=13 


TIRE  MOMENTS#  LEFT  REAR: 
FO/43# AOT#I=3# J=ll 
FO/53# AOT#I=3# J=12 
F 0/6 3# AOT  #1  =  3# J=13 


ft/.  TIRE  MOMENTS#RIGHT  REAR: 

^  F0/44#A0T#I=4#J=11 
V54#A0T  #1*4# J=12 
J/64#A0T  #1  =  4# J=13 


LEFT  FRONT  BUMP  STOP: 


tVXYVVY  ■ 


F0/1#C0#I*323# J*223#K=10.E3#L=6. 
MA/323#PA*3#9P=-5.#  0.#  36. 

MA/223#  PA=2#  9P=  16.5#  40.7#  -208.82 


I 


LEFT  REAR  BUMP  STOP: 

0/3#  CO#  1*325#  J*225#  K=10.E3#  L*6. 
\/325#  PA=3#  9P=  -5.#  0.#  -36. 
MA/225#  P A=2 #  9P=  16.5#  40.7#  -278.38 


RIGHT  FRONT  BUMP  STOP: 

F0/2#C0# 1=324# J=224#K=10.E3#L=6. 

^  MA/324#PA=4#9P*  5.#  0.#  36. 

MA/224#  PA=2#  9P=  -16.5#  40.7#  -208.82 


(J  RIGHT  REAR  BUMP  STOP: 

to  FO/4#  CO#  1=326#  J=226#  K=10.E3#  L=6. 

MA/326#  PA=4#  OP=  5.#  0.#  -36. 
g  MA/226#  PA=2#  9P=  -16.5#  40.7#  -278.38 

SPLINES  CLATERAL  FORCE  CHARACTERISTIC) 
•  SP/1#X=0.#  1.#  2.#  3.#  4.#  5.#  20.#  90. 
Js  #Y=  .0#  .16#  .30#  .40#  .49#  .57#  .75#  .85 


SYSTEM  AND  OUTPUT  CARDS: 
SYSTEM/  GC=386.088#  KGRAV=  -386.088 
OUTPUT/ SC  I ENT  I FIC#GRSAVE#SAVERE9# YPR 


33 


SLIP  ANGLES  IN  DEGREES 
RATIO:  FSIDE/FLOAD 


RE/1#  D#  1=  223#  J*10#:DISP  OF  CG  (YAW#PITCH#ROLL) 
RE/2#  V#  1=  200#  J=10#:VEL  OF  CG 
RE/3#  A#  1=  200#  J=10# : ACCEL  OF  CG 

j|(  c/41#F#I*l#:LEFT  FRONT  TIRE 
■  RE/42#F,I=2#:R1 GHT  FRONT  TIRE 

RE/43#F#I=3#:LEFT  REAR  TIRE 
ifl  RE/44#F#I=4#:RIGHT  REAR  TIRE 

RE/21#D#I=1#J=9 

»RE/22#D#I=2#J=9 
RE/23#D#I=3#J=9 
RE/24#D#I=4#J=9 


•T-C 


END  CARD: 


S 
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Ajl^a  AfoP^  fot-  rrA«»urr 

^  4^/^^TA-n. 


£ 


?r 


1 

v  HMI  OF  H198  IN  TOV  860507 

5?  TOTAL  WEIGHT  -  15780 

C.G.  X  COORD  (IN)  -  1 06084 E- 01 


C.G.  Y  COORD 

(IN)  - 

t 

C.G.  Z  COORD 

(IN)  - 

$ 

JYZ  (FT-LB-S'2) 

JXY  (FT-LB-S* 2)  - 

JXZ  (FT-LB-S* 2)  - 

* 

:v 

DESCRIPTION 

WEIGHT 

VHEELR 

263 

VHEELL 

263 

FIRING  BASE 

74 

SPADERT 

184 

SPADEL 

184 

TRAILRT 

674 

& 

TRAILL 

665 

BOTTCAR 

1555 

«*  , 

BBANDB 

5008 

*  o 

RECSYS 

2105 

CRADLE 

1020 

INTCOMP 

1276 

£< 

REMAINDER 

2509 

. 883412E-02 

21.9205 

14136.1 

2096.5 

12935.7 


X 

Y 

Z 

-46.5 

-34.34 

7 

46.5 

-34.34 

7 

0 

-8.94 

-157 

-36 

-17.14 

-96 

36 

-17.14 

-96 

-18.6 

-36.64 

-48 

18.6 

-36.64 

-48 

0 

-24.14 

6 

0 

21.86 

84 

0 

16.86 

46 

0 

17.86 

56 

0 

-3.64 

•  -2 

0 

-18.64 

-51 

* 


r*  - 

in 
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ptt-e;  MMt.bAT.  S3 

H198  -  440  INCH/SEC  ON  10"  BUMP 
FILE  M198.DAT 

INITIAL  VELOCITY  »  440  INCHES/SEC 
NOTE:  EULER  ANGLE  DISPLACEMENTS  IN  YAW#  PITCH#  ROLL 
*»A/1#GR 

A/9#  EU=  90.#  90.#  0. 

\RKER/11#EU*90#90#0  . . . .E-AXIS  PAR.V/X-AXIS 

MARKER /12#  EU=  0 #-90  # 0  . Z-AXI  S  PAR.W/Y-AXIS 

MARKER/ 13 #EU=0#0#0  . Z-AXIS  PAR.W/Z-AX1S 

GRAPHICS  OF  TERRAIN: 

MA/S001#  GP=0. #  46.5#  0. 

MA/5 Oil#  8P=2000.#  46.5#  0. 

GR/14#OU*5001#  5011 
MA/5 002#  fiP=  0. #  -46.5#  0. 

MA/5004#  «P*440.#  -46.5#  0. 

MA/5 006#  ®P=44  0«#  -46.5#  10. 

MA/5 008#  «P=470.#  -46.5#  10. 

MA/5 01 0#  GP*47 0.#  -46.5#  0. 

MA/5012#  8P=  20  00. #  -46.5#  0. 

GR/15#  OU*50 01# 5002# 5004 #5006# 500 8# 5010# 5012 #5011 
HOWITZER  BODY 

PA/2#MA=  15600. #CM=200#IP=  58.71E6#  53.14E6#  9.72E6 


MA/203, 


0.  #0. 

•  #EU=90 

#90#0 

#VX=44 

0 

«p= 

0.#  53 

.84# 

-242. 

8P= 

0.#  53 

■  84  # 

-242.# 

EU=0.# 

-90.#  -90. 

PA=1  i 

-  8  P=  - 

242.# 

0.#  5 

3 

.84 

FOR  VEHICLE 

«P* 

37.5# 

77.# 

-174. 

«P= 

37.5# 

77.# 

-286. 

QP= 

37.5# 

15.# 

-286. 

6  P= 

37.5# 

is.# 

-174. 

QP= 

37.5# 

25.# 

-174. 

«P= 

0.# 

40.# 

0. 

QP= 

0.# 

30.# 

0. 

8P= 

0.# 

75.# 

-286. 

9 

* 

© 

ii 

Z3 

LU 

i #  —167 • #  0. 

8P= 

-37.5# 

77.# 

-174. 

«P= 

-37.5# 

77.# 

-286. 

8P= 

-37.5# 

15.# 

-286. 

8P= 

-37.5# 

15.# 

-174. 

8P= 

-37.5# 

25.# 

-174. 

MA/210#  8P=  -37.5#  25.#  -174. 

GR/1#  OU=  207#201#203#205#207#213#211#209#210#211#213 
#208#202# 204# 206# 208/201# 202/ 203# 204/20 5 #206/20 7# 208 
GR/2#CY#  CM=215#  LENGTH=  138.#  RADIUS®  4.5#  SIDES*8 

S  TOWING  VEHICLE: 

PA/5#MA=10.#CM=500# IP=1«E3# 1. E3# 1 -E3 # 8G=0 #  0 #  0 
V '  #EU=90#90#0#VX=440 
mK  MA/5 00#  «P=0#0#0 

J? 

JO/3#  SP#  1=3001#  J=3002 
C-C"  MA/3001#  PA=5#  QP*0#  30#  0 
MA/3002#  PA=2#  «P=0#30#0 
GR/11#  Cl#  CM=3002#  RADIUS=3. 

JO/4#  TR#  1=4002#  J=4 00 1 
X/4001#  PA*1#  ®P=10. #  0.#  30.#  EU*90.#  90.#  0. 

^  .A/4002#  PA=5#  «P=0. #  30.#  10. 

^  GR/12#  Cl#  CM=4001#  RADIUS=3. 

GR/13#  Cl#  CM=4002#  RADIUSO.  , 


BARREL 


St 


GE/1#  CO#  J0=4 #  PAR*  440.#  0. 
LEFT  TIRE  FORCE: 


§ 


% 

dC 


3 


•>: 

ft 


|( 


% 


r 


<:* 


v. 

v 

'a 


4 


4/5. 


PA*  2# 

OP* 

46.5# 

19.5 

#  -250. 

PA*  2# 

OP* 

45.5# 

19.5 

#  -250. 

AOF# 

1*1# 

J*ll# 

PAR* 

5150.# 

10.# 

1# 

5# 

0. 

AOF# 

1*1# 

J*12# 

PAR* 

5150.# 

10.# 

1# 

5# 

0. 

AOF# 

1*1# 

J*13# 

PAR* 

5150.# 

10.# 

1# 

5# 

0. 

TIRE  GRAPHICS: 

HA/320#  PA=2#  OP*  37.5#  19.5#  -250.#  EU*90.#  90.#  0. 
GR/5#  CY#  CH*320#  LENGTH*  16.5#  RADIUS*  19.5#  SIDES=8 


HA/2# 

HA/6# 


RIGHT 

TIRE  FORCES: 

PA=2# 

QP=  -46.5# 

19.5 

#  -250. 

PA*  2# 

OP*  -47.5# 

19.5 

#  -250 

#  AOF# 

1*2#  J*ll# 

PAR* 

5150.# 

10.# 

2# 

6  9 

0. 

#  AOF, 

1*2#  J*12# 

PAR* 

5150.# 

10.# 

2# 

6  9 

0. 

#  AOF# 

1=2#  J=1 3# 

PAR* 

5150.# 

10.# 

2# 

6  9 

0. 

TIRE  GRAPHICS: 
HA/420#  PA=2#  OP*  -37.5# 
GR/7#  CY#  CM*4 20#  LENGTH* 


19.5#  -250.#  EU*90 . # 
16.5#  RADIUS*  19.5# 


-90.# 

SIDES* 


TIRE  HOHENTS# 
FO/41#  AOT  #I*1#J*11 
F0/51#A0T#I=1# J=12 
F0/61#A0T  #1*1# J*13 


LEFT: 


TIRE  HOHENTS#  RIGHT 
FO/42#AOT#I*2#J*ll 
0/5 2# AOT  #I*2#J*12 
0/62#A0T#I*2#J=13 


SPLINES  (LATERAL  FORCE  CHARACTERISTIC) 
SP/1#X*0.#  1.#  2.#  3.#  4.#  5.#  20.#  90. 
#Y*  .0#  .16#  .30#  .40#  .49#  .57#  .75#  .85 


SLIP  ANGLES  IN  DEGREES 
RATIO:  FSIDE/FLOAD 


SYSTEH  AND  OUTPUT  CARDS: 

SYSTEH/  GC=386 .088#  KGRAV*  -386. 088#HMAX=.005#ERR*. 00008 
OUTPUT /SC IENT I FIC#GRSAVE#SA VEREQ#YPR 


RE/1# 

D# 

I* 

223# 

J=10#: 

RE/2# 

V# 

1  = 

200# 

J*10#: 

RE/3# 

A  9 

I* 

200# 

J=10# : 

RE/41#F#I=1#:LEFT  FRONT  TIRE 
RE/42#F#I*2#:RIGHT  FRONT  TIRE 


Ki  RE/21#D#I=1#  J=9 
RE/22#D#I=2#J=9 


V  RE/11#V#I=1#J=9 


RE/12. 


1*2# J=9 


RE/31#A#I*1#J*9 

°E/32#A#I=2#J=9 


j§  END 


END  CARD: 
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4 


S 


h' 

/kN 


% 


a 

j 

oc 

or 


5>s 


»f  ~ 

U 

H 

u  « 

*- 

(O 

vs 

> 

« 

> 

• 

% 

in 

2 

CO 


OZ 

4  O 


OO 

o»- 


oa 


INI 

« 

« 


a 

Of 

• 

UJ 


CO 

> 

* 


a 


o 

v 

2 


0. 

UJ 

oc 

oc 

►H 

♦“ 

DO. 

a 

X 

4 

V 

UJ 

or 

H 

3 

or 

H-  X 

O 

HH 

A 

>- 

K 

w 

04 

V 

►- 

fV 

or 

o 

1 

OO 

X 

oc 

v 

o 

Z 

>  V 

M 

V 

3 

> 

UJ 

VU_ 

OC 

3 

>- 

z 

a 

ioh 

OC 

H- 

x 

o 

UJ  « 

*4H- 

V 

O 

o 

X 

KVO 

UJ 

O 

> 

h-  a 

♦ 

m> 

or 

> 

-S. 

% 

X 

>-  % 

H 

V 

A 

z 

^  4 

u. 

VO 

H- 

CL 

UJ 

*"*0 

HH 

OC 

X 

a. 

X 

►- 

ft 

x  V 

V 

•o 

V 

•M  1 

4/> 

o 

VI 

3 

UJ  CO 

. 

or 

> 

X 

••  H 

OX 

•  o 

•  > 

OC  N 

oca 

« 

N 

aocisi 

UJ  O 

W  J- 

or  ♦ 

'  j 

UJ  I 

»  N 

vac  • 

u^o 

OZ 

>o 

UJH*- 

X 

H- 

H*  vO 

OC  UJ» 

a  uj 

MJ 

IL  V 

ocm 

H 

OC 

or  uj  s 

Ooc  ft 
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PART  NUMBER:  12505714,  Traverse  Actuator 

DESCRIPTION:  TRAVERSE  ACTUATOR 

Piston  Diameter:  4,25  inches 

Rod  Diameter:  2.5  inches 

Bear  Loc  Force:  30,000  lbs 

This  actuator  was  sized  based  on  -firing  on  a  10’/.  side  slope,  with 
huncane  -force  winds,  a  muzzle  brake  imbalance,  and  the  rifling 
torque  input  into  the  gimbal,  traverse  actuator  and  platform 
assemb 1 y . 

The  reader  is  also  referred  to  5ection  C/260,  Tube  Laving 
Accuracy  for  additional  work  performed  in  finalizing  the  traverse 
cylinder  design. 

STATUS: 

Mounting  and  size  requirements  for  the  traverse  actuator  have 
been  determined  and  were  provided  to  York  to  be  finalized. 


AUTHOR:  Jeff  Ireland 
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PART  NUMBER:  12585767,  Tube  Interface 
DESCRIPTION:  TUBE  INTERFACE 

The  tube  is  a  modified  XM284  in  shape  and  was  chosen  to  provide 
the  lightest  weight  tube  possible.  The  drawing  specifies  FMC ' s 
requested  tube  weight  of  2,480  lbs. 

The  tube  OD  contains  five  raised  flat  areas  for  interfacing  with 
the  collar  sets  via  torsional  and  alignment  keys.  The  rear-most 
flat  has  a  machined-in  keyway  to  anchor  the  collars  and  rails. 
The  other  flats  do  not  have  keyways  to  allow  the  collar-rail  set 
to  "float"  a::  i  ally  when  temperature  changes  cause  thermal 
expansion.  The  locations  of  the  first  flat  (closest  to  the 
breech)  and  the  fourth  flat  allow  the  front  sides  of  the  first 
and  fourth  collar  sets  to  be  0.75  inches  in  front  of  the  front 
sides  of  the  manifolds  in  fire  position.  This  provides  a 
straight  path  for  the  transmission  of  rifling  torque  from  the 
collars  to  the  manifolds  at  the  time  of  maximum  torque.  The 
second  flat  is  positioned  to  perform  this  same  function  in  the 
loao  position  to  directly  transmit  rifling  torque  in  a  cookoff 
situation.  The  fifth  flat  position  allows  the  front-most  collar 
set  to  be  flush  with  the  front  manifold  when  the  tube  is  in  tow 
posi t i on . 

The  alignment  keys  (on  top  and  bottom  of  tube)  between  the  tube 
and  collars  maintains  side-to-side  alignment  of  the  tube.  The 
torsional  keys  transmit  the  rifling  torque  from  the  tube  to  the 
collars.  Space  between  the  alignment  and  torsional  keys  and  the 
collars  allows  the  pressure-induced  radial  tube  growth  at  the 
time  of  firing  prevents  stressing  of  the  collars 


STATUS: 

A  finalized  Tube  Interface  drawing  (TBF',  Dwg .  12585767)  has  been 

pjrovided  to  and  accepted  by  Benet. 


AUTHOR:  Scott  Dacko.  Joe  Turek.  Bart  Anderson 
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■■£2  2"  IrTI0N:  TUBE  LAYING  ACCURACY 

T  ~  e  ’-'eader  15  referred  to  the  enclosed  Gun  Laying  Accuracy  Report 
t  o'  details  o-f  the  analysis  per-formed. 


2  T AT,JS : 

T K s  tube  laving  accuracy  study  is  current  and  complete.  No 
additional  wort  in  this  area  is  anticipated  in  Phase  II. 


AUTHOR:  Sean  Marel 
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LTHD  Gun  Laying  Accuracy 

The  following  page*  tunarize  the  reaults  froa  the  analysis 
of  the  LTHD  *  s  gun  laying  accuracy.  Bxaaple  results  are 
contained  in  Appendix  A,  the  initial  model  parameter 
calculations  are  contained  in  Appendix  B  and  the  DADS*  model  is 
in  Appendix  C.  The  analysis  uses  a  DADS  mechanism  model  to 
simulate  the  gun  during  laying  exercises.  The  model  includes 
the  effects  of  the  hydraulic  cylinders,  control  valves, 
mechanism  geometry  and  friction  for  both  the  traversal  and 
elevation  systems.  Since  only  small  motions  were  simulated  the 
equilibration  system  is  modeled  as  a  constant  unbalanced  torque 
operating  on  the  gun. 


Friction 

No  major  problems  due  to  frictional  loses  in  the  system  were 
discovered  and  in  fact  some  beneficial  effects  were  found.  The 
friction  present  in  the  system  provides  a  significant  degree  of 
damping  which  will  enhance  the  system’s  controllability. 
Furthermore,  there  is  enough  friction  in  the  system  to  prevent 
typical  unbalanced  torques  from  moving  the  gun  when  the  locking 
system  is  released.  This  makes  the  timing  between  joystick 
motion  and  lock  release  less  critical. 


Valves 

The  control  valve  size  appears  to  be  the  most  critical 
factor  affecting  the  operator's  ability  to  accurately  lay  the 
gun.  For  precise  gun  laying  motion  rates  of  0.3  to  0.5  mil/sec 
are  necessary.  The  flow  rates  corresponding  to  these  velocities 
can  create  significant  lags  between  initial  motion  of  the 
joystick  and  the  motion  of  the  gun.  The  lag  is  due  to  the  tine 
it  takes  for  the  pressure  rise  in  the  cylinder  to  overcome 
unbalanced  loads  such  as  friction.  A  lag  of  1  to  3  seconds  is 
common  and  lags  up  to  8  seconds  are  possible  when  the  gun's 
direction  is  reversed  or  extended  inactivity  has  allowed  the 
cylinder  pressures  to  leak  down.  Much  smaller  tine  lags  occur 
if  notion  in  a  direction  has  just  occurred  and  resumption  of 
that  notion  in  the  same  direction  is  attempted. 


Basic  Problem 

The  time  lag  between  initial  joystick  motion  and  initial  gun 
notion  will  cause  the  operator  to  open  the  valve  more  than 
necessary.  Then  when  motion  does  occur  it  will  be  too  quick  and 
he  will  overshoot  the  desired  position  because  of  his  reaction 
tine  and  the  gun's  delayed  response  when  the  joystick  is 


returned  to  the  null  position.  The  size  of  the  tine  lag  has  a 
strong  dependence  on  the  friction  present,  the  voluae  under 
coapression,  the  effective  bulk  aodulus  and  the  flow  rate 
through  the  valve  consequently,  the  tiae  lag  cannot  be 
eliainated  and  a  way  to  deal  with  it  aust  be  found. 


Possible  Solution 

The  operator  displaces  the  joystick  enough  to  create  notion 
in  a  reasonable  tiae.  As  soon  as  the  notion  starts  he  lets  go 
of  the  joystick  so  that  the  gun's  notion  is  stopped.  The 
cylinder  is  now  pressurized  to  the  point  at  which  notion  will 
occur.  To  resuae  the  notion  the  operator  can  either  use  a 
saaller  joystick  notion  or  aore  likely  he  can  pulse  the  joystick 
until  the  desired  position  is  reached. 


Systea  requirements 

For  the  preceding  procedure  to  work  the  systea  aust  be 
configured  such  that; 

A. )  The  slowest  controlled  notion  rate  possible  with  the 
valve  is  about  0.3  to  0.5  ail/sec  when  a  saall  constant  joystick 
displacement  is  used  in  precise  gun  laying.  This  value  can  be 
doubled  if  the  joystick  is  pulsed. 

B. )  The  dynamic  flow  response  of  the  valve  should  be  nearly 
proportional  to  the  joystick  position  for  saall  stick 
displacements .  That  is  there  should  be  no  daaping  of  snail 
valve  notions  near  its’  closed  position. 


Secondary  Problems 

Flexible  hydraulic  lines  and  other  design  factors  which 
lower  the  systea’ s  effective  bulk  aodulus  coapound  the  delay 
problea.  The  tiae  it  takes  to  pressurize  a  cylinder  before 
notion  occurs  is  inversely  proportional  to  the  bulk  aodulus 
while  the  systea* s  rate  of  notion  is  independent  of  the  bulk 
aodulus.  Therefore  as  the  bulk  aodulus  is  decreased  it  takes 
longer  for  the  gun  to  start  aoving  but  once  it  does  start  its 
notion  rate  depends  aainly  on  the  joystick  position.  This  will 
add  to  the  operator’s  tendency  to  use  too  large  of  a  joystick 
displaceaent  when  initiating  notion. 

Friction  and  other  unbalanced  loads  such  as  equilibration 
also  increase  the  delay  between  joystick  and  gun  notion. 
Frictional  losses  provide  a  significant  daaping  force  to  the  gun 
and  their  elimination  could  create  other  difficulties.  Without 
this  daaping  there  is  a  auch  better  chance  that  pulsing  the 


joystick  near  the  gun’s  natural  frequency  could  lead  to  an 
instability  due  to  resonance.  Since  too  auch  friction  can  be  a 
problea  soae  optiaua  frictional  force  is  possible  but 
determination  of  that  optiaua  is  difficult. 


Valve  Sizes 

For  traversal  the  aaxiaua  desired  rate  is  80  ail/sec. 

Taking  into  account  the  losses  in  the  system  and  the  unbalanced 
loads,  a  valve  that  provides  5  gpm  of  flow  at  3000  psi  with  no 
pressure  drop  across  the  load  produces  the  desired  rate.  The 
equations  used  for  this  estimation  and  their  deviation  are 
included  as  the  next  three  pages.  To  generate  a  .3  ail/sec 
ainiaua  velocity  the  ainiaum  valve  opening  should  be 
0.25  percent  of  full  open.  Since  0.1  percent  is  possible  no 
problea  with  this  valve  size  is  apparent. 

For  elevation  the  valve  size  is  highly  dependent  on  the 
equilibration  system.  The  last  time  the  valve  was  sized  it 
needed  to  be  capable  of  delivering  40  gpm  at  3000  psi  with  no 
pressure  drop  across  the  load.  Generating  a  motion  rate  of 
0.3  ail/sec  to  0.5  ail/sec  requires  a  valve  opening  of  0.04  to 
0.06  percent  of  full  open.  This  is  a  2000  to  1  ratio  between 
aaxiaua  and  ainiaua  controlled  flow  which  could  be  outside  of 
the  valve’s  practical  limits  therefore  potential  problems  may 
exist  for  this  valve.  Either  some  degree  of  controllability 
aust  be  sacrificed  or  some  time  must  be  added  to  the  time  cycle 
and  the  valve  down  sized. 
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k'  Examples  of  Results  from  Analysis 
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In  figure  1  the  green  lines  represent  the  elevation  system's 
response  at  72  degrees  QB  to  two  constant  valve  openings  of  0.5 
and  0.05  percent  of  full  open.  The  first  vertical  red  line 
■arks  the  time  when  the  valve  is  first  opened.  The  second 
vertical  red  line  marks  the  time  when  the  0.5  percent  valve 
opening  was  closed.  The  target  position  is  0.5E-3  radians  The 
valve  is  not  closed  until  sometime  after  the  target  position  is 
reached  because  of  the  mechanical  and  human  delays  in  the 
system.  The  horizontal  red  line  and  the  third  vertical  red  line 
show  the  valve  position  when  a  0.05  constant  valve  opening  was 
used. 


Attention  should  be  paid  to  the  time  lag  between  when  the 
valve  is  opened  and  when  gun  motion  starts.  For  the  first  curve 
this  lag  is  about  0.1  seconds  and  for  the  second  curve  this  lag 
is  about  1  second.  Attention  should  also  be  paid  to  the  amount 
of  overshoot  that  is  present.  For  the  first  curve  it  is  about 
0.9E-3  radians,  (0.9  mils)  and  for  the  second  curve  the 
overshoot  is  only  about  0.15E-3  radians,  (0.15  mils).  The 
overshoot  could  be  reduced  further  with  a  smaller  valve  opening 
but  then  the  initial  time  lag  would  become  excessive. 


Figures  2  and  3  go  together,  figure  2  is  the  traversal 
response  of  the  gun  to  pulses  of  the  joystick  and  figure  3  shows 
those  pulses.  The  pulses  occur  every  0.5  seconds  with  a  0.5 
percent  aagnitude  and  a  0.1  second  duration.  In  reality  these 
pulses  could  be  saoother  in  nature  and  different  in  frequency 
and  duration  however  these  differences  should  not  have  a 
significant  effect  on  the  gun  behavior  shown  in  figure  2. 

As  shown  the  gun  takes  over  2  seconds  to  get  aoving  and 
another  5  to  6  seconds  to  traverse  0.5  ails.  In  aany  cases  this 
response  would  be  too  slow  however  this  can  be  corrected  by 
using  pulses  that  are  larger  in  aagnitude  and/or  longer  in 
duration.  The  iaportant  point  aade  by  figure  2  is  that  a  lower 
precision  valve  can  be  used  for  accurate  positioning  if  a 
pulsing  technique  is  used.  Recall  that  for  precise  gun 
traversal  the  recoaaended  lower  liait  for  a  steady  controlled 
valve  opening  is  0.25  percent  of  full  open  yet  it  appears  that 
at  least  twice  this  value  is  acceptable  when  the  valve  is 
pulsed . 


Time  5ec 


Figure  4  shows  the  gun's  traversal  response  to  a  1.0  percent 
valve  opening  when  the  effective  bulk  Modulus  is  175000  psi  and 
50,000  psi.  The  first  value  is  typical  for  a  system  that  uses 
steel  hydraulic  lines  and  has  mild  air  entrainment.  The  second 
value  is  typical  for  a  systea  that  contains  flexible  hydraulic 
lines  or  excessive  air  entrainment.  The  first  vertical  green 
line  marks  when  the  valve  was  first  opened,  the  vertical  red 
line  marks  when  the  valve  was  closed  for  the  175000  psi  bulk 
modulus  system  and  the  second  vertical  green  line  marks  when  the 
valve  was  closed  for  the  50000  psi  bulk  modulus  system.  In  both 
cases  the  target  position  was  0.5  mils. 

Two  observations  are  easily  made  from  this  figure.  First, 
as  the  bulk  modulus  decreases  the  initial  time  lag  increases 
proportionately.  Second,  the  average  elevation  rate  is 
unchanged  as  demonstrated  by  the  equal  average  slope  of  the  two 
traces.  Consequently,  as  the  bulk  modulus  decreases  the  minimum 
valve  opening  necessary  to  provide  adequate  reaction  time  is 
unchanged  but  the  gun's  response  becomes  sluggish.  This  can 
lead  to  excessive  joystick  motion  by  the  operator. 
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construction 

Type  SFL  self-lubricating  lined  Spherical  Plain  Radial  Bearings  are 
unit  assemblies,  each  consisting  of  an  inner  ring  with  a  spherical 
outer  surface  and  an  outer  ring  with  a  spherical  inner  surface.  Both 
rings  are  manufactured  from  through-hardened  steel.  The  outer  ring 
has  a  single  fracture  to  permit  assembly.  The  spherical  inner  sur¬ 
face  of  the  outer  ring  is  lined  with  a  mixture  of  solid  lubricant  TFE 
fibers  and  a  phenolic  resin  binder.  The  spherical  outer  surface  of 
the  inner  ring  is  chrome  plated.  Each  bearing  has  two  seals  made  of 
reinforced  nitrile  rubber. 

CHECK  FOR  AVAILABILITY 


dimensions 

The  controlling  dimensions  for  type  SFL  bearings  are  in  inches. 
Inch-metric  conversions  given  are  for  the  convenience  of  the  user. 

For  tolerances  of  nominal  dimensions  see  Tables  2  and  2A  on 
page  3. 

a  is  the  maximum  tilting  angle  for  effective  sealing.  If  this  angle 
Is  exceeded,  the  seal  lips  will  slide  off  the  spherical  surface. 
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Type  SFL  bearing 


D 

bearing 

desig¬ 

nation 

B, 

inner  nng 

B 

outer  nnj 

dt 

tohohcal 

load  ratings 
radial 

mast  a 

tilting 

o.d. 

width 

width 

diameter 

limit  dynamic 

angle 

(nominal) 

(nominal) 

(nominal) 

Crwf.) 

load  load 

(appr.)  (mas.) 

19.050 

.7500 

31.750 

1.2500 

750.12 

16.66 

.656 

1427 

.562 

27.43 

1.080 

21000 

10500 

.126 

60 

22225 

1750 

36.512 

1  4375 

8SF114 

19.43 

.765 

16.66 

.656 

31.95 

1.258 

29200 

14600 

.193 

60 

25.400 

1  0000 

41275 

1.6250 

105  Fill 

2222 

175 

19.05 

.750 

36.50 

1437 

38800 

19400 

276 

60 

31.750 

12500 

'  50.800 

2  0000 

I25FL20 

27.76 

1.093 

23.80 

.937 

45.59 

1.795 

61900 

31000 

.516 

6.0 

f 

34,925 

1.3750 

;  55.562 

21875 

135022 

30,15 

1.187 

t- 

26.19 

1.031  i 

4920 

1.937 

74800 

37400 

.770 

55 

'  3*.  100 

1.5000 

f  61.912 

2.4375 

15SFL24 

33.32 

1.312 

2158 

1.125  ! 

:  54,74 

2155 

89300 

44600 

934 

6.0 

TO? 

X  44,450 

1.7500 

1.-71.43* 

28125 

17SFL28 

V-  38.89 

1.531 

3322 

1212  i 

63.8* 

2515 

124000 

62200 

143 

6.0 

G-un.aftL-*50-*00 

2.0000 

pc  80. 962  a 

3.1875 

2BSFU2 

f£44,45  - 

1.750 

K  38.10 

1500 

R}u*‘ 

2.875 

157000 

78400 

207 

6.0 

•»:,  *  •*  . 

...  57.150 

2.2500 

(-5  90,4**  '. 

3.5625 

22SFU6 

L-  50.01 

1.969 

4285 

1.687 

p-O.17 

3.235 

203000 

102000 

2.92 

6.0  " 

Jhrr  \  ^ 

10-500 

2.5000 

•  100.012 

3.9375 

25SO40 

-  55,55 

2117 

'  47.62  ' 

1J75 

(■*91.19 

3  590 

251000 

126000 

4.09 

6.0 

?^TTOr*  1.  “““ 
P-iaA-fcAL-^7*-200 
-  v  *2.550 


2.7500  ;  111.125  4.3750 

3.0000  120.550  it  <  7500 

33500  ;  130,175  '  5.1250 


27SFIM  61,11 
10SF14*  fc.66,6* 
J2SR52  L  7224 


2*06  f.  52.37  2.062  !  100,33. 

2.525  P-' 57,15"*  2250  t<T09.52~ 
2*44  .  61.90  •  2437  I-  118.74 


5  38  6  0 
6.87  6  0 
8.63  6  0 


r  0.900 
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5.5000 

SOU* 
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536000 

261000 
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3  7500  K.  149225  ; 

18750 

J7SO» 

£?  0.34 

32*1  W7L42; 
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5290 
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130 

6.0 

101.600 

4  0000  -,1*158,750  .? 
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405014 
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TITLE  STATIC  ANALYSIS  RESULTS  -  Nitrogen 
SUBTITLE  ENERGY  RECOVERY  05-Feb-87  10:57  AM 

HORIZONTAL_LABEL  THETA  (degrees) 

VERTICAL  LABEL  TORQUE  (ft-lb)  _ 

Y  LEGEND  gravity  —  DVt-’eQL'l* - 

Y_LEGEND  equilibration 
Y_LEGEND  elevation 

Y__LEGEND  dyn  equil  ope<zftT»tW  , 


-  'TOtiQ'jE  _ 
E. 


C  THETA 

TGRAV 

TEQUIL 

TELEV 

DYNEQUIL 

-5.00 

81600.2 

83341.2 

31017.6 

84407.5 

-4.00 

81712.4 

83446.8 

31294.4 

83221.1 

-3.00 

81799.6 

82290.4 

31566.3 

82947.0 

-2.00 

81862.0 

82230.1 

31832.9 

82660.5 

-1.00 

81899.4 

82150.0 

32094.2 

82361.6 

0.00 

81911.9 

82050.2 

32350.1 

82050.4 

1.00 

81899.4 

81930.9 

32600.4 

81727.0 

2.00 

81862.0 

81792.3 

32845.0 

81391.5 

3.00 

81799.6 

81634.5 

33083.7 

81044.1 

4.00 

81712.4 

81457.8 

33316.4 

80684.7 

5.00 

81600.2 

81262.2 

33542.9 

79092.1 

6.00 

81463.2 

81047.9 

33763.1 

78428.6 

7.00 

81301.3 

80815.1 

33976.8 

77762.5 

8.00 

81114.7 

80563.8 

34183.8 

77093.5 

9.00 

80903.4 

80294.2 

34384.0 

76421.6 

10.00 

80667.5 

80006.4 

34577.1 

75746.8 

11.00 

80406.9 

78524.7 

34763.1 

74297.9 

12.00 

80121.9 

78087.4 

34941.7 

73265.7 

13.00 

79812.5 

77634.4 

35112.7 

72245.1 

14.00 

79478.8 

77165.9 

35276.0 

71235.5 

15.00 

79120.8 

76682.0 

35431.2 

70236.4 

16.00 

78738.8 

76183.1 

35578.4 

69247.4 

17.00 

78332.7 

75669.1 

35717.1 

68267.8 

18.00 

77902.8 

75140.3 

35847.2 

67297.3 

19.00 

77449.2 

74596.8 

35968.5 

64793.0 

20.00 

76972.0 

74038.9 

36080.8 

63443.7 

21.00 

76471.3 

73466.5 

36183.8 

62123.9 

22.00 

75947.4 

72879.9 

36277.2 

60832.4 

23.00 

75400.3 

71393.6 

36360.9 

58282.3 

24.00 

74830.2 

70680.7 

36434.5 

56485.1 

25.00 

74237.4 

69956.8 

36497.9 

54751.1 

26.00 

73621.9 

69222.1 

36550.7 

53077.1 

27.00 

72984.0 

68476.8 

36592.6 

49626.8 

28.00 

72323.9 

67721.0 

36623.5 

47225.4 

29.00 

71641.7 

66955.0 

36643.0 

48384.0 

30.00 

70937.8 

66178.8 

36650.9 

46920.3 

31.00 

70212.2 

65392.7 

36646.7 

45502.7 

32.00 

69465.2 

64596.8 

36630.3 

44129.2 

33.00 

68697.1 

63791.1 

36601.3 

42797.7 

34.00 

67908.0 

62976.0 

36559.5 

44239.7 

35.00 

67098.3 

62151.4 

36504.4 

43722.5 

36.00 

66268.1 

61317.6 

36435.8 

43197.0 

37.00 

65417.7 

59740.1 

36353.3 

42663.0 

38.00 

64547.4 

58816.1 

36256.6 

42120.7 

39.00 

63657.5 

57886.6 

36145.4 

41570.1 

40.00 

62748.1 

56951.6 

36019.3 

41011.2 

41.00 

61819.7 

56011.4 

35877.9 

40444.1 

42.00 

60872.4 

55066.0 

35721.0 

39868.9 

43.00 

59906.5 

54115.7 

35548.1 

39285.6 

44.00 

58922.5 

53160.4 

35358.9 

38694.1 

45.00 

57920.4 

52200.3 

35153.1 

38094.7 

46.00 

56900.8 

51235.7 

34930.3 

37487.3 

y.v.y* 

'  ■»  .  'L*-  *.■%  'll  *  •**  »  - 

47.00 

55863.8 

50266.5 

34690.1 

36872.0 

48.00 

54809.7 

49292.8 

34432.2 

36248.9 

49.00 

53739.0 

48314.9 

34156.1 

35618.0 

50.00 

52651.9 

47332.7 

33861.7 

34979.4 

51.00 

51548.8 

46346.4 

33548.5 

34333.2 

52.00 

50430.0 

44849.9 

33216.1 

33679.4 

53.00 

49295.8 

43814.3 

32864.3 

32527.3 

54.00 

48146.6 

42778.2 

32492.7 

31814.6 

55.00 

46982.7 

41741.5 

32101.1 

31098.5 

56.00 

45804.5 

40704.3 

31689.0 

30378.9 

57.00 

44612.4 

39666.8 

31256.3 

29656.0 

58.00 

43406.7 

38629.0 

30802.6 

28929.7 

59.00 

42187.7 

37590.9 

30327.7 

28200.0 

60.00 

40955.9 

36552.6 

29831.3 

27466.9 

61.00 

39711.6 

35514.1 

29313.4 

26730.6 

62.00 

38455.3 

34475.5 

28773.6 

25990.9 

63.00 

37187.2 

33436.8 

28211.8 

25248.0 

64.00 

35907.8 

32398.0 

27627.9 

24501.9 

65.00 

34617.4 

31359.2 

27021.8 

23752.6 

66.00 

33316.5 

30320.5 

26393.5 

23000.2 

67.00 

32005.5 

29281.7 

25742.9 

22244.8 

68.00 

30684.7 

28243.0 

25070.1 

21486.3 

69.00 

29354.5 

27204.4 

24375.1 

20724.9 

70.00 

28015.5 

26165.8 

23658.1 

19960.6 

71.00 

26667.8 

25127.4 

22919.2 

19193.6 

72.00 

25312.1 

24089.1 

22158.6 

18072.7 

c 


OJORttUb  PK€*SUKfc  -  3000  PSI 


1 

1 

1 


37 


K 


Appendix  C 


DADS  Model 


fi' 


Includes  a  description  of  the  DADS  software, 
the  input  data  for  the  aodel  and 
the  user  written  suppleaentary  code  for  the  aodel. 


3? 


Introduction  to  DAOS  Software 


Purpose  and  Use  of  DADS 

The  Dynamic  Analysis  and  Design  System  software  Is  a  set  of  general 
purpose  computer  programs  that  can  be  used  to  model  and  predict  the  motion 
of  a  variety  of  real  world  mechanical  systems.  Using  a  set  of  data  that 
describes  the  machine  to  be  modeled,  DAOS  builds  a  mathematical  model  of 
the  real  system  that  calculates  positions,  velocities,  and  accelerations 
of  the  various  parts  of  the  machine,  as  well  as  resultant  forces  that  act 
In  the  system.  By  using  such  a  computer  program  to  analyze  a  machine,  the 
designer  can  simulate  the  behavior  of  a  wide  range  of  alternate  designs 
prior  to  building  and  testing  prototypes. 


Types  of  Machines  That  Can  Be  Modeled 

DADS  contains  a  large  library  of  mechanical  elements  that  can  be  used 
to  build  a  model.  These  Include  rigid  and  flexible  bodies,  joints  and 
other  constraints,  force-  and  torque-producing  elements,  as  well  as 
control  and  hydraulic  elements.  Models  can  be  created  in  two  or  three 
dimensions.  Any  machine  whose  motion  is  entirely  planar  can  be  modeled 
more  easily  in  2D  and  can  be  analyzed  much  more  quickly  by  the  analysis 
program,  because  of  the  smaller  number  of  variables  required.  The  3D 
version  is  used  to  model  any  machine  whose  motion  Is  not  contained  within 
a  plane,  allowing  complete  generality.  If  desired,  planar  systems  can 
also  be  modeled  In  3D. 

The  control  elements  can  be  used  .In  any  rigid  or  flexible  body  model 
that  Is  analyzed  using  the  dynamic  option,  or  can  be  used  Independently  to 
create  a  control  system  model  without  any  bodies  in  It.  The  control 
elements  can  be  used  to  apply  forces,  torques,  or  constraints  to  the  rigid 
or  flexible  body  system.  The  control  elements  define  a  system  of  first 
order  differential  equations  that  are  solved  along  with  the  set  of  second 
order  differential  equations  associated  with  the  rigid  or  flexible  body 
equations  of  motion. 

The  ability  to  model  .the  flexibility  effect  allows  much  more  detailed 
mechanical  models  to  be  built.  Data  for  flexible  bodies  are  generated  by 
using  the  DADS  Intermediate  Processor  to  process  data  from  a  finite 
element  program.  Any  elements  In  the  finite  element  model  can  be  used, 
not  just  simple  beam  elements  between  nodes.  DADS  uses  modal  coordinates, 
which  are  synthesized  from  the  finite  element  data,  to  represent  the 
flexibility  effects  of  all  flexible  bodies  In  the  model.  This  approach 
allows  the  most  generality  for  modeling  flexible  bodies  within  the 
mechanical  system. 


Physical  data  entered  into  model 

Systea  geometry: 

joint  locations 
CG  locations 

Inertial  properties: 
body  weights 

rotational  inertia  about  CG 

Descriptions  of  hydraulic  coaponents 

accumulator  volume  and  charge  pressure 
control  valve  size 
amount  of  lapping  in  the  valve 
actuator  size 

trapped  volumes  on  either  side  of  the  piston 
area  of  each  side  of  the  piston 
connectivity  of  hydraulic  circuit 

Hydraulic  systea  parameters 
oil  bulk  modulus 
supply  pressure 
return  line  pressure 
valve  position 

Miscellaneous 

gravity  vector 

initial  values  of  state  variables 
error  tolerances  for  simulation 
initial  integration  time  step 
integration  technique 


Output  available  from  the  model 


Forces 

reaction  forces  at  joints 
forces  acting  on  body  CGs 
forces  produced  by  actuators 

Kinematic  information 

body  position  and  orientation 

linear  and  angular  velocity  of  bodies 

linear  and  angular  acceleration  of  bodies 

Hydraulic  systea  information 

actuator  pressures  and  pressure  rates 
actuator  motion 

state  variables  and  their  derivatives 


Input  data  for  the  DADS  aodel  of  the  traversal  system 


CREATE  HEADER 


TRAVERSAL  SYSTEM 

This  is  a  model  of  the  LTHD  traversal  system.  It  contains  the  main 
valve,  the  traversal  cylinder,  energy  accumulator,  the  gun  body 
the  friction  in  the  gimb  1  bearing  and  locking  mechanisms  and  a  model 
of  the  man  which  is  a  first  order  system  representing  the  delay  due  to 
his  reaction  time. 


ANALYSIS 

CREATE  SYSTEM. DATA 
UNITS 

ANALYSIS. TYPE 
STARTING. TIME 
ENDING. TIME 
PRINT. INTERVAL 
GRAVITY . SEA . LEVEL 

X.  GRAVITY 

Y.  GRAVITY 

Z.  GRAVITY 

SCALE. GRAVITY. COEF 
MATRIX . OPERATIONS 
REDUNDANCY. CHECK 
LU . TOL 

ASSEMBLY. TOL 
BYPASS. ASSEMBLY 
OUTPUT. FILE 
REFERENCE. FRAME 
DEBUG. FLAG 
UP 

CREATE  DYNAMIC. DATA 
REACTION. FORCES 
FORCE. COORDINATES 
PRINT. METHOD 
MAX. INT. STEP 
SOLUTION. TOL 
INTEGRATION. TOL 
UP 
UP 

JOINTS 


' INCH-LB-SEC' 
'DYNAMIC' 
'0.0' 

'3' 

'  .01' 

'  386  ' 

'0.0' 

'0.0' 

'-1.0' 

'1.0' 

'SPARSE' 

'TRUE' 

' 1 . OD-12  ' 

'  1 . OD-3 ' 

'FALSE' 

'BINARY' 

'GLOBAL' 

'TRUE' 


'FALSE' 
'GLOBAL' 
'ACTUAL' 
'  .01' 
'0.001' 
'0.0001' 


CREATE  REVOLUTE. JOINT 
NAME 

BODY. 1 .NAME 
BODY. 2. NAME 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 

Q.  ON. BODY. 2 

R.  ON. BODY. 1 
R. ON. BODY. 2 
NODE . 1 
NODE.  2 

UP 

UP 

CONTROLS 

CREATE  ACCUMULATOR 
NAME 

INPUT. NODE 
GAMMA 

INITIAL. PRESSURE 
VOLUME 

CHARGE. PRESSURE 
ATMOS. PRESSURE 


-  ' GIMBAL_REV' 

-  ' GND  BOD' 


-  'GUN  BOD' 


( 

0/ 

"-16, 

0 

) 

( 

0, 

-16, 

0 

) 

( 

0, 

-16, 

1 

) 

( 

0, 

-16, 

1 

) 

( 

1, 

-16, 

0 

) 

( 

1, 

-16, 

0 

) 

-  '0' 
-  >  o  • 


-  ' ENERGY_ACCUMULATOR 

-  'ACC  PSI' 

-  '  1 . 4t 
.  '3000' 

.  '4500' 

-  '2500' 

-  '14.7' 


UP 

CREATE  AMPLIFIER 
NAME 

INPUT. NODE 
OUTPUT. NODE 
TYPE 
GAIN 

CURVE. NAME 
UP 

CREATE  DOUBLE. ACTUATOR 
NAME 

INPUT. NODE. 1 
INPUT. NODE. 2 
OUTPUT. NODE 
BODY. 1 .NAME 
BODY. 2. NAME 
PRESSURE. 1 
PRESSURE. 2 
AREA.  1 
AREA. 2 
VOLUME. 1 
VOLUME. 2 
COULOMB. COEFF 
VISCOUS. COEFF 
P. ON. BODY. 1 
P. ON. BODY. 2 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
UP 

CREATE  FIRST. ORDER 
NAME 

INPUT. NODE 
OUTPUT. NODE 
D. FIRST. ORDER. COEF 
D. CONSTANT. TERM 
N. FIRST. ORDER. COEF 
N. CONSTANT. TERM 
S . ZERO 
UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE . NAME 
TYPE 

BODY . 1 . NAME 
BODY . 2 . NAME 
FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE . NAME 
JOINT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UN ITS 

UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY. 1. NAME 
BODY. 2. NAME 


-  ' AMP_SWITCH' 

-  ' SWITCH_VALUE '  \ 

-  ' AMP_OUTPUT ' 

-  'CONSTANT' 

-  ' .01' 

«*  'NONE' 


-  ' TRAVERS E_CYL' 

-  ' EXT_PSI ' 

-  ' RET_PSI ' 

-  ' CYL_FORCE ' 

-  ' GND_BOD ' 

=  'GUN  BOD' 

-  'ISO'S' 

-  '2000' 

-  '14.186' 

-  '9.276' 

-  '85.118' 

-  '49.708' 

-  '0.0' 

-  '10' 

-  (  27.5,  -28.938,  0  ) 

-  (  0,  -30,  0  ) 

-  '0' 

-  '0' 


-  ' MAN_MODEL ' 

-  ' GUN_POSITION' 

-  ' MAN_OUTPUT ' 

-  '  .3' 

-  '1' 

.  '0' 

-  '1.0' 

-  '0.0' 


-  'VELOCITYJFEEDBACK' 

-  'GUN_VELOCITY' 

-  'ZL. OMEGA' 

-  'GUN  BOD' 


'NONE' 

(  0.0, 

0.0,  0.0, 

(  o.o. 

0.0,  0.0 

) 

(  0.0, 

0.0,  0.0 

) 

(  1.0, 

0.0,  0.0 

) 

(  1.0, 

0.0,  0.0 

) 

'NONE' 

'NONE' 

-  '0' 

-  '0' 

-  'DEGREES' 


-  'TANK  PRESSURE' 

-  ' TANK~PS I ' 

-  'POLYNOMIAL' 

-  'NONE' 

-  ’NONE' 


1 


I 

B 


to 


I 

1 

s 

a 

A 


i 

I 


FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JOINT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UNITS 

UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY. 1 .NAME 
BODY. 2 .NAME 
FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JOINT. NAME 
FLEXIBLE. NODE.  1 
FLEXIBLE. NODE. 2 
ANGULAR. UNITS 

UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY . 1 . NAME 
BODY. 2 .NAME 
FUNCTION. PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JO I NT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UN ITS 

UP 

CREATE  INTEGRATOR 
NAME 

INPUT. NODE 
OUTPUT. NODE 
S . ZERO . 

UP 

CREATE  OUTPUT 
NAME 

OUTPUT. NODE 
TYPE 

BODY. 1. NAME 
BODY. 2 .NAME 
P. ON. BODY. 1 
P. ON. BODY. 2 
JO I NT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 


(  300,  0,  0,  0  ) 

(  0.0,  0.0,  0.0  ) 

(  0.0,  0.0,  0.0  ) 

(  1.0,  0.0,  0.0  ) 

(  1.0,  0.0,  0.0  ) 

'NONE' 

'NONE' 

'0' 

'0' 

'DEGREES' 


' VALVE_INPUT' 

' VALVE_INPUT ' 
'STEP' 

'NONE' 

'NONE' 

(  .5,  .01,  0,  0  ) 
(  0.0,  0.0,  0.0  ) 

(  0.0,  0.0,  0.0  ) 

(  1.0,  0.0,  0.0  ) 

(  1.0,  0.0,  0.0  ) 

'NONE' 

'NONE' 

'0' 

'0' 

'DEGREES' 


' DUMMY_INPUT ' 

' DUMMY_OUTPUT ' 
'STEP' 

'NONE' 

'NONE' 

(  0.0,  0.0,  0 
(  0.0,  0.0,  0 
(  0.0,  0.0,  0 
(  1.0,  0.0,  0 


0,  0.0  ) 


1.0,  0.0,  0.0 

NONE' 

'NONE' 

■0' 

’  0  ' 

'DEGREES' 


'  POS ITION_FEEDBACK ' 

'GUN_VELOCITY' 

'GUN_POSITION' 

'0.0' 


’ DUMMY_OUTPUT ' 

' DUMMY_OUTPUT ' 

' FORCE. DIFF' 

' GND_BOD ' 

'GUN  BOD' 

(  2775,  -28.938,  0  ) 
(  0,  -30,  0  ) 

'NONE' 

'0' 

'0' 


CREATE  PARAMETER 
NAME 

BULK. MODULUS 
UP 

CREATE  SUMMER 
NAME 

OUTPUT. NODE 
INPUT. NODE. 1 
INPUT. NODE. 2 
INPUT. NODE. 3 
COEFFICIENT. 1 
COEFFICIENT. 2 
COEFFICIENT. 3 
UP 

CREATE  SWITCH 
NAME 

INPUT. NODE 
OUTPUT. NODE 
ON. VALUE 
UP 

CREATE  VALVE 
NAME 

INPUT. NODE 
TYPE 

CHAMBER. 1 .NODE 
CHAMBER. 2. NODE 
CHAMBER. 3 .NODE 
CHAMBER. 4 .NODE 
FLOW. COEFFICIENT 
MAX. SPOOL 
LAP. SPOOL 
UP 
UP 

CREATE  BODY 
NAME 

CENTER . OF . GRAVITY 
TYPE . ANGULAR . COORD 
ANGLE. 1 
ANGLE. 2 
ANGLE. 3 

FIXED. TO. GROUND 
MASS 

INERTIA. XXL 
INERTIA. YYL 
INERTIA. ZZL 
INERTIA. XYL 
INERTIA. XZL 
INERTIA. YZL 
XG. FORCE 
YG. FORCE 
ZG. FORCE 
XL. TORQUE 
YL. TORQUE 
ZL. TORQUE 
CURVE. XGF 
CURVE. YGF 
CURVE. ZGF 
CURVE. XLT 
CURVE. YLT 
CURVE. ZLT 
SIGN. EO 


BULK  MODULUS' 

1750^0' 


'SUMMER' 

' VALVE_POS ITION ' 
' VALVE_ INPUT ' 

' AMP_OUTPUT ' 
'NONE' 


' SWITCH' 

' MAN_OUTPUT ' 

' SWITCH_VALUE ' 

' .0005' 


' MAIN_VALVE ' 

' VALVE_POS ITION ' 
' 4WAY ' 

' ACC_PS I ' 

' EXT_PS I ' 

' RET_PS I ' 

' TANK_PS I ' 

'  .4' 

'1' 

'0.0' 


' GND  BOD' 

(  O.ff,  0.0,  0.0  ) 

'EULER' 

'0.0' 

'0.0' 

'0.0' 

'TRUE' 


■0.0' 

’NONE' 

’  NONE' 
’NONE' 

’  NONE ' 
NONE' 
’NONE' 
’POSITIVE1 
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ANGULAR. UNITS 

FLEXIBLE 

SUPERELEMENT 

UP 

CREATE  BODY 
NAME 

CENTER . OF . GRAVITY 
TYPE .ANGULAR . COORD 
ANGLE. 1 
ANGLE. 2 
ANGLE. 3 

FIXED. TO. GROUND 
MASS 

INERTIA. XXL 
INERTIA. YYL 
INERTIA. ZZL 
INERTIA. XYL 
INERTIA. XZL 
INERTIA. YZL 
XG. FORCE 
YG. FORCE 
ZG. FORCE 
XL. TORQUE 
YL. TORQUE 
ZL. TORQUE 
CURVE. XGF 
CURVE. YGF 
CURVE. ZGF 
CURVE. XLT 
CURVE. YLT 
CURVE. ZLT 
SIGN.EO 
ANGULAR. UNITS 
FLEXIBLE 
SUPERELEMENT 

UP 

CREATE  INITIAL. CONDITION 
NAME 

BODY. 1 .NAME 
BODY. 2 .NAME 
ELEMENT. NAME 
TYPE. INITIAL. COND 
INITIAL. VALUE 
TIME. DERIVATIVE 
OMEGA.  Y 
OMEGA.  Z 
P. ON. BODY. 1 
P. ON. BODY. 2 
EXTRA. COORD 
ANGULAR. UNITS 


-  'DEGREES' 

-  'FALSE' 

-  'FALSE' 


-  ' GUN_BOD' 

-  (  0,  177.3,  0  ) 

-  'EULER' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  'FALSE' 

-  '16.34' 

«  '1.0' 

-  '1.0' 

-  '69250' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

=  '0.0' 

=  '0.0' 

-  'NONE' 

=  'NONE' 

=  'NONE' 

=  'NONE' 

*  'NONE' 

*  'NONE' 

=  'POSITIVE' 

-  'DEGREES' 

-  'FALSE' 

-  'FALSE' 


-  ' GUN_X_IC ' 

-  ' GUN_BOD ' 

-  'NONE' 

-  'NONE' 

-  'X' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  '0.0' 

-  (  0.0,  0.0,  0.0  ) 
-  (  0.0,  0.0,  0.0  ) 
-  '0' 

-  'DEGREES' 


UP 


CREATE  HEADER 


ELEVATION  SYSTEM 

is  a  model  of  the  LTHD  elevation  system.  It  contains  the  main 
,  the  elevation  cylinder,  the  energy  accumulator,  the  gun  body  the 
in  the  trunnion  bearings  and  locking  mechanism,  a  constant 

represents  the  inbalance  between  equilibration  and  gravity 
of  elevation  and  a  model  of  the  man  which  is  a  first  order 
represents  the  delay  due  to  his  reaction  time. 


This 
valve 
friction 
torque  which 
at  a  point 
system  that 

ANALYSIS 
CREATE  SYSTEM. DATA 
UNITS 

ANALYSIS. TYPE 
STARTING. TIME 
ENDING. TIME 
PRINT. INTERVAL 
GRAVITY . SEA . LEVEL 

X.  GRAVITY 

Y.  GRAVITY 

Z.  GRAVITY 

SCALE. GRAVITY. COE F 
MATRIX. OPERATIONS 
REDUNDANCY. CHECK 
LU.TOL 

ASSEMBLY. TOL 
BYPASS .ASSEMBLY 
OUTPUT. FILE 
REFERENCE. FRAME 
DEBUG. FLAG 
UP 

CREATE  DYNAMIC. DATA 
REACTION. FORCES 
FORCE. COORDINATES 
PRINT. METHOD 
MAX. INT. STEP 
SOLUTION. TOL 
INTEGRATION. TOL 
UP 
UP 

JOINTS 

CREATE  REVOLUTE. JOINT 
NAME 

BODY. 1 .NAME 
BODY. 2 .NAME 

P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 

Q.  ON. BODY. 2 

R.  ON. BODY. 1 

R . ON . BODY . 2 
NODE . 1 
NODE. 2 

UP 

UP 

CONTROLS 

CREATE  ACCUMULATOR 
NAME 

INPUT. NODE 
GAMMA 

INITIAL. PRESSURE 
VOLUME 

CHARGE. PRESSURE 


-  'INCH-LB-SEC' 

-  'DYNAMIC' 

.  '0.0' 

-  '3' 

-  '  .01' 

-  '386' 

-  '1' 

-  '0.0' 

-  '0' 

-  '1.0' 

-  'SPARSE' 

*  'TRUE' 

*  ' 1 . OD-12 ' 

-  ' 1 . 0D-3 ' 

»  'FALSE' 

-  'BINARY' 

-  'GLOBAL' 

»  'TRUE' 


-  'FALSE' 

-  'GLOBAL' 

-  'ACTUAL' 

*  '  .01' 

-  '0.001' 

«  '0.0001' 


' TRUNNION’ 
' GND_BOD ' 
'GUN  BOD' 


(  o,' 

“0, 

0  ) 

(  0, 

0, 

0  ) 

(  1, 

0, 

0  ) 

(  1, 

0, 

0  ) 

(  0, 

1, 

0  ) 

(  0, 
'0' 

'  0' 

1, 

0  ) 

' ENERGY_ 

ACCUMULATOR 

'ACC 

PSI 

’  f 

-  '  1 . 4~ 

-  '3000' 
.  '4500' 

-  '2500' 


ATMOS. PRESSURE 


14.7 


1 

fig 


UP 

CREATE  AMPLIFIER 
NAME 

INPUT. NODE 
OUTPUT. NODE 
TYPE 
GAIN 

CURVE. NAME 
UP 

CREATE  DOUBLE. ACTUATOR 
NAME 

INPUT. NODE. 1 
INPUT. NODE. 2 
OUTPUT. NODE 
BODY. 1 .NAME 
BODY. 2. NAME 
PRESSURE. 1 
PRESSURE. 2 
AREA. 1 
AREA. 2 
VOLUME. 1 
VOLUME . 2 
COULOMB. COEFF 
VISCOUS. COEFF 
P. ON. BODY. 1 
P. ON. BODY. 2 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
UP 

CREATE  FIRST. ORDER 
NAME 

INPUT. NODE 
OUTPUT. NODE 
D. FIRST. ORDER. COEF 
D. CONSTANT. TERM 
N. FIRST. ORDER. COEF 
N. CONSTANT. TERM 
S . ZERO 


'  AMP_SWITCH' 

' SWI TCH_VALUE ' 
' AMP_OUTPUT ' 
'CONSTANT' 

'  .01' 

'NONE' 


' ELEVATION_CYL ' 
' EXT_PSI ' 

' RET_PSI ' 

'  CYL_FORCE ' 

’ GND_BOD ' 

'GUN  BOD' 

'1007' 

'2237.6' 

'6.492' 

'4.725' 

'289.7' 

'23.3' 

'0.0' 

'10' 

(  0,  -10,  35  ) 

(  0,  96,  18.5  ) 
'0' 

'0' 


'  MAN_MODEL ' 
'GUN_POSITION' 
'  MAN_OUTPUT ' 

'  .3' 

'1' 

'0' 

'1.0' 

'0.0' 


UP  # 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY. 1 .NAME 
BODY. 2 .NAME 
FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JOINT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UNITS 


.  ' VELOCITY_FEEDBACK ' 

-  ' GUN_VELOCITY' 

-  'XL. OMEGA' 

-  'GUN  BOD' 


'NONE' 

(  0.0, 

0.0, 

0.0, 

(  0.0, 

0.0, 

0.0  ) 

(  0.0, 

0.0, 

0.0  ) 

(  1.0, 

0.0, 

0.0  ) 

(  1.0, 
'NONE' 

'  NONE ' 

0.0, 

0.0  ) 

.  '0' 

-  '0' 

-  'DEGREES' 


UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY. 1 .NAME 


'  TANK_PRESSURE ' 
'  TANK_PSI ' 
'POLYNOMIAL' 
'NONE' 


i 


1 


I 

1 

I 

1 

8 


1 

I 

ff 


8 

I 

I 


BODY. 2. NAME 
FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JO I NT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UN ITS 

UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY. 1 .NAME 
BODY. 2 .NAME 
FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JO I NT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UNITS 

UP 

CREATE  INPUT. FUNCTION 
NAME 

NODE. NAME 
TYPE 

BODY. 1 .NAME 
BODY. 2 .NAME 
FUNCTION . PARAMETERS 
P. ON. BODY. 1 

P.  ON. BODY. 2 

Q.  ON. BODY. 1 
Q. ON. BODY. 2 
CURVE. NAME 
JOINT. NAME 
FLEXIBLE. NODE. 1 
FLEXIBLE. NODE. 2 
ANGULAR. UNITS 

UP 

CREATE  INTEGRATOR 
NAME 

INPUT. NODE 
OUTPUT. NODE 
S . ZERO 
UP 

CREATE  OUTPUT 
NAME 

OUTPUT. NODE 
TYPE 

BODY. 1. NAME 
BODY. 2. NAME 
P. ON. BODY. 1 
P. ON. BODY. 2 
JOINT. NAME 
FLEXIBLE. NODE. 1 


-  'NONE' 

-  (  300,  0,  0,  0  ) 

-  (  0.0,  0.0,  0.0  ) 
-  (  0.0,  0.0,  0.0  ) 
-  (  1.0,  0.0,  0.0  ) 
-  (  1.0,  0.0,  0.0  ) 

-  'NONE' 

-  'NONE' 

-  '0' 

-  '0' 

-  'DEGREES' 


'VALVE 

INPUT' 

'  VALVE' 

“INPUT' 

' STEP' 

'NONE' 

'  NONE ' 

(  .5, 

o 

o 

rH 

o 

) 

{  0.0, 

0.0,  0.0 

) 

(  0.0, 

0.0,  0.0 

) 

(  i.o, 

0.0,  0.0 

) 

(  1.0, 

0.0,  0.0 

) 

'NONE' 

'NONE' 

>  » 0 ' 

-  '0' 

-  'DEGREES' 


'DUMMY_ INPUT' 

' DUMMY_OUTPUT ' 
'STEP' 

'NONE' 


'NONE' 
(  0.0, 
(  0.0, 
(  0.0, 
(  1.0, 
(  1.0, 
'NONE' 


0.0,  0.0  ) 

0.0,  0.0  ) 

0.0,  0.0  ) 

0.0,  0.0  ) 


'NONE' 


'0' 

'0' 

'DEGREES' 


-  ' POSITION_FEEDBACK ' 

-  ' GUN_VELOCITY ' 

-  'GUN  POSITION' 

-  '0.0T 


-  ' DUMMY_OUTPUT ' 

-  ' DUMMY_OUTPUT ' 

-  'FORCE. DIFF' 

-  ' GND  BOD' 

-  'GUN~BOD' 

-  (  0,~-10,  35  ) 

-  (  0,  96,  18.5  ) 

-  'NONE' 

•  '0' 


s 


50 


B 

I 


I 

B 

1 


I 


FLEXIBLE. NODE. 2 
UP 

CREATE  PARAMETER 
NAME 

BULK. MODULUS 
UP 

CREATE  SUMMER 
NAME 

OUTPUT. NODE 
INPUT. NODE. 1 
INPUT. NODE. 2 
INPUT. NODE. 3 
COEFFICIENT. 1 
COEFFICIENT . 2 
COEFFICIENT . 3 
UP 

CREATE  SWITCH 
NAME 

INPUT. NODE 
OUTPUT. NODE 
ON. VALUE 
UP 

CREATE  VALVE 
NAME 

INPUT. NODE 
TYPE 

CHAMBER. 1. NODE 
CHAMBER. 2. NODE 
CHAMBER. 3. NODE 
CHAMBER. 4 .NODE 
FLOW. COEFFICIENT 
MAX. SPOOL 
LAP. SPOOL 
UP 
UP 

CREATE  BODY 
NAME 

CENTER. OF. GRAVITY 
TYPE . ANGULAR . COORD 
ANGLE.  1 
ANGLE . 2 
ANGLE.  3 

FIXED. TO. GROUND 
MASS 

INERTIA. XXL 
INERTIA. YYL 
INERTIA. ZZL 
INERTIA. XYL 
INERTIA. XZL 
INERTIA. YZL 
XG. FORCE 
YG. FORCE 
ZG. FORCE 
XL. TORQUE 
YL. TORQUE 
ZL. TORQUE 
CURVE. XGF 
CURVE . YGF 
CURVE. ZGF 
CURVE. XLT 
CURVE. YLT 
CURVE. ZLT 


'0' 


'BULK  MODULUS' 
'  1750^0 ' 


'SUMMER' 

'  VALVE_POSITION ' 
'  VALVE_INPUT ' 

'  AMP_OUTPUT ' 
'NONE' 

'  +  ' 


'  +  ' 


'SWITCH'  ' 

'  MAN_OUTPUT ' 

'  SWI TCH_VALUE ' 
' .0005' 


'  MAIN_VALVE ' 

' VALVE_POSITION' 
'  4WAY ' 

'  ACC_PSI ' 

'  EXT_PSI ' 

'  RET_PSI ' 

'  TANK_PSI ' 

'  .  4  ' 

'1' 

'0.0' 


' GND  BOD' 

(  0.75,  0.0,  0.0  ) 

'EULER' 

'0.0' 

'0.0' 

'0.0' 

'TRUE' 

'1.0' 

'1.0' 

'1.0' 

'1.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'  NONE  ’ 

'NONE' 

'  NONE  ’ 

'  NONE ' 

'  NONE ' 

'NONE' 


I 

i 

I 

I 

I 

1 


I 

B 

I 

8 

I 


SIGN. EO 
ANGULAR. UNITS 
FLEXIBLE 
SUPERELEMENT 

UP 

CREATE  BODY 
NAME 

CENTER. OF. GRAVITY 
TYPE . ANGULAR . COORD 
ANGLE. 1 
ANGLE. 2 
ANGLE. 3 

FIXED. TO. GROUND 
MASS 

INERTIA. XXL 
INERTIA. YYL 
INERTIA. ZZL 
INERTIA. XYL 
INERTIA. XZL 
INERTIA. YZL 
XG. FORCE 
YG. FORCE 
ZG. FORCE 
XL. TORQUE 
YL. TORQUE 
ZL. TORQUE 
CURVE. XGF 
CURVE. YGF 
CURVE. ZGF 
CURVE. XLT 
CURVE. YLT 
CURVE. ZLT 
SIGN. EO 
ANGULAR. UN ITS 
FLEXIBLE 
SUPERELEMENT 

UP 

CREATE  INITIAL. CONDITION 
NAME 

BODY. 1 .NAME 
BODY. 2. NAME 
ELEMENT. NAME 
TYPE. INITIAL. COND 
INITIAL. VALUE 
TIME. DERIVATIVE 
OMEGA. Y 
OMEGA.  Z 
P. ON. BODY. 1 
P. ON. BODY. 2 
EXTRA. COORD 
ANGULAR. UNITS 


'POSITIVE' 

'DEGREES' 

'FALSE' 

'FALSE' 


'GUN  BOD' 

(  0, “161. 75,  0  ) 
•EULER' 

'0.0' 

'0.0' 

'0.0' 

'FALSE' 

'15.74' 

'83700' 

'1.0' 

'1' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

'120000' 

'0.0' 

'0.0' 

•NONE' 

'NONE' 

'NONE' 

'NONE' 

' NONE ' 

'NONE' 

'POSITIVE' 

'DEGREES' 

' FALSE' 

' FALSE' 


'GUN  Z_IC' 
'GUN_BOD' 

'NONE' 

'NONE' 

'Z' 

'0.0' 

'0.0' 

'0.0' 

'0.0' 

(  0,  10,  0  ) 

(  0.0,  0.0,  0.0  ) 

'0' 

'DEGREES' 


1 


I 


nnonnnonnnn 


The  following  routine  calculates  forces  defined  by  the 
DADS  user.  The  top  part  of  the  routine  is  the  basic  template 
provided  to  the  user,  the  bottom  part  is  the  code  necessary 
for  the  friction  calculations 

'----------Author  and  Modification  Section-— —————————— 

Author:  Chuck  Mead 

Date  written:  December  5,  1985 
Modifications : 

Copyright  (c)  CADSI  1985 


C-»— — — — — Algo  r  i  thm  Descr  ipt ion-—-——— ——————————— 

C 

C  Purpose  and  use: 

C  This  subroutine  is  used  to  calculate  the  components  of  any  user- 

C  defined  forces  defined  in  the  model.  Any  of  the  various  input 

C  variables  can  be  used  in  these  calculations ,  with  the  results 

C  being  added  to  the  generalized  coordinate  force  array,  FRC .  In 
C  addition,  the  subroutine  SPL32  can  be  used  to  reference  curve 

C  data  elements  for  use  in  the  calculations. 

C  The  various  input  arguments  are  described  and  declared  in  three 

C  separate  groups:  read-only  arguments,  arguments  read  and  modi- 
C  fied,  and  return-only  arguments. 

C  Read-only  arguments  are  passed  in  with  meaningful  values  which 

C  should  not  be  changed  within  this  subroutine.  This  constitutes 
C  the  largest  group  of  arguments  passed.  Consult  the  documentation 
C  on  user-defined  force  elements  for  a  description  of  the  structure 

C  and  contents  of  the  various  variables  and  arrays  passed  in. 

C  Arguments  read  and  modified  are  also  passed  in  with  meaningful 

C  values  which  may  be  changed  within  this  subroutine.  This  group 

C  includes  the  user-defined  force  element  data  arrays,  UDF  and  IUDF, 

C  and  the  two  generalized  coordinate  force  arrays,  FRC  and  XFRC. 

C  The  final  group  of  arguments,  return-only  arguments,  are  passed 

C  in  with  unknown  values  and  should  therefore  be  assigned  some  value 
C  within  this  subroutine  before  they  are  used  in  any  other  computa- 

C  tion.  The  only  variable  in  this  group  is  ERRCOD.  Since  this 

C  variable's  value  is  used  in  the  calling  routine  to  determine  the 

C  success  of  the  computations  performed  in  FRC32 ,  it  is  suggested 

C  that  ERRCOD  be  initialized  to  zero  before  beginning  any  other 
C  computations. 

C 

C  Error  conditions: 

C  Error  conditions  are  reported  via  the  variable  ERRCOD.  A  zero 

C  value  indicates  no  error;  a  positive  value  indicates  some  fatal 
C  error,  i.e.,  files  will  be  closed  and  execution  aborted  upon 
C  return  to  the  calling  routine;  negative  values  are  reserved  for 

C  future  use. 

C 

C  Machine  dependencies:  none  known 
C 


c 


SUBROUTINE  FRC32  (  UDF ,  IUDF,  NUDF ,  NPTRS ,  MPTRS ,  TIME, 


4  Q,  QD,  QDD,  XQ,  XQD,  XQDD,  NXCRD,  FRC , 

4  XFRC,  RB,  NRB ,  NPRB ,  NTYPES ,  NELEMS, 

4  A,  NMAX,  IA,  MMAX,  IPRDCT,  INFOF ,  ERRCOD  ) 

C«-»««------<«Variable  Descr ipt ions 

C 

C - Read-only  arguments - 

C 

C  NUDF . Number  of  user-defined  force  elements  in  the  input  data. 

C  (This  variable  has  a  minimum  value  of  1  for  array  dec- 

C  laration  purposes.) 

C  NPTRS ... .Number  of  real  (double  precision)  data  stored  for  each 
C  user-defined  force  element.  (This  variable  has  a  minimum 

C  value  of  1  for  array  declaration  purposes.) 

C  MPTRS ... .Number  of  integer  data  stored  for  each  user-defined  force 
C  element.  (This  variable  has  a  minimum  value  of  1  for 

C  array  declaration  purposes.) 

C  TIME . Current  time  in  the  model  simulation. 

C  Q . Array  of  displacement  values  for  all  rigid  body  coor- 

C  dinates  in  the  model. 

C  QD . Array  of  velocities  of  all  rigid  body  coordinates. 

C  QDD . Array  of  accelerations  of  all  rigid  body  coordinates. 

C  XQ . Array  of  displacement  values  for  all  extra  generalized 

C  coordinates  in  the  model. 

C  XQD . Array  of  velocities  of  all  extra  generalized  coordinates. 

C  XQDD . Array  of  accelerations  of  all  extra  generalized  coor- 

C  dinates. 

C  NXCRD. .. .Number  of  extra  generalized  coordinates  in  the  model. 

C  (This  variable  has  a  minimum  value  of  1  for  array  dec- 

C  laration  purposes.) 

C  RB . Array  of  real  (double  precision)  data  for  each  rigid  body 

C  in  the  model. 

C  NRB . Number  of  rigid  bodies  in  the  model. 

C  NPRB . Number  of  real  (double  precision)  data  stored  for  each 

C  rigid  body  in  the  model. 

C  NTYPES ... Numbe r  of  different  element  types  available  in  DADS. 

C  NELEMS ... Numbe r  of  occurrences  of  each  element  type  in  the  model. 

C  A . Array  of  all  real  (double  precision)  data  used  in  the 

C  analysis. 

C  NMAX . Length  of  the  A  array. 

C  IA . Array  of  all  integer  data  used  in  the  analysis. 

C  MMAX . Length  of  the  IA  array. 

C  IPRDCT. .. Flag  used  to  identify  the  type  of  integration  step  cur- 
C  rently  being  performed: 

C  If  IPRDCT  -  0,  this  is  a  corrector  step. 

C  If  IPRDCT  >■  1,  this  is  a  predictor  step. 

C  INFOF ....  FORTRAN  file  unit  number  for  the  DADS  information  file. 

C  Error,  warning,  or  other  informational  messages  can  be 

C  written  to  this  file  if  so  desired. 

C 

INTEGER  NUDF,  NPTRS,  MPTRS,  NXCRD,  NRB,  NPRB,  NTYPES, 

4  NELEMS ( NTYPES ) ,  NMAX,  MMAX,  IA(0:MMAX),  IPRDCT, 

4  INFOF 

DOUBLE  PRECISION  TIME,  Q ( 7 , NRB ) ,  QD( 7 , NRB ) ,  QDD( 7 ,NRB ) , 

4  XQ( NXCRD),  XQD( NXCRD),  XQDD( NXCRD), 

4  RB ( NPRB , NRB ) ,  A(0:NMAX) 


onnnn  nnnnnnonoooo  ononn  oooonnnnnnon 


- Arguments  read  and  modified- 


UDF . Array  of  real  (double  precision)  data  for  all  user- 

defined  force  elements  in  the  model. 

IUDF . Array  of  integer  data  for  all  user-cr  . -fined  force  elements 

in  the  model. 

FRC . Array  of  forces  for  each  rigid  body  coordinate  in  the 

model . 

XFRC . Array  of  forces  for  each  extra  generalized  coordinate  in 

the  model. 

INTEGER  IUDF ( MPTRS , NUDF ) 

DOUBLE  PRECISION  UDF ( NPTRS , NUDF ) ,  FRC ( 7 , NRB ) ,  XFRC(NXCRD) 


Return-only  arguments- 


ERRCOD. .. Error  condition  flag. 
INTEGER  ERRCOD 


COMMON  blocks- 
none 


Local  variables- 


Functions  and  subroutines- 


SPL32 ....  Function  used  to  return  various  values  associated  with  a 
a  point  on  a  particular  curve  data  element. 

EXTERNAL  SPL32 

DOUBLE  PRECISION  SPL 3 2 


DATA  statements- 


C  «  »  #  «  »  *  *  «  4  *  *  t  t  *  *  I  I  *  t  *  4  «  *  t  «  *  *  t  *  *  *  »  t  * 

Remaining  Code  Supplied  By  User 

C  *»***t»##l»«###*#*##*i**t**#*#*»*« 


common  /fric/  friction,  velocity,  accel,  ext_torque 
common  /hook/  act_force,  pi,  p2 
integer  body_axis,  body  number 

double  precision  act_force,  pl(3),  p2(3) 

double  precision  friction,  dt,  fudge_f actor ,  inertia,  torque, one 
double  precision  last  vel,  last_time,  step_flag,  velocity 
double  precision  max_Friction,  stopping  torque,  accel , ext_torque 
double  precision  euler_params( 4 ) ,  euler~torques( 4 ) , 

eule r~accel ( 4 ) ,  jt_locT3),  torque_arm( 3 ) 
double  precision  euler~velocity( 4 ) 
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c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 
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double  precision  local_torques ( 3 ) ,  local_velocity( 3 ) , 

local  accel(3) 

double  precision  forc(7),  unit(3), 

ptecml(3),  pterm2(3), 

torq2 { 4 ) ,  dist,  posl(3),  pos2(3), 

value,  tterml(3),  tterm2(3) 


save  /hook/ 
save  /fric/ 


############  VARIABLES  »#*»****##*«* 


accel 
act  force 


body_axis 

body_number 

dist 

dt 

ext_torque 

friction 
fudge_f actor 


inertia 

last_time 

last_vel 

max_f riction 
moment  arm 
step_fTag 

stopping_torque 

torque 

velocity 


--  Acceleration  of  body  about  joint  axis 

--  Force  produced  by  the  actuator,  value  entered 
into  FRC32  3D  through  a  common  block,  value 
in  common  Block  set  by  modified  FR3520  routine 

—  Axis  number  about  which  the  torque  is  applied 

—  Number  of  the  body  of  interest 

—  Distance  between  the  actuator's  end  points 

—  Si2e  of  simulation's  calculation  time  step 

—  Total  external  torque  excluding  friction  that 
is  acting  on  the  body  about  the  joint  axis 

—  Value  used  for  friction 

—  Value  used  to  adjust  the  stopping  torque  when 
the  friction  torque  is  sufficient  to  stop  the 
inertia  in  a  time  step 

—  Inertia  of  the  body  in  motion  in  the 
direction  of  the  torque 

—  Value  of  time  at  last  completed  simulation 
time  step 

—  Value  of  the  body's  velocity  at  the  last 
completed  time  step 

—  Maximum  torque  due  to  friction 

—  Distance  from  body  CG  to  pivot  point 

—  Flag  used  to  determine  if  a  time  step  has 
been  completed,  («  -1  if  yes,  «  time  if  no) 

—  Approximate  torque  necessary  to  stop  body's 
motion  in  the  next  time  step 

—  Torque  applied  to  the  body  about  the  friction 
axis 

—  Velocity  of  the  body 


ARRAYS  t  ############  # 


eule  r_accel 
euler_params 
euler_torques 
fore 

j  t_loc 

local_accel 

local~torques 

local~velocity 

posl 

pos2 

pterml 

pterm2 

torque_arm 


—  Angular  acceleration  of  body  in  euler  c-sys. 

—  Euler  parameters  for  body  of  interest 

—  Global  euler/cartesian  torques  acting  on  body 

—  Force  vector  in  global  coordinates  from  end  1 
to  end  2  of  actuator 

—  Global  location  of  the  joint  axis 

—  Acceleration  of  body  2  in  its  coordinate  system 

—  Body  fixed  torques  acting  on  the  body 

—  Velocity  of  body  2  in  its  coordinate  system 

—  Global  location  of  actuator's  first  end 

—  Global  location  of  actuator's  second  end 

—  Location  of  actuator's  first  end  on  body  one 
with  respect  to  that  body's  CG 

—  Location  of  actuator's  second  end  on  body  two 
with  respect  to  that  body's  CG 

—  Vector  from  joint  location  to  actuator's 
second  end  on  body  2,  in  global  coordinates 


nnnnnnn 


torq2 


tte  rm2 
unit 


One  half  of  the  resultant  torque  acting  on 
body  2  because  of  the  actuator,  in  Euler 
parameter  coordinate  system  • 

Resultant  torque  acting  on  body  2  because  of 
the  actuator,  in  global  coordinate  system 
Unit  vector  along  actuator  from  end  one  to 
end  two,  in  global  coordinates 


C 


Process  Block 


C - First  zero  out  the  error  condition  flag  to  indicate  that  no  errors 

C  have  occurred  yet. 

ERRCOD  -  0 


i f  ( time  . eq .  0 ) 

max_f riction 
fudge_f actor 
body_numbe  r 
body_axis 
moment_arm 
j  t_loc ( 1 ) 
jt_loc( 2 ) 
j  t_loc ( 3 ) 


then 

-  udf(l,l) 

-  udf ( 2 , 1 ) 

-  int( udf ( 3 , 1 )  ) 

-  int ( udf (4,1)) 

-  udf ( 5 , 1 ) 

-  udf ( 7 , 1 ) 

-  udf  (  8 , 1 ) 

-  udf ( 9 , 1 ) 


one  time  calcs 


C  ***  Calculation  of  inertia  about  joint  axis 


1 


inertia 


RB{ body_axis+l ,  body_number)  + 

RB(1,  body_number ) *moment_arm**2 


step  flag  -  0 
dt  “  -  1 
one  «  1 


endi  f 


C  ***  If/then/else  statement  used  to  capture  simulation  time  and  body's 
C  ***  velocity  about  joint  axis  at  last  COMPLETED  time  step 

if  (time  .eq.  step_flag)  then 
step_flag  -  -1 
last_vel  -  velocity 
last  time  ■  time 


else 

step_flag  -  time 
dt  -  time  -  last_time 
if  (time  .eq.  0)  dt  -  1 

endi  f 


do  i  -  1,  4 
j  -  i  +  3 

euler_params( i )  -  Q(j,  body  number) 

euler_torques( i )  »  FRC ( j ,  bo3y_number) 
euler_accel ( i )  ■  QDD(j,  body_number) 


enddo 


C - Transform  pi  and  p2  from  local  to  global  coordinates. 

call  asp{  RB ( 38 , 1 )  ,  pl(l),  pterml(l)  ) 

call  asp(  RB{ 38 , body_number ) ,  p2(l),  pterm2(l)  ) 

posl(l)  =  Q( 1 , 1 )  +  pterml(l) 
posl(2)  -  Q( 2 , 1 )  +  pterml(2) 
posl(3)  -  Q(  3 , 1 )  +  pterml(3) 

pos2(l)  -  Q  ( 1 , body_numbe r )  +  pterm2(l) 
pos2(2)  -  Q{ 2 , body_numbe r )  +  pterm2(2) 
pos2(3)  -  Q( 3 , body_numbe r )  +•  pterm2(3) 


C  ***  Calculate  actuator  length 

dist  -  (  pos2(l)  -  posl(l)  )**2  +  (  pos2(2)  -  posl(2)  )**2 
&  +  (  pos2{3)  -  posl(3)  )**2 

if  (  dist  .gt.  l.d-10  )  then 
dist  -  sqrt(dist) 
else 

write( INFOF, 1000 ) 

1000  formate/,'  Error  from  the  control  output  element.', 

&  /, '  The  distance  is  near  zero,  a  value  of', 

&  '  l.D-5  will  be  used.  '  ) 

dist  -  l.d-5 
endif 


C - Define  the  unit  vector  from  PI  to  P2  in  global  space. 

unit(l)  =  (  pos2{l)  -  posl(l)  )  /  dist 
unit(2)  -  (  pos2(2)  -  posl(2)  )  /  dist 
unit(3)  -  (  pos2(3)  -  posl(3)  )  /  dist 


C  ***  Calculate  the  force  vector  acting  in  body  2  due  to  the  actuator 

forc(l)  -  act_force  *  unit(l) 
forc(2)  =  act_force  *  unit(2) 
forc(3)  ■  act  force  *  unit(3) 


q  ***  Calculate  the  torque  acting  on  body  2  due  to  the  actuator  and 
C  ***  convert  to  the  euler  parameter  c-sys 

do  i  ■  1,  3 

torque_arm( i )  -  pos2(i)  -  jt_loc(i) 
enddo 

call  WCS(  torque  arm(l),  forc(l),  tterm2(l)) 
call  ETS (  Q( 4 , bo3y_number ) ,  tterm2(l),  torq2(l)  ) 


C  ***  Add  the  estimate  of  the  actuator  torque  to  the  other  torques 
C  ***  already  acting  on  the  body 

do  i  -  1,  4 


euler_torques  (  i  )  >■  euler_torques( i )  +  2.d0  *  torq2(i) 
enddo 


C  ***  Calculate  the  torques  and  accelerations  of  the  body  in  its  c-sys 


call  GV4 ( eule r_pa r ams ,  euler_accel,  local_accel) 
call  GV4 ( euler_params ,  euler_torques,  local_torgues ) 


do  i  *  1 ,  3 

local^accel ( i )  =  local_accel ( i )/2 
local'torques ( i )  -  local_torques ( i )/2 
local_velocity( i )  -  RB(34+i,  body_number) 

enddo 


C  ***  Set  the  variables  to  be  used  in  the  friction  calculation 


velocity  *  local_velocity(body_axis ) 
accel  «  local  accel { body_axis )~ 
torque  «  locaT_torques ( body_axis ) 
ext_torque  -  torque 


Calculation  of  Friction 


if  (abs(torque)  .ge.  max  friction)  then 


if  (velocity  .eq.  0)  then 


C  =*=*  Set  friction  to  its  max  value  and  in  the  opposite  direction  of  the 
C  net  external  torque  if  the  body  is  stopped  otherwise  set  it  in 

C  opposition  to  the  body's  motion 


friction  -  -signCone,  torque)*max  friction 


else 


friction  -  -sign(one,  veloci ty ) *max  friction 


endi  f 


if  (velocity  .eq.  0)  then 
friction  =  -torque 


Cancel  out  the  external  torq 


C  ***  Estimate  the  torque  necessary  to  stop  the  body  in  one  time  step 


stopping_torque  «  -torque  -  velocity  * 

inertia*f udge_f actor/dt 


C  ***  Set  the  value  of  friction  to  that  torque  value  enless  it  exceeds 
C  ***  friction's  maximum  possible  value 


if  (abs( stopping_torque)  .ge.  max_f riction)  then 
friction  ■  -sign(one,  velocity ) *max_f riction 


else 


friction  ■  stopping  torque 


to 


endi  f 
endi  f 
endi  f 


C  ***  Set  the  local  torque  vector  to  be  equal  to  the  friction  present 

do  i  -  1,  3 

local_torques ( i )  -  0 

enddo 

local_torques ( body_axi s )  *  friction 

C  ***  Convert  the  frictional  torque  to  the  euler  parameter  c-sys  and 
C  ***  add  it  to  the  external  torques  already  acting  on  the  body 

call  GTS  ( euler_params ,  local_torques ,  euler_torques ) 

do  i  -  1,  4 
j  -  i  +  3 

frc(j,  body_number)  ■  2*euler_torques ( i )  + 

*  frc(j,  body  number) 

frc(j,  1)  *  -2*euler_torques ( i )  +  frc(j  ,  1) 
enddo 


re  tu  rn 
end 
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non  n n n n n n n n n n n n n n n n n  n n n n n n n n n n n n n n n n n n 


This  routine  is  used  to  apply  forces  calculated  by  DADS  control 
module  elements  to  bodies  in  the  model.  Typically  it  is  not 
modified  by  the  user.  To  model  friction  it  was  modified  only  to 
gain  access  to  some  model  variables  that  where  needed  and  that 
where  not  available  to  the  user  written  force  calculation 
routine.  The  variables  were  made  availible  through  the  addition 
of  a  common  block  to  the  two  routines. 


C  FR3520:  Transfer  forces  due  the  control  elements  to  DADS. 

--Author  and  Modification  Section------------------------ 

Author:  Dick  Kading 

Date  written:  April  1986 
Written  on:  MicroVax,  VMS 

Modi f i cat ions : 

1)  6/2/86  Changed  the  TORQ1.2  vectors  dimension  from  3  to  4. 

Dick  Kading 

2)  7/7/86  Changed  this  routine  so  it  could  be  used  or  modified 
by  customers  wanting  special  access  to  control  forces. 
Dick  Kading 

Copyright  (c)  CADSI  1986 
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:*««, *««»»=» Algor i thm  Description— 

Purpose  and  use: 

This  routine  appends  the  forces  calculated  by  the  control  elements 
to  the  generalized  force  vector.  There  are  8  possible  options  in 
this  element  that  apply  forces,  some  to  one  body,  some  to  both. 

This  routine  has  the  common  blocks  locally  declared,  not  "INCLUDED", 
so  that  it  can  be  edited,  recompiled  and  linked  as  a  "user"  routine. 
It  provides  access  to  all  rigid  body  position,  velocities, 
accelerations,  and  forces  along  with  the  control  node  and  state 
variables . 

Error  conditions:  none 

Machine  dependencies:  none 


SUBROUTINE  FR3520  ( ENS , IC , RC , Q , QD , QDD , FRC , NEL , RB , RVLT , CYL , S , SD ) 

---••------•Variable  Descriptions------ ------ ----------------------- 

- Arguments  passed - 


non 


C 

C  NEL . Number  of  outputs  to  the  rigid  body  program. 

C  IC . Integer  array  of  output  element  data. 

C  ENS . Vector  of  node  variable  values. 

C  RC . Real  array  of  control  element  data. 

C  Q . Real  array  of  rigid  body  generalized  coordinates. 

C  QD . Real  array  of  rigid  body  velocities. 

C  QDD . ...Real  array  of  rigid  body  accelerations. 

C  FRC . Real  array  of  generalized  forces. 

C  RB . Array  of  rigid  body  real  data. 

C  RVLT . Array  of  revolute  joint  real  data. 

C  CYL . Array  of  cylindrical  joint  real  data. 

C  S . Vector  of  state  variables. 

C  SD . Vector  of  state  variable  time  derivatives. 

C 

INTEGER  NEL, IC(NEL, 1 ) 

DOUBLE  PRECISION  ENS(l),  RC(NEL,1),  Q(7,l),  FRC(7,1),  RB(65,1), 

&  RVLT( 36,1),  CYL ( 36,1) ,  QD(7,1),  QDD(7,1),  S(l), 

&  SD( 1 ) 

C 


C - COMMON  blocks - 

C 

C  INPUTF . . .File  unit  the  input  data  is  in. 

C  INFOF . . . . Fi le  unit  cf  the  information  output  file. 

C  OUTPF . . . . Fi le  unit  of  the  standard  output  file. 

C  BINRYF . . . Fi le  unit  of  the  unformatted  or  "binary"  output  file. 
C  CRTIN . . . . File  unit  for  reading  from  the  user's  terminal. 

C  CRTOUT . . . Fi le  unit  for  writing  to  the  user's  terminal. 

C  NXTREC . . .Next  record  available  in  the  binary  output  file. 


INTEGER  INPUTF,  INFOF,  OUTPF,  BINRYF,  CRTIN,  CRTOUT,  NXTREC 

COMMON  /  IOUNIT  /  INPUTF,  INFOF,  OUTPF,  BINRYF,  CRTIN,  CRTOUT, 
&  NXTREC 

SAVE  /IOUNIT/ 

- Controls  COMMON  Block. 

COMMON  /STAT  3  5/  I EVAL , I MODUL , I  BLOCK , I FN 
INTEGER  IEVAL, IMODUL, IBLOCK, IFN 
SAVE  /STAT3  5/ 

C 

C - Local  variables - 


C 

C  I . Loop  index  . 

C  I  OUT . Output  node  number. 

C  TYPE . Type  of  coordinate  force  is  appended  to, 

C 

C  TYPE  -  1,  X 

C  -  2,  Y 

C  -  3,  Z 

C  «  4,  Force  difference 

C  -5,  Torque  difference 

C  -  6,  X  torque,  local 

C  -  7,  Y  torque 

C  -  8,  Z  torque 

C 

C  BOD  Yl . Body  1. 

C  BODY2 . Body  2. 


no  non 


C  FORCE . Value  of  the  force  output  from  the  control  system. 

C  FORC . Vector  of  forces  in  global  space. 

C  UNIT . Unit  vector. 

C  SLP1 . Vector  locating  Pi  on  body  1,  in  local  coordinates. 

C  SLP2 .  "  P2  on  body  2. 

C  PTERMl . Global  value  of  SLP1. 

C  PTERM2 .  "  SLP2. 

C  T0RQ1 . Torque  on  body  1,  transformed  into  the  Euler  parameter 

C  coordinate  system. 

C  T0RQ2 .  "  body  2. 

C  DIST . Distance  from  Pi  on  bodyl  to  P2  on  body  2. 

C  POS 1 . Global  position  of  PI  on  body  1. 

C  P0S2 .  "  P2  on  body  2. 

C  VALUE . Temporary  value  used  to  test  the  SQRT  argument. 

C  TTERMl . Torque  term,  result  of  cross  product  of  Force  and 

C  moment  arm. 

C  TTERM2 .  "  body  2. 

C  AXIS1L . Revolute  or  cylindrical  axis  about  which  the  actuator 

C  torque  is  applied  on  body  1. 

C  AXIS2L .  "  body  2. 

C  JNTNUM . Joint  number,  either  a  revolute  or  cylindrical. 

C  JNTTYP . Type  of  joint,  revolute  -  1,  cylindrical  »  3. 

INTEGER  I,  IOUT ,  B0DY1 ,  BODY2 ,  TYPE,  JNTNUM,  JNTTYP 

DOUBLE  PRECISION  FORCE,  FORC(3),  UNIT(3),  AXIS1L(3),  AXIS2L(3), 
&  SLPl ( 3  )  ,  SLP2 ( 3 ) ,  PTERMl ( 3 ) ,  PTERM2(3), 

&  T0RQ1 ( 4 ) ,  TORQ2 ( 4 ) ,  DIST,  POSl(3),  POS2(3), 

&  VALUE,  TTERMl ( 3 ) ,  TTERM2 ( 3 ) 


- Functions  and  subroutines - 

INTRINSIC  ABS,  SQRT 
EXTERNAL  ASP,  WCS ,  ETS ,  GTS 


Cl  «  |  «  I  |  I  |  |  |  «  4  «  t  *  t  «  t  I  »  I  I  «  «  I  I  I  I  »  I  »  *  « 

USER  MODIFICATIONS 
common  /hook/  act_force,  pi,  p2 
double  precision  act_force,  p 1 ( 3 ) ,  p 2 ( 3 ) 
save  /hook/ 

C*  H  t  «  I  I  M  I  U  t  t  I  «  t  M  U  I  M  U  I  t  I  . . 

C  «  U  «  i  U  n  M  U  I  I  M  M  M  M  II  U  U  M  «  <  t  t 

C-----«- -----Process  Block ■  —  - 

C - If  this  isn't  the  right  analysis  phase  return. 

IF  (  IEVAL  .NE.  4  )  RETURN 

C###*####*#**#i***#«#*#*tf 


4>5 

USER  MODIFICATIONS 

act  force  -  ens(13) 
do  I  -  1 ,  3 

pi ( i )  -  rc(  1 ,  i ) 
p2( i )  -  rc( 1 , i  +  3  ) 

enddo 


C - Go  through  the  loop  for  each  output  required. 

DO  100  I  -  l.NEL 

C - Define  local  variables  for  convenience. 


IOUT  -  I C ( 1,1) 

TYPE  -  ICC  1,2) 

BODYl  -  ICCI.3) 

BODY2  -  I C ( I , 4  ) 

JNTNUM  -  ICC  I , 7 ) 

JNTTYP  -  ICC  1,8) 

C - Define  the  vectors  locating  points  Pi  and  P2 . 

.SLP1  (  1  )  -  RC(  I  ,  1  ) 

SLPl ( 2)  -  RC (1,2) 

SLP1 (  3  )  -  RC( I , 3 ) 

SLP2 ( 1  )  -  RC (1,4) 

SLP2C2)  -  FCCI.5) 

SLP2C  3 )  -  RC( I ,6 ) 

C - Define  the  output  value  of  the  node. 

FORCE  -  ENSC IOUT) 

C - Transform  SLPl  and  2  from  local  to  global  coordinates. 

CALL  ASP  C  RBC  38 , BODYl ) ,  SLP1C1),  PTERMl ( 1 )  ) 

CALL  ASPC  RB ( 38 , BODY2 ) ,  S L P 2 ( 1  )  ,  PTERM2 ( 1 )  ) 


i 

s 

« 


FORC(l)  -  FORCE 
FORC ( 2 )  -  0.0 
FORC ( 3)  -  0.0 

C - Call  WCS  to  calc,  the  cross  product  of  thte  moment  arm  and  the 

C  force.  Use  ETS  to  transform  the  the  torque  to  the  Euler  parameter 

C  coordinate  system. 

CALL  WCS(  PTERM1 (  1  )  ,  FORC ( 1  )  ,  TTERMl ( 1 )  ) 

CALL  ETS (  Q{ 4, BODYl),  TTERMl ( 1 ) ,  TORQ1 ( 1 )  ) 

C - Update  the  generalized  force  vector. 


s 

s 


FRC ( 1 , BODY1 ) 
FRC ( 4 , BODYl ) 
FRC ( 5, BODYl ) 
FRC ( 6 , BODYl ) 
FRC ( 7 , BODYl ) 


FRC ( 1 , BODYl  ) 
FRC ( 4 , BODYl  ) 
FRC ( 5, BODYl  ) 
FRC ( 6 , BODYl  ) 
FRC( 7 , BODYl  ) 


+  FORCE 
♦  2. DO  * 
+  2. DO  * 
t-  2.  DO  * 
f  2 . DO  * 


TORQ1 ( 1  ) 
TORQ1 ( 2 ) 
TORQ1 ( 3  ) 
TORQ1 ( 4  ) 


C - This  is  the  Y  of  point  PI  option. 


ELSE  IF  (  TYPE  .EQ.  2  )  THEN 

FORC ( 1 )  -  0.0 
FORC ( 2 )  -  FORCE 
FORC ( 3 )  -  0.0 


CALL  WCS(  PTERMl(i),  FORC(l),  TTERMl) 1)  ) 

CALL  ETS (  Q( 4 .BODYl),  TTERMl ( 1 ) ,  TORQl(l)  ) 


FRC ( 2 , BODYl  ) 

M 

FRC ( 2 , 

BODYl  ) 

+ 

FORCE 

FRC ( 4 , BODYl ) 

• 

FRC  S 4 , 

BODYl  ) 

+ 

2. DO  * 

TORQl ( 1 ) 

FRC ( 5 , BODYl  ) 

m 

FRC ( 5, 

BODYl  ) 

4- 

2  .DO  * 

TORQl ( 2 ) 

FRC ( 6 , BODYl ) 

m 

FRC ( 6 , 

BODYl  ) 

+ 

2. DO  * 

TORQl ( 3 ) 

FRC (  7 .BODYl  ) 

m 

FRC ( 7 , 

BODYl  ) 

4 

2  .DO  * 

TORQl ( 4 ) 

C - This  is  the  Z  of  point  PI  option. 


ELSE  IF  (  TYPE  .EQ.  3  )  THEN 


FORC ( 1 )  -  0.0 
FORC ( 2 )  -  0.0 
FORCt  3 )  -  FORCE 


ft 

ft 


& 


A' 


J  ‘ 


CALL  WCS  (  PTERMl  (  1  )  ,  FORC'D,  TTERMl  (  1  )  ) 

CALL  ETS (  Q( 4, BODYl),  TTERMl ( 1 ) ,  TORQ1 ( 1  )  ) 

FRC(3, BODYl)  -  FRC(3, BODYl)  +  FORCE 

FRC(4, BODYl)  -  FRC(4, BODYl)  +  2. DO  *  TORQl ( 1 ) 

FRC(5, BODYl)  -  FRC(5, BODYl)  +  2. DO  *  TORQl ( 2 ) 

FRC(6, BODYl)  -  FRC(6, BODYl)  +  2. DO  *  TORQl ( 3 ) 

FRC(7, BODYl)  -  FRC(7, BODYl)  +  2. DO  *  TORQl ( 4 ) 

C - Update  generalized  force  with  force  difference  between  the  2 

C  bodies  connected.  This  is  to  apply  force  at  PI  on  body  1  and 

C  P2  on  body  2. 

ELSE  IF  (  TYPE  .EQ.  4  )  THEN 

VALUE  -  (  POS2 ( 1 )  -  POSl(l)  ) *  *  2  (  POS2(2)  -  POS 1(2)  )**2 

&  +  (  POS  2 ( 3 )  -  POS1 ( 3  )  )  *  *  2 


& 


Eft! 


ru  v  ..  .  <c+  c 


.  <. 


IF  (  VALUE  .GT.  l.D-10  )  THEN 
DIST  -  SQRT! VALUE ) 

ELSE 

WRITE( INFOF, 1000 ) 

1000  FORMAT ( / , '  Error  from  the  control  output  element.', 

6  /,'  The  distance  is  near  zero,  e  value  of', 

a  '  l.D-5  will  be  used.  '  ) 

DIST  -  l.D-5 
END  IF 


C - Define  the  unit 

UNIT( 1 ) 
UNIT( 2 ) 
UNIT( 3 ) 

FORC ( 1  ) 
FORC ( 2 ) 
FORC ( 3  } 


vector  from  PI  to 

-  (  POS2 ( 1 )  -  POS1 

-  (  POS2 ( 2 )  -  POS1 

-  (  POS2 ( 3)  -  POS1 

-  FORCE  *  UNIT(l) 

-  FORCE  *  UNIT( 2 ) 

-  FORCE  *  UNIT( 3 ) 


P2  in  global  space. 

(1)  )  /  DIST 

(2)  )  /  DIST 

(  3 )  )  /  DIST 


C - Calculate  the  torques  and  convert  to  the  Euler  parameter  system. 


CALL  WCS(  PTERMl(l),  FORC(l),  TTERMl ( 1 )  ) 

CALL  ETS (  Q ( 4 , BODYl ) ,  TTERMl ( 1 ) ,  TORQ1 ( 1 )  ) 

CALL  WCS(  PTERM2 ( 1 ) ,  FORC(l),  TTERM2 ( 1 )  ) 

CALL  ETS (  Q  (  4  , BODY2 ) ,  TTERM2 ( 1 ) ,  TORQ2(l)  ) 


the  generalized 


C - Update 


FRC ( 1 , BODYl ) 
FRC(2,B0DY1) 
FRC ( 3, BODYl ) 
FRC ( 4 , BODYl ) 
FRC (5, BODYl) 
FRC ( 6 , BODYl ) 
FRC( 7 .BODYl ) 


force  vector. 

FRC ( 1 , BODYl ) 
FRC( 2 .BODYl ) 
FRC ( 3 , BODYl ) 
FRC{ 4 .BODYl ) 
FRC ( 5 , BODYl ) 
FRC( 6 .BODYl ) 
FRC! 7 .BODYl ) 


FORC ( 1  ) 

FORC! 2  ) 

FORC! 3  ) 

2. DO  *  TORQl(l) 
2. DO  *  TORQ1 ( 2 ) 
2. DO  *  TORQ1 ( 3 ) 
2. DO  *  TORQ1 ( 4 ) 


FRC! 1 , BODY2 ) 
FRC! 2 .BODY2 ) 
FRC! 3 , BODY2 ) 
FRC! 4 .BODY2 ) 
FRC! 5.BODY2 ) 
FRC ( 6 , BODY2 ) 
FRC ( 7 , BODY2 ) 


FRC ( 1 , BODY2 )  ♦ 
FRC ( 2 , BODY2  )  + 
FRC ( 3 , BODY2 )  ♦ 
FRC ( 4 , BODY2  )  + 
FRC ( 5 , BODY2  )  + 
FRC ( 6 , BODY2  )  + 
FRC ( 7 , BODY2  )  + 


FORC! 1 ) 

FORC! 2  ) 

FORC! 3 ) 

2. DO  *  TORQ2 ( 1 ) 
2. DO  *  TORQ2 ( 2 ) 
2. DO  *  TORQ2 ( 3 ) 
2. DO  *  TORQ2 ( 4 ) 


C - Update  generalized  force  with  rotational  actuator  force.  This 

C  option  applies  torque  on  the  2  bodies  connected  by  revolute  or 

C  cylindrical  joints.  The  torque  is  applied  about  the  rev  or  cyl 
C  joint  Z  axes. 

ELSE  IF  (  TYPE  .EQ.  5  )  THEN 


C - If  this  is  a  revolute  joint,  use  the  RVLT  data  to  define  the 

C  local  joint  Z  axis.  The  force  term  is  included  at  this  step; 
C  it  could  be  applied  at  the  last  step  equivalently. 


IF  (  JNTTYP  .EQ.  1  )  THEN 

AXIS1L(1)  -  RVLT f  2  5  , JNTNUM )  *  FORCE 
AXISlL! 2  )  -  RVLT! 26, JNTNUM)  *  FORCE 
AXIS1L! 3 )  -  RVLT! 27, JNTNUM)  *  FORCE 


iX  V 


vacw- 
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AXIS2L(1)  -  RVLT{ 34, JNTNUM)  *  FORCE 
AXIS2L(2)  -  RVLT< 35, JNTNUM)  *  FORCE 
AXIS2L(3)  -  RVLT( 36 , JNTNUM)  *  FORCE 

If  this  is  a  cylindrical  joint,  use  the  CYL  data. 

ELSE  IF  (  JNTTYP  . EQ .  3  )  THEN 

AXIS  1 L ( 1 )  -  CYL( 2  5, JNTNUM)  *  FORCE 
AXIS1L(2)  -  CYL{ 26, JNTNUM)  *  FORCE 
AXIS1L{ 3 )  -  CYL( 27, JNTNUM)  *  FORCE 

AX I S 2 L ( 1 )  -  CYL ( 34, JNTNUM)  *  FORCE 
AXIS2L(2)  -  CYL{ 35 , JNTNUM)  *  FORCE 
AXIS2L ( 3 )  -  CYL( 36, JNTNUM)  *  FORCE 

END  IF 

Transform  the  torque  to  the  Euler  parameter  system. 

CALL  GTS (  Q( 4 , BODYl ) ,  AXISlL(l),  TORQ1 ( 1 )  ) 

CALL  GTS (  Q  (  4  , BODY2 ) ,  AXIS2L(1),  TORQ2(l)  ) 

FRC(4, BODYl)  -  FRC(4, BODYl)  -  2. DO  *  TORQ1 ( 1 ) 

FRC(5, BODYl)  -  FRC(5, BODYl)  -  2. DO  *  TORQ1 ( 2 ) 

FRC(6, BODYl)  -  FRC{6, BODYl)  -  2. DO  *  TORQ1 ( 3 ) 

FRC ( 7 , BODYl )  -  FRC(7,BODYl)  -  2 . DO  *  TORQ1 ( 4 ) 

FRC ( 4 , BODY2 )  -  FRC(4,BODY2)  ♦  2. DO  *  TORQ2(l) 

FRC ( 5 , BODY2 )  -  FRC ( 5 , BODY2 )  +  2 . DO  *  TORQ2(2) 

FRC ( 6 , BODY2 )  -  FRC(6,BODY2)  +  2. DO  *  TORQ2(3) 

FRC  (  7  ,  BODY2  )  -  FRC(7,BODY2)  ♦  2.  DO  *  TORQ2  (  4  ) 

This  option  applies  a  torque  about  the  local  X  axis  directly 

ELSE  IF  (  TYPE  . EQ .  6  )  THEN 

TTERMl(l)  -  FORCE 
TTERMl ( 2 )  -  0.0 
TTERMl ( 3 )  -  0.0 

CALL  GTS  (  Q{ 4, BODYl),  TTERMl ( 1  )  ,  TORQ1 ( 1 )  ) 

FRC ( 4 , BODYl )  -  FRC{4,BODYl)  +  2. DO  *  TORQ1 ( 1 ) 

FRC(5, BODYl)  -  FRC(5, BODYl)  ♦  2. DO  *  TORQ1 ( 2 ) 

FRC(6, BODYl)  -  FRC(6, BODYl)  +  2. DO  *  TORQ1 ( 3 ) 

FRC(7, BODYl)  -  FRC(7, BODYl)  +  2. DO  *  TORQ1 ( 4 ) 

This  option  applies  a  torque  about  the  local  Y  axis  directly 

ELSE  IF  (  TYPE  . EQ .  7  )  THEN 

TTERMl ( 1  )  -  0.0 
TTERMl ( 2  )  -  FORCE 
TTERMl ( 3  )  -  0.0 

CALL  GTS (  Q( 4, BODYl),  TTERMl ( 1 ) ,  TORQ1 ( 1 )  ) 

FRC(4, BODYl)  -  FRC(4, BODYl)  ♦  2. DO  *  TORQ1 ( 1 ) 

FRC<5, BODYl)  -  FRC!5, BODYl)  ♦  2. DO  *  TORQ1 ( 2 ) 

FRC (6, BODYl)  -  FRC (6, BODYl)  +  2. DO  *  TORQ1 ( 3 ) 


FRC  (  7  ,  BODYl )  -  FRC (7, BODYl)  +  2.  DO  *  TORQ1  (  4  ) 


— This  option  applies  the  torque  about  the  local  Z  axis  directly. 
ELSE  IF  (  TYPE  .EQ.  8  )  THEN 


TTERMl ( 1 )  -  0.0 
TTERMl ( 2  )  -  0.0 
TTERMl ( 3 )  -  FORCE 


CALL  GTS (  Q( 4, BODYl),  TTERMl ( 1 ) ,  TORQl(l)  ) 


FRC ( 4 , BODYl ) 
FRC { 5 , BODYl ) 
FRC ( 6 , BODYl ) 
FRC { 7 , BODYl ) 


-  FRC(4, BODYl)  + 

-  FRC(5, BODYl)  + 

-  FRC<6, BODYl)  + 

-  FRC(7, BODYl)  + 


2. DO  *  TORQl(l) 
2. DO  *  TORQ1 ( 2 ) 
2. DO  *  TORQl ( 3 ) 
2. DO  *  TORQl ( 4 ) 


END  IF 


100  CONTINUE 


RETURN 

END 


nnononnnnnnonoonnonnnnononnnon  nnnnnnnnono 


This  routine  allows  the  user  to  input  parameters  that  are 
to  be  used  with  their  force  calculations.  Most  of  this  routine 
is  a  vendor  supplied  template,  the  modifications  done  by  the 
appear  near  its  end  and  are  designated  as  such. 


C  IN32:  Reads  and  stores  input  data  for  user-defined  forces  (3D). 
-----■---^Author  and  Modification  Section----— 

Author:  Chuck  Mead 

Date  written:  December  5,  1985 
Modifications : 


Copyright  (c)  CADSI  1985 


--------.-Algorithm  Descr iption------------------------------------- 

Purpose  and  use: 

This  subroutine  reads  the  user-defined  force  input  data  from  the 
internal  file  HLDINP ,  checks  it  for  errors,  and  stores  it  in  the 
user-defined  force  data  areas:  UDF  for  real  (double  precision) 
data  and  IUDF  for  integer  data.  Note  that  the  first  record  of 
each  user-defined  force  element  must  have  the  string  "USERFORCE  " 
in  columns  1  through  10  in  order  to  be  identified  as  such. 

The  number  of  records  of  data  required  for  each  user-defined 
force  element  is  defined  by  the  array  element  REC3D(32).  The 
number  of  real  and  integer  data  stored  for  each  user-defined  force 
element  are  defined  by  the  array  elements  NPTRS3(32)  and 
MPTRS3(32),  respectively.  These  three  values  are  initialized  in 
the  block  data  routine  BLOCKD.  For  use  within  this  routine,  the 
array  elements  NPTRS3(32)  and  MPTRS3(32)  are  passed  in  via  the 
argument  list  as  the  variables  NPTRS  and  MPTRS ,  respectively. 

The  various  input  arguments  are  described  and  declared  in  three 
separate  groups:  read-only  arguments,  arguments  read  and  modi¬ 
fied,  and  return-only  arguments. 

Read-only  arguments  are  passed  in  with  meaningful  values  which 
should  not  be  changed  within  this  subroutine..  This  constitutes 
the  largest  group  of  arguments  passed.  Consult  the  documentation 
on  user-defined  force  elements  for  a  description  of  the  structure 
and  contents  of  the  various  variables  and  arrays  passed  in. 

Arguments  read  and  modified  are  also  passed  in  with  meaningful 
values  which  may  be  changed  within  this  subroutine.  There  are  no 
read-and-modi fy  arguments  passed  to  IN32. 

The  final  group  of  arguments,  return-only  arguments,  a  e  passed 
in  with  unknown  values  and  should  therefore  be  assigned  some  value 


within  this  subroutine  before  they  are  used  in  any  other  computa¬ 
tion.  The  variables  in  this  group  are  the  two  user-defined  force 
element  data  arrays,  UDF  and  IUDF,  and  the  error  condition  flag 
ERRCOD.  Since  the  value  of  ERRCOD  is  used  in  the  calling  routine 
to  determine  the  success  of  the  computations  performed  in  IN32,  it 
is  suggested  that  ERRCOD  be  initialized  to  zero  before  beginning 
any  other  computations. 


Error  conditions: 

Error  conditions  are  reported  via  the  variable  ERRCOD.  A  zero 
value  indicates  no  error;  a  positive  value  indicates  some  fatal 
error,  i.e.,  files  will  be  closed  and  execution  aborted  upon 
return  to  the  calling  routine;  negative  values  are  reserved  for 
future  use. 


Machine  dependencies:  none  known 


SUBROUTINE  IN32  (  UDF,  IUDF,  NUDF ,  NPTRS ,  MPTRS ,  HLDINP,  MAXLNS , 

&  HLDLNK ,  STARTU,  INFOF ,  BODMOD,  CRVMOD ,  ERRCOD  ) 


-Variable  Descriptions- 


: - Read-only  arguments- 


NUDF. 


NPTRS . . 


MPTRS . . 


HLDINP. 


MAXLNS . 
HLDLNK. 


STARTU. 


INFOF .  . 


BODMOD , 


CRVMOD. 


..Number  of  user-defined  force  elements  in  the  input  data. 
(This  variable  has  a  minimum  value  of  1  for  array  dec¬ 
laration  purposes.) 

..Number  of  real  (double  precision)  data  stored  for  each 
user-defined  force  element.  (This  variable  has  a  minimum 
value  of  1  for  array  declaration  purposes.) 

..Number  of  integer  data  stored  for  each  user-defined  force 
element.  (This  variable  has  a  minimum  value  of  1  for 
array  declaration  purposes.) 

..Array  of  character  data  read  into  memory  from  the  DADS 
input  file.  This  is  used  as  a  FORTRAN  internal  file  in 
this  subroutine  in  order  to  read  any  user-defined  force 
element  data  and  store  it  in  the  appropriate  data  arrays. 

..The  maximum  number  of  lines  storable  in  HLDINP. 

..Array  of  pointers  used  to  link  together  the  locations  of 
various  occurrences  of  each  element  type  in  the  input 
data . 

..Pointer  to  the  line  in  HLDINP  at  which  the  data  for  the 
first  user-defined  force  element  starts. 

..FORTRAN  file  unit  number  for  the  DADS  information  file. 
Error,  warning,  or  other  informational  messages  can  be 
written  to  this  file  if  so  desired. 

..Element  type  number  of  the  rigid  body  module.  Can  be 
used  to  map  body  names  to  body  numbers  via  the  GETNUM 
function . 

..Element  type  number  of  the  curve  data  module.  Can  be 
used  to  map  curve  data  element  names  to  curve  data 
element  numbers  via  the  GETNUM  function. 


INTEGER  NUDF,  NPTRS,  MPTRS,  MAXLNS,  HLDLNK ( MAXLNS ) ,  STARTU, 
4  INFOF,  BODMOD,  CRVMOD 


CHARACTER  *  8  0  HLDINP ( MAXLNS ) 
CHARACTER* 2 0  UDFNAM 


ms 


C - Arguments  read  and  modified - 

C  none 

C 

C - Return-only  arguments - 

C 

C  UDF . Array  of  real  (double  precision)  data  for  all  user- 

C  defined  force  elements  in  the  model. 

C  IUDF . Array  of  integer  data  for  all  user-defined  force  elements 

C  in  the  model. 

C  ERRCOD . . .Error  condition  flag. 

C 

INTEGER  IUDF ( MPTRS , NUDF ) ,  ERRCOD 


DOUBLE  PRECISION  UDF ( NPTRS , NUDF ) 


C 

C - COMMON  blocks - 

C  none 
C 

C - Local  variables - 

C 

C  PNTR . Pointer  used  to  step  through  the  internal  file  HLDINP  and 

C  point  to  the  location  of  the  first  record  for  each  user- 

C  defined  force  element  in  the  input  data. 

C  IOSTAT ...  FORTRAN  I/O  error  flag.  A  zero  value  indicates  no  error 
C  encountered  in  the  execution  of  a  I/O  statement. 

C  I . Loop  counter  . 

C 


INTEGER  PNTR,  IOSTAT,  I 
CHARACTER* 60  MNAME , MNAMEX 
C 

C - Functions  and  subroutines - 

C 

C  GETNUM. . .Given  the  element  type  number  (i.e.,  module  number)  and 
C  an  element  name,  this  function  returns  the  element  num- 

C  ber.  A  zero  function  value  means  either  the  named  ele- 

C  ment  doesn't  exist  or  the  input  values  are  invalid. 

C 

EXTERNAL  GETNUM 

INTEGER  GETNUM 
C 

C - DATA  statements - 

C  none 
C 


C=======»====Process  Block  =  =  =  =  «»  =  *=  =  =  =  =  =  »*  =  =  ”  =  »»=»»»»«»»-»«»»»»«»»» 

C - No  errors  have  been  detected  yet. 

ERRCOD  =  0 

C - If  one  or  more  extra  output  files  will  be  used  to  report  user- 

C  defined  force  output  data,  they  could  be  opened  at  this  point. 

C - Set  the  pointer  to  the  first  line  of  data  for  the  first  user- 

C  defined  force  element. 

PNTR  -  STARTU 

C - Loop  to  read  in  and  process  the  data  for  each  user-defined  force 


•  *•>  r. 
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C  element  in  the  input  data. 


74- 


A. 


USER  MODIFICATIONS 


PI 


(START  AT  THIS  POINT  AND  APPEAR  AS  lower  case  CODE.) 


C  ***  OPEN  A  FILE  FOR  ADDITIONAL  OUTPUT  SPECIFIED  BY  THE  USER  IN 
C  ***  ROUTINE  RPT32_3D 

write  (*,'(lx,  a,  $)')  'EXTRA  DATA  FILE  NAME?  ' 
read  ( *  ,  ' ( a ) ' )  mname 

nn  *  index ( mname , '  ') 
mname ( nn : nn+4 )  «  '.pll' 

open(unit  =  11,  file  =  mname, status  -  'unknown',  iostat-ios) 
rewind( 11 ) 

DO  1000  1=1,  NUDF 

C - Read  the  data  on  the  first  record  for  this  user-defined  force 

C  element. 


R 


ft 


'5 

*  * 


?  • 


$ 


s 


READ  (  HLDINP ( PNTR ) ,  15,  ERR=100,  IOSTAT-IOSTAT  )  udfnam 
15  FORMAT  ( BN, 10X, A20 ) 

C - If  there  was  a  READ  error,  report  it  and  return  to  the  calling 

C  routine. 

100  CONTINUE 

if  (  iostat  .ne.  0  )  then 

write  (  infof,  115  )  i,  udfnam 

115  format  (/, 'Error  in  reading  data  for  a  user-defined  ', 

&  'force  element  #  ',12,'  ',A,/) 

ERRCOD  =  31032 
RETURN 
ENDIF 

C - Read  the  data  on  the  second  record  for  this  user-defined  force 

C  element. 


C  1  1  1  I  I  I  I  1  !  1  I  1  1  REMEMBER  TO  INCREMENT  PNTR  POINTER  FOR  THE  NEXT  RECORD  1  1 


C  ***  READING  OF  ADDITIONAL  DATA 


1 

1 

1 

125 


read  (  hldinpt pntr  +  1 )  ,  125,  err  =  200, 
( udf ( j ,  i  )  ,  j  =  1,  3) 
read  (  hldinp( pntr  +  2 )  ,  125,  err  =  200, 
( udf ( j , i ) ,  j  =  4,  6) 
read  (  hldinpt pntr  +  3 )  ,  125,  err  =  200, 
( udf ( j , i ) ,  j  =  7,  9) 
FORMAT  ( BN , 3 ( g20 . 1 4 ) ) 


iostat=iostat  ) 
iostat-iostat  ) 
iostat-iostat  ) 


C - If  there  was  a  READ  error,  report  it  and  return  to  the  calling 


c 


routine  . 


CONTINUE 

if  (  iostat  .ne.  0  )  then 

write  (  infof,  215  )  i,  udfnam 

FORMAT  (/.'Error  in  reading  data  for  a  user-defined  ', 
'force  element 12 , '  A,  /) 

ERRCOD  -  31032 
RETURN 
ENDIF 


END  OF  USER  MODIFICATIONS 

c#  ################################### 


C - Continue  to  read  data  from  however  many  records  exist  for  each 

C  user-defined  force  element.  Following  this,  any  input  data  error 
C  checking  can  be  performed. 

C - Set  the  pointer  to  the  first  record  for  the  next  user-defined 

C  force  element. 

PNTR  -  HLDLNK  ( PNTR ) 

1000  CONTINUE 


200 

215 


RETURN 

END 


BLOCKD. FOR 


& 


V 


This  block  of  code  sets  up  how  data  is  stored  in  DADS  for  its 
model  elements,  it  must  be  modified  by  the  user  to  define  how 
data  is  stored  in  the  input  file  for  his  force  calculation 
routine  FRC32_3D.  The  modification  appears  at  the  end  of  the 
block. 


BLOCKD :  Initialization  for  common  block  data. 


■=»Author  and  Modification  Section* 


Author : 


Chuck  Mead 


Date  written:  December  14,  1984 


Written  on:  PRIME  750,  PRlMOS  Rev.  19.1,  compiled  by  F77  Rev.  19.1 


Copyright  (c)  CADSI  1986 


■Algorithm  Description* 


Purpose  and  use: 

This  ‘routine  is  used  to  initialize  any  common  block  variables. 


Error  conditions: 


Machine  dependencies:  none  known 


BLOCKDATA  BLOCKD 


■Variable  Descriptions* 


C - COMMON  blocks- 

C 

C  Harwel . blk 


LIRN. . 
LICN. . 
LIRN2 . 
LICN2 . 
MREGEl 
ISCAL1 
MREGE2 
ISCAL2 
IDISP1 


IDISP2 
U . 


.array  length  of  IRN  for  position  and  velocity  analysis 

.array  length  of  ICN  for  position  and  velocity  analysis 

.array  length  of  IRN  for  acceleration  analysis 

.array  length  of  ICN  for  acceleration  analysis 

.matrix  regeneration  for  position  and  velocity  analysis 
.rescaling  of  the  matrix  for  position  and  velocity  analysis 
.matrix  regeneration  for  acceleration  analysis 
.rescaling  of  the  matrix  for  acceleration  analysis 
.position  of  free  space  of  ICN  for  position  and  velocity 
analys i s 

.position  of  free  space  of  ICN  for  acceleration  analysis 
.pivoting  value  of  sparce  matrix  recomended  value  0.1  -  0.25 
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INTEGER  LIRN,  LICN,  LIRN2,  LICN2,  MREGEl ,  ISCALI, 
&  MREGE2,  ISCAL2,  IDISP1 ( 2 ) ,  IDISP2(2) 

DOUBLE  PRECISION  U 


COMMON  /HARWEL/  LIRN,  LICN,  LIRN 2 ,  LICN2,  U,  MREGEl,  ISCALI, 
&  MREGE2,  ISCAL2,  IDISP1,  IDISP2 

SAVE  /HARWEL/ 


Integv . blk 


IFAIL...Flag  for  prediction  failure 
NORDER . . . Orde r  of  integration 

IDIVG...Flag  for  integration  step  controller  from  DYNAM2 
IPOSIT. . Position  iteration  counter 
IVELOC . . Veloci ty  analysis  counter 
IACCEL . .Acceleration  analysis  counter 

COMMON  /INTEGV/  I FAI L , NORDER , I DI VG , I POS I T , I VELOC , I ACCEL 
INTEGER  I  FAIL, NORDER, I DIVG, IPOSIT, IVELOC, IACCEL 
SAVE  /INTEGV/ 


§ 

c 

I  ouni  t . 

blk 

c 

c 

INPUTF  . 

. . File 

\V 

,  * 

c 

INFOF.  . 

. . File 

c 

OUTPF.  . 

. .File 

/ 

c 

BINRYF  . 

. .File 

( 

c 

NXTREC . 

. . Next 

INTEGER  INPUTF,  INFOF ,  OUTPF,  BINRYF,  NXTREC 

COMMON  /  IOUNIT  /  INPUTF,  INFOF,  OUTPF,  BINRYF,  NXTREC 


SAVE  /IOUNIT/ 


Ptrs . blk 


NPTRS2 , 


NPTRS3 

MPTRS2 


MPTRS3 . 


.Vector  containing  the  number  of  real  data  required  for  an 
element  of  each  module  in  2D. 

.Same  vector  for  3D  modules. 

.Vector  containing  the  number  of  integer  data  required  for 
an  element  of  each  module  in  2D. 

.Same  vector  for  3D  modules. 


INTEGER  NPTRS  2 (50),  NPTRS3(50),  MPTRS2(50),  MPTRS3 ( 50 ) 
COMMON  /  PTRS  /  NPTRS2,  NPTRS 3 ,  MPTRS2 ,  MPTRS 3 


SAVE  /  PTRS  / 


Steps . blk 


NS . Number  of  steps  taken  with  size  H 

NORND. . .Logical  value  for  too  small  error  tolerance  in  STEP 
PHASE1 .. Logical  value  of  predicton  failure  in  STEP 
ALPHA. .. Integration  parameter 
BETA. ...  Integration  parameter 


C  SIG . Integration  parameter 

C  G . Integration  parameter 

C 

INTEGER  NS 

DOUBLE  PRECISION  ALPHA(12),  BETA (12),  SIG(13),  G(13),  V(12), 
&  W  ( 1 2  ) 

LOGICAL  NORND ,  PHASEl 


COMMON  /  STEPS  /  ALPHA,  BETA,  SIG,  G,  V,  W,  NS,  NORND,  PHASEl 
SAVE  /  STEPS  / 

C 

C  Timing. blk 
C 

C  DOTIME ....  Logical  variable  that  is  .TRUE,  if  the  user  wants  to  call 
C  the  system  timing  routines. 

C 

LOGICAL  DOTIME 
COMMON  /  TIMING  /  DOTIME 
C 

C  Rcdnum.blk 

C 

C  REC2D. . . .Vector  of  the  number  of  records  for  each  2D  module. 

C  REC3D. .. .Vector  of  the  number  of  records  for  each  3D  module. 

C 

INTEGER  REC2D ( 50),  REC3D(50) 

COMMON  /  RCDNUM  /  REC2D ,  REC3D 

SAVE  /  RCDNUM  / 

C 

C  Multip.blk 
C 

COMMON  /MULTIP/  1 1 N ( 4 0 ) , I OU ( 4 0 ) , I LN ( 4 0 ) 

INTEGER  I  IN , IOU , I LN 
SAVE  /  MULTIP  / 

C  Crcrds . blk 
C 

COMMON  /  CRCRDS  /  CREC2D ,  CREC3D 
INTEGER  CREC2D( 29  )  ,  CREC3D(29) 

SAVE  /  CRCRDS  / 


C - Local  variables - 

C 

C  I . Loop  index  . 

C 

INTEGER  I 
C 

C - DATA  statements - 

C 

C - Initialization  of  variables  in  the  HARWEL  common  block. 

DATA  U  /  0.25  / 

C - Initialization  of  variables  in  the  INTEGV  common  block. 

DATA  I  POSIT ,  IVELOC,  IACCEL  /  0,  0,  0  / 

C - Added  data  for  MULTIP.  This  is  control  element  data,  1/17/86. 


non 


do 


DATA  I IN  70,3,1, 1,1, 1,1, 1,1,1, 1,1, 1,1, 1,1. 1,0, 1,1, 
&  0,0,0, 1,0, 0,0, 2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,5/ 

DATA  IOU  /l, 1,1, 1,1, 1,1, 1,1, 1,1. 1,1, 1.1, 1,1, 1,2,0, 
&  1,0, 1,0, 3, 1,0, 1,0, 0,0,0, 0,0, 0,0, 0,0, 0,5/ 

DATA  ILN  /l, 1,1, 2, 2, 2, 2, 1,1, 1,1, 1,1, 1,1, 1,1, 1,3,1, 

&  1,1, 1,1, 5, 3, 2, 3, 2, 0,0, 0,0, 0,0, 0,0, 0,0, 3/ 


C - Define  the  number  of  records  per  control  element,  input  data, 


C - 


DATA 

CREC2D  / 

6, 

2, 

2, 

3, 

3, 

& 

2, 

0, 

2, 

2, 

2, 

& 

2, 

0, 

2, 

0, 

0, 

& 

0, 

0, 

0, 

5. 

3, 

& 

0, 

1, 

5, 

3, 

5, 

& 

2, 

2, 

3 

,2 

DATA 

CREC3D  / 

8, 

2, 

2, 

3, 

3, 

& 

2, 

0, 

2, 

2, 

2, 

& 

2, 

0, 

2, 

0, 

0, 

& 

0, 

0, 

0, 

6, 

5, 

& 

0, 

1, 

6, 

3, 

6, 

& 

2, 

2, 

3 

,2 

Ini tialization 

of 

va  riables 

in 

DATA 

INPUTF , 

INFOF 

,  OUTPF 

DATA 

BINRYF 

/ 


/  100,  102,  101  / 

/  103  / 


-Initialization  of  variables  in  the  PTRS  common  block,  except  for 
array  elements  number  32,  defined  below  for  the  user-defined 
force  element  module. 


DATA 

( NPTRS2 ( I ) 

,1-1,31) 

/ 

o, 

0, 

20, 

6, 

7, 

& 

10, 

14, 

11, 

2, 

19, 

& 

9, 

2, 

11, 

1, 

6, 

& 

10, 

11, 

13, 

10, 

0, 

& 

7, 

9, 

4, 

8, 

6, 

& 

0, 

0, 

0, 

0, 

0, 

& 

0 

DATA 

( NPTRS2 ( I ) 

,1*33,50) 

/ 

16* 

0,  14 

,  6 

/ 

DATA 

( NPTRS  3 ( I ) 

,1-1,31) 

/ 

0, 

0, 

65, 

19, 

21, 

& 

36, 

36, 

37, 

36, 

43, 

& 

26, 

3, 

37, 

1, 

36, 

& 

36, 

37, 

36, 

37, 

0, 

& 

39, 

37, 

36, 

36, 

0, 

& 

0, 

0, 

0, 

0, 

0, 

& 

30 

DATA 

( NPTRS  3 ( I ) 

,1-33,50) 

/ 

11* 

0,  62 

,  4 

*0, 

40, 

DATA 

( MPTRS2 ( I ) 

,1-1,31) 

/ 

o. 

Or 

10, 

3, 

2, 

& 

5, 

5, 

3, 

3, 

7, 

& 

6  , 

3, 

3, 

1. 

4, 

& 

4, 

4, 

3, 

5, 

0, 

& 

3, 

3, 

3, 

4, 

3, 

& 

0, 

0, 

0, 

0, 

21, 

& 

0 

DATA 

( MPTRS2 ( I ) 

,1-33,50) 

/ 

16* 

0,  6, 

4  / 

DATA 

( MPTRS3 ( I ) 

,1-1,31) 

/ 

0, 

0, 

10, 

3, 

2, 

& 

5, 

5, 

5, 

5, 

7, 

& 

7, 

3, 

5, 

1, 

5, 

& 

5, 

6, 

3, 

3,  0, 

& 

4, 

3, 

3, 

3,  0, 

& 

0, 

0, 

0, 

0,  21, 

& 

DATA  ( MPTRS3 ( I ) , 1-33 , 50 ) 

7 

/  11*0 

,  6, 

4*0, 

8, 

-Initialization  of  variables  in  the  RCDNUM  common  block,  except  for 
array  elements  number  32,  defined  below  for  the  user-defined 
force  element  module. 


DATA 

( REC2D ( 1  )  , 1-1 , 31  ) 

/ 

0, 

4  , 

5, 

2, 

& 

3, 

3, 

3, 

4, 

& 

2  , 

5, 

3, 

2, 

& 

4  , 

0, 

3, 

2, 

& 

4, 

4, 

3, 

0, 

& 

2, 

3, 

2, 

3, 

& 

3, 

2, 

1 , 

1  , 

& 

1, 

0, 

0 

/ 

DATA 

( REC2D ( I ) , 1-33 , 50 ) 

/ 

16 

*0, 

5, 

3 

/ 

DATA 

( REC  3  D ( I  )  ,  1-1 , 31  ) 

/ 

0, 

4  , 

9  , 

4  , 

& 

5, 

7, 

7  , 

8, 

& 

7, 

9, 

3, 

2, 

& 

8, 

0, 

7, 

7, 

& 

8, 

7  , 

8  , 

0, 

& 

8  , 

8, 

7, 

7  , 

& 

0, 

2, 

1, 

1, 

& 

1  , 

0, 

6 

/ 

DATA 

( REC  3  D ( I ) , 1-3  3 , 50  ) 

/ 

2, 

10* 

O 
» — ♦ 

o 

4 

*0, 

Initialization  of  variables 

in 

the 

STEPS 

common 

DATA 

G ( 1 ) ,  G ( 2 !  /  1 .0, 

0.5 

/ 

DATA 

S I G  (  1  )  /  1 . 0 

/ 

■Initialization  of  variables  in  the  TIMING  common  block. 
DATA  DOTIME  /  .FALSE.  / 


-Initialization  of  variables  required  for  implementation  of  a  user- 
defined  force  element  module. 


NPTRS 2 (32)  defines  the  number  of  real  (double  precision)  data 
stored  for  each  user-defined  force  element  in  DADS-2D . 

DATA  NPTRS  2 ( 32 )  /  0  / 

NPTRS 3(32)  defines  the  number  of  real  (double  precision)  data 
stored  for  each  user-defined  force  element  in  DADS-3D . 


#######################11######## 
USER  MODIFICATION 
The  value  of  NPTR32(32)  for  this  problem  has  been  changed  from 


!,t*!*l*!*!  .VJVn 


DATA  NPTRS3 ( 32 )  /  9  / 

C - MPTRS 2 (32)  defines  the  number  of  integer  data  stored  for  each  user- 

C  defined  force  element  in  DADS-2D. 

DATA  MPTRS2 ( 32  )  /  0  / 

C - MPTRS  3 ( 32 )  defines  the  number  of  integer  data  stored  for  each  user- 

C  defined  force  element  in  DADS- 3D . 

DATA  MPTRS 3 (  32 )  /  0  / 

C - REC2D(32)  defines  the  number  of  records  in  the  input  file  for  each 

C  user-defined  force  element  in  DADS-2D. 

DATA  REC2D ( 32 )  /  1  / 

C - REC 3D ( 3  2 )  defines  the  number  of  records  in  the  input  file  for  each 

C  user-defined  force  element  in  DADS-3D. 

c  i  i  m  i  >  i  i  i  ■  i  i  i  t  i  u  i  «  u  «  i  n  «  u  i  i  i  ii  i  it 

USER  MODIFICATION 

The  value  of  REC3D(32)  for  this  problem  has  been  changed  from 
a  1  to  a  4  . 

DATA  REC  3  D (  32 )  /  4  / 


RPT32  3D. FOR 


'5 


IWSSV 


RWWl 
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This  routine  allows  the  user  to  output  additional  data  that  is 
calculated  by  their  force  calculation  routine  FRC32_3D.  Most  of 
this  routine  is  a  vendor  supplied  template,  the  user's 
modifications  appear  at  its  end  and  are  designated  as  such. 


RPT 3 2 :  Reports  user-defined  force  data  at  each  output  time  step  (3D) 

‘-----■•“---Autho r  and  Modification  Sect i on*»»*»”»*»«»«»»— »—»»»«««' 

Author:  Chuck  Mead 

Date  written:  December  5,  1985 
Modi f i cat ions : 

Copyright  (c)  CADSI  1985 


■•Algorithm  Description* 


Purpose  and  use: 

This  subroutine  reports  any  data  associated  with  the  user- 
defined  force  elements  in  the  r.odel.  Any  of  the  various  input 
arguments  for  this  subroutine  can  be  reported  or  used  to  calculate 
reported  values.  The  output  data  can  be  written  to  the  standard 
DADS  ASCII  output  file  (FORTRAN  file  unit  OUTPF)  or  to  some  other 
user-defined  output  file(s).  The  function  GETNAM  can  be  used  to 
convert  an  element  number  into  the  corresponding  element  name  for 
writing  out  to  the  output  file(s).  This  function  can  be  used  to 
retrieve  the  name  of  each  user-defined  force  element  processed  if 
the  first  input  record  of  each  user-defined  force  element  contains 
such  a  name  in  columns  11  through  30. 

The  various  input  arguments  are  described  and  declared  in  three 
separate  groups:  read-only  arguments,  arguments  read  and  modi¬ 
fied,  and  return-only  arguments. 

Read-only  arguments  are  passed  in  with  meaningful  values  which 
should  not  be  changed  within  this  subroutine.  This  constitutes 
the  largest  group  of  arguments  passed.  Consult  the  documentation 
on  user-defined  force  elements  for  a  description  of  the  structure 
and  contents  of  the  various  variables  and  arrays  passed  in. 

Arguments  read  and  modified  are  also  passed  in  with  meaningful 
values  which  may  be  changed  within  this  subroutine.  There  are  no 
read-and-modi fy  arguments  passed  to  RPT32 . 

The  final  group  of  arguments,  return-only  arguments,  are  passed 
in  with  unknown  values  and  should  therefore  be  assigned  some  value 
within  this  subroutine  before  they  are  used  in  any  other  computa¬ 
tion.  The  only  variable  in  this  group  is  ERRCOD.  Since  this 
variable's  value  is  used  in  the  calling  routine  to  determine  the 
success  of  the  computations  performed  in  RPT32 ,  it  is  suggested 
that  ERRCOD  be  initialized  to  zero  before  beginning  any  other 
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computations . 

Error  conditions: 

Error  conditions  are  reported  via  the  variable  ERRCOD .  A  zero 
value  indicates  no  error;  a  positive  value  indicates  some  fatal 
error,  i.e.,  files  will  be  closed  and  execution  aborted  upon 
return  to  the  calling  routine;  negative  values  are  reserved  for 
future  use. 

Machine  dependencies:  none  known 


SUBROUTINE  RPT 3 2  (  UDF,  IUDF,  NUDF,  NPTRS ,  MPTRS,  TIME, 

&  Q,  QD,  QDD ,  XQ,  XQD ,  XQDD ,  NXCRD ,  FRC , 

&  XFRC ,  RB,  NRB ,  NPRB ,  NTYPES,  NELEMS, 

&  INFOF,  OUTPF,  BODMOD,  ERRCOD  ) 


c  P 

C - Read-only  arguments - 

n 

C  UDF . Array  of  real  (double  precision)  data  for  all  user- 

C  defined  force  elements  in  the  model. 

C  IUDF . Array  of  integer  data  for  all  user-defined  force  elements 

C  in  the  model. 

C  NUDF . Number  of  user-defined  force  elements  in  the  input  data. 

C  (This  variable  has  a  minimum  value  of  '1  for  array  dec- 

C  laration  purposes.) 

C  NPTRS .... Numbe r  of  real  (double  precision)  data  stored  for  each 
C  user-defined  force  element.  (This  variable  has  a  minimum 

C  value  of  1  for  array  declaration  purposes.) 

C  MPTRS .... Numbe r  of  integer  data  stored  for  each  user-defined  force 
C  element.  (This  variable  has  a  minimum  value  of  1  for 

C  array  declaration  purposes.) 

C  TIME . Current  time  in  the  model  simulation. 

C  Q . Array  of  displacement  values  for  all  rigid  body  coor- 

C  dinates  in  the  model. 

C  QD . Array  of  velocities  of  all  rigid  body  coordinates. 

C  QDD . Array  of  accelerations  of  all  rigid  body  coordinates. 

C  XQ . Array  of  displacement  values  for  all  extra  generalized 

C  coordinates  in  the  model. 

C  XQD . Array  of  velocities  of  all  extra  generalized  coordinates. 

C  XQDD . Array  of  accelerations  of  all  extra  generalized  coor- 

C  dinates. 

C  NXCRD .... Numbe r  of  extra  generalized  coordinates  in  the  model. 

C  (This  variable  has  a  minimum  value  of  1  for  array  dec- 

C  laration  purposes.) 

C  FRC . Array  of  forces  for  each  rigid  body  coordinate  in  the 

C  model. 

C  XFRC . Array  of  forces  for  each  extra  generalized  coordinate  in 

C  the  model. 

C  RB . Array  of  real  (double  precision)  data  for  each  rigid  body 

C  in  the  model . 

C  NRB . Number  of  rigid  bodies  in  the  model. 

C  NPRB . Number  of  real  (double  precision)  data  stored  for  each 

C  rigid  body  in  the  model. 

C  NTYPES ... Numbe r  of  different  element  types  available  in  DADS. 

C  NELEMS ... Numbe r  of  occurrances  of  each  element  type  in  the  model. 

C  INFOF ....  FORTRAN  file  unit  number  for  the  DADS  information  file. 


nnonn  onoonnnnoonooo  onnonnnn  nnonnono 


Error,  warning,  or  other  informational  messages  can  be 
written  to  this  file  if  so  desired. 

OUTPF ....  FORTRAN  file  unit  number  for  the  standard  DADS  ASCII 
output  file. 

BODMOD ...  Element  type  number  of  the  rigid  body  module.  Can  be 
used  to  map  body  numbers  to  body  names  via  the  GETNAM 
function . 

INTEGER  NUDF ,  NPTRS ,  MPTRS ,  IUDF ( MPTRS , NUDF ) ,  NXCRD ,  NRB ,  NPRB , 

&  NTYPES ,  NELEMS ( NTYPES ) ,  INFOF ,  OUTPF,  BODMOD 

DOUBLE  PRECISION  UDF ( NPTRS , NUDF ) ,  TIME,  Q ( 7 , NRB ) ,  QD ( 7 , NRB ) , 

&  QDD ( 7 , NRB )  ,  XQ(NXCRD),  XQD( NXCRD ) ,  XQDD ( NXCRD  )  , 

&  FRC ( 7 , NRB ) ,  XFRC ( NXCRD ) ,  RB ( NPRB , NRB ) 

- Arguments  read  and  modified - 

none 

- Return-only  arguments - - - 

ERRCOD. .. Error  condition  flag. 

INTEGER  ERRCOD 

- COMMON  blocks - 

none 

- Local  variables - 

none 

- Functions  and  subroutines - 

GETNAM ...  Given  the  element  type  number  (i.e.,  module  number)  and 
the  element  number,  this  function  returns  the  element 
name.  A  blank  function  value  means  that  either  or  both 
of  the  input  values  are  invalid. 

EXTERNAL  GETNAM 

CHARACTER*  2  0  GETNAM 

- DATA  statements - 

none 


USER  MODI  FACATION 
common  /trie/  friction,  velocity,  accel,  ext_torque 
double  precision  friction,  velocity,  accel,  ext_torque 
save  /trie/ 


■Process  Block' 


C - First  zero  out  the  error  condition  flag  to  indicate  that  no  errors 

C  have  occurred  yet. 

ERRCOD  -  0 


write  (11,  ' ( 5 ( 2 x ,  pel2.5) ) ' ) 

1  time,  friction,  velocity,  accel,  ext_torque 

C#  ################################## 


RETURN 

END 


FART  NUMBER: 


12585721-5724.  Walking  Beam  (Wheel;  Actuators 
BE SEP IF  TION:  WALKING  BEAU  ACTUATORS 


All  functional  and  dimensional  requirements  have  been  determined 
bv  FMC  for  each  of  the  four  wheel  (walking  beam)  actuators  and 
based  on  this  information,  V  or  k  has,  completed  preliminary 
drawings  (TDP,  Dwgs.  12585721-12585724'. 

It  should  be  noted  that  the  rear  actuators  are  designed  to  not 
on 1 v  lift  the  rear  of  the  LTHD  off  the  ground  during  routine 
disci  a cement  but  can  lift  a  "stuck"  spade  out  of  muddv  ground  as 

Also  note  that  the  wheel  actuators  are  designed  to  act  as  "shock 
absorbers''  in  the  LTHD '  s  travel  mode  and  will  allow  the  walking 
beams  to  walk-  freely.  If  the  LTHD  encounters  a  large  obstacle  at 
a  high  speed,  however,  the  limited  flow  rate  of  the  cylinders  mav 
cause  the  sheet  absorbing  system  tc  be  unnecessarily  stiff  and 
decrease  towing  stability.  Possible-  solutions  to  this  problem 
ar e  1 i si ed  in  th i s  sec 1 1  on . 
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SHEARING  RESISTANCE  AND  STRENGTH 


Shearing  forct 


.mrnmtm. 

— — 

for  Direct  Shear  Test 

of  failure  of  the  sample  along  the 
-_,uired  to  produce  continuous  move- 


length  of  the  specimen  under  the 
‘  iS  test. 


i  maximum  applied  shearing  force 
't  the  shear  plane  gives  the 
™Si  .  failure  of  the  sample.  These 
igram  like  that  in  Fig.  18-6.  The 
.^*nt  on  the  graph  representing  the 
,,th.  In  order  to  obtain  sufficient 
Itional  tests  must  be  performed 
^ates  of  the  first.  The  procedure 
n  the  first,  except  that  a  different 

VVrrp npri  .  1 1  1  • 


-►generalization,  the  plotted  points 
ure  of  the  various  specimens  will 
.V there  will  be  deviations  due  to* 
Uniformity  of  the  specimens.  Also 
e  clay  soils  in  the  undrained  test, 
;en  normal  stress  and  shearing 
•-■dally  at  low  values  of  normal 
"practice  to  draw  the  best  straight 
the  Coulomb  graph  and  equation 
-“.iction  angle  and  cohesion  can  be 


t  pa  riving  data  were  obtained  in  tests 
.  glacial  drift.  Determine  the  cohesion 


e- 180  psf 

0-32* 


180  +N  tan  32*! 


Normal  Pressure,  in  Pounds  per  sq  ft 
Tig.  18-8.  Shear  Diagram  for  Example  18-1 


Specimen 

No. 


Failure  Stresses,  in  Pounds  per  sq  ft 


Normal  Load 


Shearing  Force 


Solution:  The  data  are  plotted  on  the  shear  diagram  in  Fig.  18-8.  A  straight 
line  through  the  plotted  points  makes  an  angle  of  32°  with  the  horizontal,  and 
the  intercept  on  the  vertical  axis  represents  180  psf.  Therefore,  the  shearing 
strength  properties  of  the  soil  arc  designated  as  <f=  180  psf  and  £  =  32°. 
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Fig.  18-9.  Schematic  Diagram  of  Triaxial  Compres¬ 
sion  Device 
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